MATLAB please

Write a 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. 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 = h (3R - h) b-height of cap=R(1-cos(O)) is an angle formed by h and R in the cap. sin(0) /R HV-Veap nr 2 NOTES: * Make sure you can solve this problem using your calculator first. If you don't know how to do it with your calculator, you won't be able to write the script file. * 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=200,000 cubic feet Sketch of silo