DEL RÍO ALCÁNTARA XIMENA 3AM1
clc
syms x
f=input('Intresa la funcion a integrar: ');
b=input('Limite superior: ');
a=input('Limite inferior: ');
k=input('Aproximar cuando k= ');
fx=inline(f);
real=int(f,a,b)
fa=fx(a);
fb=fx(b);
if k==1
I11=(b-a)*((fa+fb)/2);
rombergfinal=I11;
'Valor Calculado por Método de Romberg'
I11
else
end
if k==2
n=2
h=(b-a)/n;
sint=fx(h);
I11=(b-a)*((fa+fb)/2);
I21=(b-a)*((fa+(2*sint)+fb)/(2*n));
I12=((4*I21)-I11)/3;
rombergfinal=I12;
'Valor Calculado por Método de Romberg'
I12
else
end
if k==3
n=4
ha=(b-a)/2;
h=(b-a)/n;
x1=a+h;
x2=a+2*h;
x3=a+3*h;
fx1=fx(x1);
fx2=fx(x2);
fx3=fx(x3);
sint=fx(ha);
si=fx1+fx2+fx3;
I11=(b-a)*((fa+fb)/2);
I21=(b-a)*((fa+(2*sint)+fb)/(2*2));
I31=(b-a)*((fa+(2*si)+fb)/(2*n));
I12=((4*I21)-I11)/3;
I22=((4*I31)-I21)/3;
I13=((16*I22)-I12)/15;
rombergfinal=I13;
'Valor Calculado por Método de Romberg'
I13
