Write your own script file in MATLAB for the following problem:

A cylindrical silo with radius r has a spherical cap roof with radius R. The height of the cylindrical portion is H. The height of the cap is h. (See sketch below)

Write a program that determines the cylindrical height H for given values of r, R and the total volume of the silo V.

You will find the following equations useful:

Vcap=Volume of cap=(1/3)pi(h^2)(3R-h)

h=height of cap=R(1-cos(theta))

sin(theta)=r/R theta is an angle formed by h and R in the cap.

H=(V-Vcap)/(pi*r^2)

Use the input command to enter the values for r, R, and V.

TEST your program with the following values: r=30 feet, R=45 feet, and a volume V=120,000 cubic feet

DELIVERABLES:

All spreadsheets and graphs showing data. Matlab script file and output run