ep2_p2.m

来自「此代码展示了Euler计算方法与差分方法」· M 代码 · 共 28 行

%  Ep2_p2: << 实验二 >>  数值积分    Simpson 算式（2.10） 
%  Designed by FGH

% m= 4;
a= 30.6;
b= 29.6;
R= 30;
h= 1/(2*m);
k= 2*pi/(2*m);
e= 0 : k : 2*pi;
t= (0 : h : 1)';

%     算式(2.13)
f= sqrt(t.^2*ones(size(t))' + R^2*(1-t.^2)*((cos(e)/a).^2+(sin(e)/b).^2));

clear Iij;
for j= 2:2:2*m
   for i= 2:2:2*m
%     算式(2.10)
      Iij(i,j)= k*h/9*( f(i-1,j-1)+f(i+1,j-1)+f(i-1,j+1)+f(i+1,j+1)...
              + 4*(f(i,j-1)+f(i-1,j)+f(i+1,j)+f(i,j+1))...
              + 16*f(i,j) );
   end
end
I= sum(sum(Iij));
S= a*b*I;
sprintf(' m = %d ,      S = %7.2f\n',m,S)