Breaking Ice with Air cushion Vehicle - Find minimum pressure with Newton-Raphson method

Objective: CODE : import matplotlib.pyplot as plt import math import numpy as np beeta=0.5 r=40 sigma=150 def func(p,h): term1=pow(p,3)*(1-pow(beeta,2)) term2=(0.4*h*pow(beeta,2)-sigma*pow(h,2)/pow(r,2))*pow(p,2) term3=pow(sigma,2)*pow(h,4)/3/pow(r,4)*p term4=pow(sigma*pow(h,2)/3/pow(r,2),3) eq=term1+term2+term3+term4 return eq def…