Newton-Raphson method
Objective: PROGRAM: 1 \"\"\" Newton - Raphson method \"\"\" import matplotlib.pyplot as plt 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)*p/(3*pow(r,4))-pow((sigma*pow(h,2))/(3*pow(r,2)),3) def fprime(p,h): return 3*pow(p,2)*(1-pow(beta,2))+(0.4*h*pow(beta,2)-sigma*pow(h,2)/pow(r,2))*2*p+pow(sigma,2)*pow(h,4)/(3*pow(r,4))…