newton raphson method

Objective: import matplotlib.pyplot as pyplot beta=0.5; r=40; sigma=150; h=0.6 def f(p): return (pow(p,3)*(1-pow(beta,2)))+(0.4*h*pow(beta,2))+(((sigma*pow(h,2))/(pow(r,2)))*p*p)+(p*((pow(sigma,2)*pow(h,4))/(3*pow(r,4))))-(pow(((sigma*pow(h,2))/(3*pow(r,2))),3)) def fprime(p): return ((pow(p,2)*(1-pow(beta,2)))*3) + (((0.4*h*pow(beta,2))-((sigma*pow(h,2))/(pow(r,2))))*2*p)…