code.m

1% Evaluation of Volume of solid of revolution
2clear
3clc
4syms x
5f(x)=sqrt(x); % Given function
6yr=1; % Axis of revolution y=yr
7I=[0,4]; % Interval of integration
8a=I(1);b=I(2);
9vol=pi*int((f(x)-yr)^2,a,b);
10disp('Volume of solid of revolution is: ');
11disp(vol); % Visualization of solid of revolution
12fx=matlabFunction(f);
13xv = linspace(a,b,101); % Creates 101 points from a to b
14[X,Y,Z] = cylinder(fx(xv)-yr);
15Z = a+Z.*(b-a); % Extending the default unit height of the
16% cylinder profile to the interval of integration
17surf(Z,Y+yr,X) % Plotting the solid of revolution about y=yr
18hold on;
19plot([a b],[yr yr],'-r','LineWidth',2); % Plotting the line y=yr
20view(22,11); % 3-D graph viewpoint specification
21xlabel('X-axis');ylabel('Y-axis');zlabel('Z-axis');