exn541.m

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

%《MATLAB及其在理工课程中的应用指南》(十一五规划版）第五章例5-4-1程序
% 求曲线所围面积
% 西安电子科技大学出版社出版  陈怀琛编著 2007年7月
%
clear, close all
for dx=[2,1,0.5,0.1]   						% 设不同步长
    x=0:.1:10;y=-x.*x+115;              % 取较密的函数样本
    plot(x,y),hold on               	% 画出被积曲线，并保持
    x1=0:dx:10;y1=-x1.*x1+115; 				% 求取样点上的y1
    n=length(x1);s=sum(y1(1:n-1))*dx; 	% 用欧拉法求积分，注意末尾去掉一点
    q=trapz(y1)*dx; 						% 用梯形法求积分
    stairs(x1,y1),plot(x1,y1)			% 画出欧拉法及梯形法的积分区域
    [dx,s,q],pause,hold off		
end