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

Objective: \"\"\" Breaking ice with air cushion \"\"\" import matplotlib.pyplot as plt import math import numpy as np #function f(x,y) 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)) #derivative…