def NR(p,q,r,s,t,u,v,w,x,y): n=0 p=p+0.1-0.1 q=q+0.1-0.1 r=r+0.1-0.1 s=s+0.1-0.1 t=t+0.1-0.1 u=u+0.1-0.1 a=input("first approximation=") while n<200: y=p+(q*a)+(r*((a)**2))+(s*((a)**3))+(t*((a)**4))+(u((a)**5))+(v*((a)**6))+(w*((a)**7))+(x*((a)**8))+(y*((a)**9)) dydx=q+(2*r*a)+(3*s*(a**2))+(4*t*(a**3))+(5*u*(a**4))+(6*v*(a**5))+(7*w*(a**6))+(8*x*(a**7))+(9*y*(a**8)) root=a-(y/dydx) a=root n=n+1 print ("root="),root