ex544.m

来自「《MATLAB及其在理工课程中的应用指南》书上的程序集」· M 代码 · 共 21 行

%《MATLAB及其在理工课程中的应用指南》第五章例5－4－4程序
% 三重积分
% 西安电子科技大学出版社出版  陈怀琛编著 1999年10月初版，2004年11月第二版
%
clear,format compact
dx=input('dx=(书上为0.1及0.01) ');dy=dx;dz=dx;
x=0:dx:1;
for k=1:length(x)
    x1=(k-1)*dx;
    y=0:dy:x1;
    for j=1:length(y)
        y1=(j-1)*dy;
        z1=0:dz:x1*y1;
        f=x1.*y1.^2.*z1.^3;
        s1(j)=trapz(f)*dz;
    end
    s2(k)=trapz(s1)*dy;
end
s=trapz(s2)*dx
%dx=0.1;s=0.0030;
%dx=0.01;s=0.0027