Este documento presenta tres funciones (trapecio, simpson13, simpson38) para aproximar el cálculo numérico de integrales definidas utilizando los métodos del trapecio, Simpson 1/3 y Simpson 3/8 respectivamente. Cada función solicita la entrada de la función, los límites inferior y superior de integración, el número de trapecios o subdivisiones a considerar, y devuelve el valor aproximado de la integral y genera una gráfica de la función.
0 calificaciones0% encontró este documento útil (0 votos)
400 vistas3 páginas
Este documento presenta tres funciones (trapecio, simpson13, simpson38) para aproximar el cálculo numérico de integrales definidas utilizando los métodos del trapecio, Simpson 1/3 y Simpson 3/8 respectivamente. Cada función solicita la entrada de la función, los límites inferior y superior de integración, el número de trapecios o subdivisiones a considerar, y devuelve el valor aproximado de la integral y genera una gráfica de la función.
Este documento presenta tres funciones (trapecio, simpson13, simpson38) para aproximar el cálculo numérico de integrales definidas utilizando los métodos del trapecio, Simpson 1/3 y Simpson 3/8 respectivamente. Cada función solicita la entrada de la función, los límites inferior y superior de integración, el número de trapecios o subdivisiones a considerar, y devuelve el valor aproximado de la integral y genera una gráfica de la función.
Este documento presenta tres funciones (trapecio, simpson13, simpson38) para aproximar el cálculo numérico de integrales definidas utilizando los métodos del trapecio, Simpson 1/3 y Simpson 3/8 respectivamente. Cada función solicita la entrada de la función, los límites inferior y superior de integración, el número de trapecios o subdivisiones a considerar, y devuelve el valor aproximado de la integral y genera una gráfica de la función.
Descargue como DOCX, PDF, TXT o lea en línea desde Scribd
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