%《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