Breaking Ice with Air Cusion Vehicle

Objective: Case 1 : For relaxation factor = 1 import matplotlib.pyplot as plt def f(p,h): x1 = pow(p,3)*(1-pow(b,2)) x2 = (0.4*h*pow(b,2)) - ((s*pow(h,2))/pow(r,2)) x3 = (pow(s,2))*(pow(h,4))/(3*pow(r,4)) x4 = s*pow(h,2)/(3*pow(r,2)) x5 = pow(x4,3) function = x1 + pow(p,2)*x2 + p*x3 - x5 return function def fprime(p,h): x11…