Matlab Project – 4

Name :- Gineeth.R

Registration Number :- 24BAI1017

Matlab code for Volume of solid of revolution :-
clc

clearvars

syms x ;

f = input('Enter the function:');

fL = input('Enter the interval on which the function is defined:');

yr = input('Enter the axis of rotation y=c (Enter only c value):');

iL = input('Enter the integration limits:');

Volume = pi*int((f-yr)^2,iL(1),iL(2));

disp(['Volume is:', num2str(double(Volume))])

fx = inline(vectorize(f));

xvals = linspace(fL(1),fL(2),201);

xvalsr = fliplr(xvals);

xivals = linspace(iL(1),iL(2),201);

xivalsr = fliplr(xivals);

xlim = [fL(1) fL(2)+0.5];

ylim = fx(xlim);

figure('Position',[100 200 560 420])

subplot(2,1,1)
hold on;

plot(xvals,fx(xvals),'-b','Linewidth',2);

fill([xvals xvalsr],[fx(xvals) ones(size(xvalsr))*yr],[0.8 0.8 0.8], 'FaceAlpha',0.8)

plot([fL(1) fL(2)], [yr yr], '-r','LineWidth',2);

legend('Fuction Plot','Filled Region','Axis of Rotation','Location','Best');

title('Function y=f(x) and Region');

set(gca,'XLim',xlim)

xlabel('x-axis');

ylabel('y-axis');

subplot(2,1,2)

hold on;

plot(xivals,fx(xivals),'-b','Linewidth',2);

fill([xivals xivalsr],[fx(xivals) ones(size(xivalsr))*yr],[0.8 0.8 0.8],'FaceAlpha',0.8)

fill([xivals xivalsr],[ones(size(xivals))*yr -fx(xivalsr)+2*yr],[1 0.8 0.8],'FaceAlpha',0.8)

plot(xivals,-fx(xivals)+2*yr,'-m','Linewidth',2);

plot([iL(1) iL(2)],[yr yr],'-r','Linewidth',2);

title('Rotated Region in xy-Plane');

set(gca,'XLim',xlim)

xlabel('x-axis');

ylabel('y-axis');

[X,Y,Z] = cylinder(fx(xivals)-yr,100);

figure('Position',[700 200 560 420])

Z = iL(1) + Z.*(iL(2)-iL(1));

surf(Z,Y+yr,X,'EdgeColor','none','FaceColor','flat','FaceAlpha',0.6);
hold on;

plot([iL(1) iL(2)],[yr yr],'-r','Linewidth',2);

xlabel('X−axis');

ylabel('Y−axis');

zlabel('Z−axis');

view(22,11);

Q1) Visualize and find the volume of the region in the first quadrant bounded
above by the line y=sqrt(x), below by the curve y=sec(x)tan(x), and on the left by
the y-axis, about the line y=sqrt(2).

Ans)
Graph :-
In Matlab :-

In Command Window :-
Enter the function:

sec(x)*tan(x)

Enter the interval on which the function is defined:

[0 pi/4]

Enter the axis of rotation y=c (Enter only c value):

sqrt(2)

Enter the integration limits:

[0 pi/4]

Volume is:2.3014

In Figure Window :-
Q2) Visualize and find the volume of the solid generated by revolving region
bounded by curve y=sin(x), 0<=x<= π, about the line y=c, 0<=c<=1 by taking
c=0,0.2,0.4,0.6,0.8,1. Can you identify the range/exact value of c that maximize
and minimize the volume of the solid?

Ans)
Graph :-

C=0 :-

C=0.2 :-
C=0.4 :-

C=0.6 :-
C=0.8 :-

C=1 :-
C=0.7 :-
C=0.65 :-

In Matlab :-

In Command Window :-
i) c=0 :-

Enter the function:

sin(x)

Enter the interval on which the function is defined:

[0 pi]

Enter the axis of rotation y=c (Enter only c value):

Enter the integration limits:

[0 pi]

Volume is:4.9348
In Figure Window :-

In Command Window :-
ii) c=0.2 :-

Enter the function:

sin(x)
Enter the interval on which the function is defined:

[0 pi]

Enter the axis of rotation y=c (Enter only c value):

0.2

Enter the integration limits:

[0 pi]

Volume is:2.8163

In Figure Window :-

In Command Window :-
iii) c=0.4 :-

Enter the function:

sin(x)

Enter the interval on which the function is defined:

[0 pi]

Enter the axis of rotation y=c (Enter only c value):

0.4

Enter the integration limits:

[0 pi]

Volume is:1.4874

In Figure Window :-

In Command Window :-
iv) c=0.6 :-

Enter the function:

sin(x)

Enter the interval on which the function

is defined:

[0 pi]

Enter the axis of rotation y=c (Enter only

c value):

0.6

Enter the integration limits:

[0 pi]
Volume is:0.94804

In Figure Window :-
In Command Window :-
v) c=0.8 :-

Enter the function:

sin(x)

Enter the interval on which the function is defined:

[0 pi]

Enter the axis of rotation y=c (Enter only c value):

0.8

Enter the integration limits:

[0 pi]

Volume is:1.1983

In Figure Window :-

In Command Window :-
vi) c=1 :-

Enter the function:

sin(x)

Enter the interval on which the

function is defined:

[0 pi]
Enter the axis of rotation y=c (Enter only c value):

Enter the integration limits:

[0 pi]

Volume is:2.238

In Figure Window :-
In Command Window :-
vii) c=0.7 :-

Enter the function:

sin(x)

Enter the interval on which the function is defined:

[0 pi]

Enter the axis of rotation y=c (Enter only c value):

0.7

Enter the integration limits:

[0 pi]

Volume is:0.97445

In Figure Window :-

In Command Window :-
viii) c=0.65 :-

Enter the function:

sin(x)

Enter the interval on which the

function is defined:

[0 pi]

Enter the axis of rotation y=c (Enter only c value):

0.65

Enter the integration limits:

[0 pi]

Volume is:0.93657

In Figure Window :-
Q3) Modify the above code appropriately to visualize and find the volume of the
solid of revolution by the curve y=tan(πy/4), 0<=y<=1 about the y-axis.
Ans)
Graphs :-
In Matlab :-

In Command Window :-
Enter the function:

4/pi*(atan(x))

Enter the interval on which the function is defined:

[0 1]

Enter the axis of rotation y=c (Enter only c value):

Enter the integration limits:

[0 1]

Volume is:1.2492

In Figure Window :-

---THE
END---