exn546.m

来自「基天子matlab的双缝干涉从单缝(线光源)发出的单色光射到与之平行的双缝上」· M 代码 · 共 19 行

%《MATLAB及其在理工课程中的应用指南》(十一五规划版）第五章例5-4-6程序
% 三重积分
% 西安电子科技大学出版社出版  陈怀琛编著 2007年7月
%
clear, close all
dx=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;	    %z1数组
      f=x1.*y1.^2.*z1.^3;   %f(z1)
      s1(j)=trapz(f)*dz;    %沿z1积分
   end
   s2(k)=trapz(s1)*dy;		% 沿y1积分
end,
s=trapz(s2)*dx		% 沿x1积分