Optimization of Ice Pressure using Python

Objective: import numpy import scipy import matplotlib.pyplot as pyplot def f(p,h): t1 = pow(p,3)*(1 - pow(b,2)) t2 = pow(p,2)*((0.4*h*pow(b,2) - sigma*pow(h,2)/pow(r,2))) t3 = (pow(sigma,2)*pow(h,4)/(3*pow(r,4)))*p t4 = sigma*pow(h,2)/(3*pow(r,2)) return t1+t2+t3+t4 def df(p,h): dt1 = 3*pow(p,2)*(1 - pow(b,2)) dt2 = 2*p*((0.4*h*pow(b,2)…