Codigo en Matlab de Integracion Numerica

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

function trapecio

f=input('ingrese la funcion a integrar f(x)= ','s');


a=input('ingrese el limite inferior:');
b=input('ingrese el limite superior:');
fprintf('ingrese numero de trapecios a \n');
n=input('cosiderar en la integracion:');
xmin=a-1;xmax=b+1;
h=(b-a)/n;
x=a:h:b;
fx=eval(f);y=abs(fx);
A=y(1)+y(n+1);
B=2*sum(y(2:n));
integral=(h/2)*(A+B);
fprintf('el area es:%10.9f\n',integral);
%grafica
xp=xmin:0.2:xmax;
x=xp;
yp=eval(f);
plot(xp,yp,'g');
hold on
x=a:0.05:b;
y=eval(f);
bar(x,y,'r');
grid
function simpson13
f=input('ingrese la funcion a integrar f(x)= ','s');
a=input('ingrese el limite inferior:');
b=input('ingrese el limite superior:');
fprintf('ingrese numero de trapecios a \n');
n=input('cosiderar en la integracion:');
n=2*n;xmin=a-1;xmax=b+1;
h=(b-a)/n;
x=a:h:b;
fx=eval(f);y=abs(fx);
suma1=y(1)+y(n+1);
suma2=4*sum(y(2:2:n));
suma3=2*sum(y(3:2:n-1));
suma=suma1+suma2+suma3;
integral=(h/3)*suma;
fprintf('el area es:%10.9f\n',integral);
%grafica
xp=xmin:0.2:xmax;
x=xp;
yp=eval(f);
plot(xp,yp,'g');
hold on
x=a:0.05:b;
y=eval(f);
bar(x,y,'r');
grid
function simpson38
f=input('ingrese la funcion a integrar f(x)= ','s');
a=input('ingrese el limite inferior:');
b=input('ingrese el limite superior:');
fprintf('ingrese numero de trapecios a \n');
n=input('cosiderar en la integracion:');
n=3*n;xmin=a-1;xmax=b+1;
h=(b-a)/n;
x=a:h:b;
fx=eval(f);y=abs(fx);
suma1=y(1)+y(n+1);
suma2=3*sum(y(2:3:n-1));
suma3=3*sum(y(3:3:n));
suma4=2*sum(y(4:3:n-2));
suma=suma1+suma2+suma3+suma4;
integral=(3/8)*h*suma;
fprintf('el area es:%10.9f\n',integral);
%grafica
xp=xmin:0.2:xmax;
x=xp;
yp=eval(f);
plot(xp,yp,'g');
hold on
x=a:0.05:b;
y=eval(f);
bar(x,y,'r');
grid

También podría gustarte