Matlab Assignment Solution on Complex Loops

• 5th Aug, 2022
• 15:20 PM

function A10Prob1_noloop_sphereVol_kang401(radius)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% ENGR 132 
% Program Description 
% The code computes volume of a sphere as a function of height of fluid
% The input to function call is radius
% The function does not return any variable
% variable h and vol are vectors of height and volume
% volume is computed as a function of each height
% isempty() function is used to check if vector vol is empty
% If vol is empty, then do not plot
% The code will plot if radius is >=100 and <=500 in cm
%
% Function Call
?lls function A10Prob1_noloop_sphereVol_kang401(radius) with input as radius
%
% Input Arguments
% Takes radius as input
%
% Output Arguments
% The function does not return any output.
% However, it displays plot of volume as a function of height
%
% Assignment Information
%   Assignment:     A10, Problem 1
%   Author:         Drew Kang
%   Team ID:        ###-##
%   Academic Integrity:
%     [] I worked with one or more peers but our collaboration
%        maintained academic integrity.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%% ____________________
%% VARIABLE INITIALIZATION
format longe
increment = 2*radius/100;
height = 0:increment:2*radius;
volume = []; % Initialize
%% ____________________
%% NESTED STRUCTURE

if radius>=100 && radius<=500
volume = pi*radius*height.^2 - (pi/3)*height.^3; % compute volume
disp('----------------------------------------------------')
disp('       Heignt (cm)                 Volume(cm^3)')
disp('----------------------------------------------------')
end

if isempty(volume) == 0 % isempty() is a Matlab function to check vol is empty
disp([height' volume'])
plot(height,volume)
xlabel('Height, cm')
ylabel('Volume,cm^3')
title('Filled volume of a sphere')
else
    display('Give radius between 100 and 500 cm')
end
end

