Newton Raphson method for ice breaker

Objective: import matplotlib.pyplot as plt import math import xlwt from xlwt import Workbook def f(p,beta,r,sigma,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)*p/(3*pow(r,4))-pow((sigma*pow(h,2)/(3*pow(r,2))),3) def fprime(p,beta,r,sigma,h): return 3*pow(p,2)*(1-pow(beta,2))…