中quad.m函数,我把它简化了,[Q]=(,a,b)f=();h=0.*(b-a);x=[aa+ha+2*h(a+b)/2b-2*hb-hb];y=f(x);[Q(1)]=...(f,x(1),x(3),y(1),y(2),y(3));[Q(2)]=...(f,x(3),x(5),y(3),y(4),y
2016-03-30
中quad.m函数,我把它简化了,
[Q]=(,a,b)
f=();
h=0.*(b-a);
x=[aa+ha+2*h(a+b)/2b-2*hb-hb];
y=f(x);
[Q(1)]=...
(f,x(1),x(3),y(1),y(2),y(3));
[Q(2)]=...
(f,x(3),x(5),y(3),y(4),y(5));
[Q(3)]=...
(f,x(5),x(7),y(5),y(6),y(7));
Q=sum(Q);
[Q]=(f,a,b,fa,fc,fb)
tol=1.e-6;
h=b-a;
c=(a+b)/2;
x=[(a+c)/2,(c+b)/2];
y=f(x);
fd=y(1);
fe=y(2);
%'srule.
Q1=(h/6)*(fa+4*fc+fb);
%Five'srule.
Q2=(h/12)*(fa+4*fd+2*fc+4*fe+fb);
%Onestepof.
Q=Q2+(Q2-Q1)/15;
ifabs(Q2-Q)
【下载地址】magnet推荐使用utorrent、BitComet等bt客户端下载
为了防止迅雷扫版权风险,本片迅雷链接,在2个小时后自动放出!