romberg2

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

deff('y=f(x)','y=(x^2)*(exp(x^2))')
a=0
b=2
h=b-a
R=zeros(6,6)
R(1,1) = h*(f(a) + f(b))/2
for k = 2:6
 R(k,1) = R(k-1,1)/2
 for i=1:2^(k-2)
 R(k,1) = R(k,1) + h*f(a + (2*i - 1)*(h/(2^(k-1))))/(2^(k-1))
 end
end
for j=2:6
 for k=j:6
 R(k,j) = R(k,j-1) + (R(k,j-1) - R(k-1,j-1))/(4^(j-1) - 1)
 end
end
disp(R)