编写程序，使用牛顿迭代法求方程在x附近的一个实根defsolution(a,b,c,d):x=1.5x0=xf=a*x0**3+b*x0**2+c*x0+dfd=3*a*x0**2+2*b*x0+ch=f/fdx=x0-hwhileabs(x-x0)>=1e-5:x0=x