mygauss.m

自己编写的matlab程序 数值分析里的程序

clc
a=1;b=6;
x0=[-sqrt(15)/5 0 sqrt(15)/5];
% (b-a)/(1-(-1))=(b-x)/(1-x0)
x=b-(b-a)*(1-x0)/2;
ygauss=(5/9*fun(x(1))+8/9*fun(x(2))+5/9*fun(x(3)))*(b-a)/2
y=quadl('fun(x)',a,b)
h=(b-a)/2;
ysimpson=h/3*(fun(a)+4*(fun(a/2+b/2))+fun(b))
if abs(ygauss-y)<abs(ysimpson-y)
    disp('Gauss Algorithm is better than Simson Algorithm');
elseif abs(ygauss-y)<abs(ysimpson-y)
    disp('Simson Algorithm is better than Gauss Algorithm');
else
    disp('Same')
end