python算法3

牛顿迭代法求取方程根

def solution(a,b,c,d):
    x=1.5
    x0=x;
    f=a*x0*x0*x0+b*x0*x0+c*x0+d
    fd=3*a*x0*x0+2*b*x0+c
    h=f/fd
    x=x0-h
    while abs(x-x0)>=1e-5:
        x0=x
        f=a*x0*x0*x0+b*x0*x0+c*x0+d
        fd=3*a*x0*x0+2*b*x0+c
        h=f/fd
        x=x0-h
    return x
if __name__=="__main__":
    print("请输入方程的系数：")
    a,b,c,d=map(float,input().split())
    print("方程的系数为:",a,b,c,d)
    x=solution(a,b,c,d)
    print("方程的根为:x=%.6f"%x)