ep2_p0.m

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

%  Ep2_p0: << 实验二 >>  教堂顶部面积的计算
%  Designed by FGH

%  球面面积公式为：S = 4*PI*R^2

%  故哈桑计算的教堂顶部面积为

  s1 = 2*pi*30^2;
  sprintf('     %7.2f',s1)
  
  pause

%  加上损耗后的面积为

  s2 = (1+0.015)*s1;
  sprintf('     %7.2f',s2)
  
  pause

%   从而盈余金箔(平方米)

  s3 = 5750 - s2

  pause
  
%   由算式 (2.13) 直接计算二重积分得
  
  s4 = dblquad('ep2_f0',0,1,0,2*pi);
  sprintf('     %7.2f\n',s4)
  
  pause
  
%  加上损耗后的面积为
  
  s5 = (1+0.015)*s4;
  sprintf('     %7.2f\n',s5)
  
  pause
  
%  实际亏空金箔(平方米)
  
  5750-s5