计算方法实验截图

function t=agui_trapz(fname,a,b,n)

h=(b-a)/n;

fa=feval(fname,a);

fb=feval(fname,b);

f=feval(fname,a+h:h:b-h+0.001*h); t=h*(0.5*(fa+fb)+sum(f));

format long

t=agui_trapz(inline('x.*exp(x)'),1,2,8)

function s=agui_simpson(fname,a,b,n)

h=(b-a)/n;

fa=feval(fname,a);

fb=feval(fname,b);

s=fa-fb;

x=a;

for i=1:n

x=x+h/2;s=s+4*feval(fname,x);

x=x+h/2;s=s+2*feval(fname,x); end

s=s*h/6;

s=agui_simpson(inline('x.*exp(x)'),1,2,4)

function r=agui_rbg(fname,a,b)

e=0.5e-7;

i=1;j=1;h=b-a;

T(i,1)=h/2*(feval(fname,a)+feval(fname,b));

T(i+1,1)=T(i,1)/2+sum(feval(fname,a+h/2:h:b-h/2+0.001*h))*h/2;

T(i+1,j+1)=4^j*T(i+1,j)/(4^j-1)-T(i,j)/(4^j-1);

while abs(T(i+1,i+1)-T(i,i))>e

i=i+1;h=h/2;

T(i+1,1)=T(i,1)/2+sum(feval(fname,a+h/2:h:b-h/2+0.001*h))*h/2;

for j=1:i;

T(i+1,j+1)=4^j*T(i+1,j)/(4^j-1)-T(i,j)/(4^j-1);

end

end

T

r=T(i+1,j+1);

agui_rbg(inline('x.*exp(x)'),1,2)