exn544.m

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

%《MATLAB及其在理工课程中的应用指南》(十一五规划版）第五章例5-4-4程序
% 数值方法求回转体体积
% 西安电子科技大学出版社出版  陈怀琛编著 2007年7月
%
clear,clf
dx=input('dx= ') ;x=0:dx:pi; 
g=x.*sin(x).^2; plot(x,g),hold on
fill([1.5,1.6,1.6,1.5,1.5],[0,0,[1.6,1.5].*sin([1.6,1.5]).^2,0],'r')
pause
v=trapz(2*pi*x.*g*dx)       % 绕y轴回转所得回转体体积
pause
v1=trapz(pi*g.^2*dx)       % 绕x轴回转所得回转体体积
pause
syms x, g=x*sin(x)^2;
v1t= int(pi*g^2,0,pi),double(v1t)      % 绕x轴回转所得回转体体积
hold off