Nothing Special   »   [go: up one dir, main page]

La Regla Del Trapecio en Matlab

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 2

Practica de la Regla del Trapecio (realizado en clase)

CODIGO:
clc

clear

close

disp('REGLA DEL TRAPECIO')

syms x

f=input('Ingrese la Función fx: ');

a=input('Ingrese el límite inferior del intervalo a: ');

b=input('Ingrese el límite superior del intervalo b: ');

n=input('Ingrese el número de intervalos n: ');

grid on

ezplot(f,(-2:-1))

if n==1

disp('Aplicación Simple de la Regla del Trapecio')

%Evaluar el valor real

x=b;

fb=eval(f);

x=a;

fa=eval(f);

Ireal= (b-a)*(fa+fb)/2;

%Evaluar el valor aproximado

I=int(f);

x=a;

fa=eval(I);

x=b;

fb=eval(I);

Iaprox=(fb-fa);

% Calcular el error real

Ereal = abs(Ireal-Iaprox);
else if n>1

disp('Aplicación Múltiple de la Regla del Trapecio')

% Evaluar el valor real

x=b;

fb=eval(f);

x=a;

fa=eval(f);

Ireal= (b-a)*(fa+fb)/2;

%Evaluar el valor aproximado

f=inline(f);

h=(b-a)/n;

Iaprox=fa+fb;

for i=1:n-1

x=a+i*h;

Iaprox=Iaprox+2*f(x);

end

Iaprox=(h/2)*Iaprox; a=0;

Ereal=abs(Ireal-Iaprox);

end

end

fprintf('El Valor Real de la integral es: %9.5f\n',Ireal)

fprintf('El Valor Aproximado de la Integral es: %9.5f\n',Iaprox)

fprintf('El Error Real es: %9.5f\n',Ereal)

También podría gustarte