%% ____________________
%% RESULTS
% >>A10Prob1_noloop_sphereVol_kang401(250)
% ----------------------------------------------------
%        Heignt (cm)                 Volume(cm^3)
% ----------------------------------------------------
%                          0                         0
%     5.000000000000000e+000    1.950405439103663e+004
%     1.000000000000000e+001    7.749261878854822e+004
%     1.500000000000000e+001    1.731802950291373e+005
%     2.000000000000000e+001    3.057816849494065e+005
%     2.500000000000000e+001    4.745113903859583e+005
%     3.000000000000000e+001    6.785840131753953e+005
%     3.500000000000000e+001    9.172141551543200e+005
%     4.000000000000000e+001    1.189616418159335e+006
%     4.500000000000000e+001    1.495005404027043e+006
%     5.000000000000000e+001    1.832595714594046e+006
%     5.500000000000000e+001    2.201601951696947e+006
%     6.000000000000000e+001    2.601238717172348e+006
%     6.500000000000000e+001    3.030720612856853e+006
%     7.000000000000000e+001    3.489262240587064e+006
%     7.500000000000000e+001    3.976078202199582e+006
%     8.000000000000000e+001    4.490383099531011e+006
%     8.500000000000000e+001    5.031391534417953e+006
%     9.000000000000000e+001    5.598318108697010e+006
%     9.500000000000000e+001    6.190377424204788e+006
%     1.000000000000000e+002    6.806784082777885e+006
%     1.050000000000000e+002    7.446752686252905e+006
%     1.100000000000000e+002    8.109497836466452e+006
%     1.150000000000000e+002    8.794234135255128e+006
%     1.200000000000000e+002    9.500176184455534e+006
%     1.250000000000000e+002    1.022653858590427e+007
%     1.300000000000000e+002    1.097253594143795e+007
%     1.350000000000000e+002    1.173738285289316e+007
%     1.400000000000000e+002    1.252029392210652e+007
%     1.450000000000000e+002    1.332048375091462e+007
%     1.500000000000000e+002    1.413716694115407e+007
%     1.550000000000000e+002    1.496955809466147e+007
%     1.600000000000000e+002    1.581687181327341e+007
%     1.650000000000000e+002    1.667832269882651e+007
%     1.700000000000000e+002    1.755312535315737e+007
%     1.750000000000000e+002    1.844049437810259e+007
%     1.800000000000000e+002    1.933964437549876e+007
%     1.850000000000000e+002    2.024978994718251e+007
%     1.900000000000000e+002    2.117014569499042e+007
%     1.950000000000000e+002    2.209992622075910e+007
%     2.000000000000000e+002    2.303834612632515e+007
%     2.050000000000000e+002    2.398462001352517e+007
%     2.100000000000000e+002    2.493796248419578e+007
%     2.150000000000000e+002    2.589758814017355e+007
%     2.200000000000000e+002    2.686271158329512e+007
%     2.250000000000000e+002    2.783254741539708e+007
%     2.300000000000000e+002    2.880631023831601e+007
%     2.350000000000000e+002    2.978321465388854e+007
%     2.400000000000000e+002    3.076247526395125e+007
%     2.450000000000000e+002    3.174330667034076e+007
%     2.500000000000000e+002    3.272492347489367e+007
%     2.550000000000000e+002    3.370654027944659e+007
%     2.600000000000000e+002    3.468737168583611e+007
%     2.650000000000000e+002    3.566663229589881e+007
%     2.700000000000000e+002    3.664353671147134e+007
%     2.750000000000000e+002    3.761729953439029e+007
%     2.800000000000000e+002    3.858713536649223e+007
%     2.850000000000000e+002    3.955225880961379e+007
%     2.900000000000000e+002    4.051188446559158e+007
%     2.950000000000000e+002    4.146522693626218e+007
%     3.000000000000000e+002    4.241150082346220e+007
%     3.050000000000000e+002    4.334992072902826e+007
%     3.100000000000000e+002    4.427970125479694e+007
%     3.150000000000000e+002    4.520005700260484e+007
%     3.200000000000000e+002    4.611020257428858e+007
%     3.250000000000000e+002    4.700935257168476e+007
%     3.300000000000000e+002    4.789672159662998e+007
%     3.350000000000000e+002    4.877152425096084e+007
%     3.400000000000000e+002    4.963297513651395e+007
%     3.450000000000000e+002    5.048028885512589e+007
%     3.500000000000000e+002    5.131268000863329e+007
%     3.550000000000000e+002    5.212936319887274e+007
%     3.600000000000000e+002    5.292955302768083e+007
%     3.650000000000000e+002    5.371246409689418e+007
%     3.700000000000000e+002    5.447731100834940e+007
%     3.750000000000000e+002    5.522330836388308e+007
%     3.800000000000000e+002    5.594967076533183e+007
%     3.850000000000000e+002    5.665561281453222e+007
%     3.900000000000000e+002    5.734034911332090e+007
%     3.950000000000000e+002    5.800309426353445e+007
%     4.000000000000000e+002    5.864306286700947e+007
%     4.050000000000000e+002    5.925946952558257e+007
%     4.100000000000000e+002    5.985152884109034e+007
%     4.150000000000000e+002    6.041845541536939e+007
%     4.200000000000000e+002    6.095946385025634e+007
%     4.250000000000000e+002    6.147376874758779e+007
%     4.300000000000000e+002    6.196058470920028e+007
%     4.350000000000000e+002    6.241912633693050e+007
%     4.400000000000000e+002    6.284860823261501e+007
%     4.450000000000000e+002    6.324824499809040e+007
%     4.500000000000000e+002    6.361725123519333e+007
%     4.550000000000000e+002    6.395484154576030e+007
%     4.600000000000000e+002    6.426023053162801e+007
%     4.650000000000000e+002    6.453263279463303e+007
%     4.700000000000000e+002    6.477126293661197e+007
%     4.750000000000000e+002    6.497533555940139e+007
%     4.800000000000000e+002    6.514406526483795e+007
%     4.850000000000000e+002    6.527666665475822e+007
%     4.900000000000000e+002    6.537235433099879e+007
%     4.950000000000000e+002    6.543034289539632e+007
%     5.000000000000000e+002    6.544984694978735e+007

%% ____________________
%?ADEMIC INTEGRITY STATEMENT
% I have not used source code obtained from any other unauthorized
% source, either modified or unmodified. I have not provided
?cess to my code to anyone in any way. The program I am 
% submitting is my own original work.