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

Objective: import matplotlib.pyplot as pl import math import numpy as np # Function def f(p,h): return (pow(p,3)*(1-pow(beta,2))) + (((0.4*h*pow(beta,2))-(sigma*pow(h,2)/pow(r,2)))*pow(p,2)) + (((pow(sigma,2)*pow(h,4))/(3*pow(r,4)))*p)-(pow((sigma*pow(h,2)/(3*pow(r,2))),3)) #Derivation of Function def fprime(p,h): return (3*(pow(p,2)*(1-pow(beta,2))))…