exn545.m

来自「MATLAB及在电子信息课程中的应用.感觉很有用,本来也不花钱,有想用的来拿吧」· M 代码 · 共 18 行

%《MATLAB及其在理工课程中的应用指南》(十一五规划版）第五章例5-4-5程序
% 二重积分
% 西安电子科技大学出版社出版  陈怀琛编著 2007年7月
%
clear, close all
format compact
fill([0,1,1,0],[0,0,1,0],'y'),hold
fill([0.55,0.6,0.6,0.55,0.55],[0,0,0.6,0.55,0],'r')
dx=input('dx=  ');dy=dx;
x=0:dx:1;lx=length(x);
for k=1:lx
   x1=(k-1)*dx;
   y1=0:dy:x1;
   f=x1.^2+y1.^2;
   s1(k)=trapz(f)*dy;
end
s=trapz(s1)*dx