中quad.m函数,我把它简化了, [Q] = (,a,b)f = ();h = 0.*(b-a);x = [a a+h a+2*h (a+b)/2 b-2*h b-h b];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 = [a a+h a+2*h (a+b)/2 b-2*h b-h b];

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);

% 's rule.

Q1 = (h/6)*(fa + 4*fc + fb);

% Five 's rule.

Q2 = (h/12)*(fa + 4*fd + 2*fc + 4*fe + fb);

% One step of .

Q = Q2 + (Q2 - Q1)/15;

if abs(Q2 - Q)