disp('Metodo de Simpson');
a=input('valor de a:');
b=input('valor de b:');
n=input('valor de n:');
Ar=0;
h=(b-a)/n;
for i=0:n-1
xi=a+i*h;
xii=a+(1+i)*h;
xm=(xi+xii)/2;
S=1/(1-xi)+4*(1/(1-xm))+1/(1-xii);
if (S<0)
S=(-1)*S;
else
S=S*1;
end
Ar=Ar+ h/6*(S);
end
fprintf('el valor del area es %.8f',Ar);
No hay comentarios:
Publicar un comentario