勒让德多项式(Python实现)

下面是使用了sympy,其中x为设置的symbols。
point是一个整数的参数,我把这个放到了全局变量中,方便修改
也可以放到函数参数中。
代表的是有多少个零点的勒让德多项式。

def P():
    if point == 0:
        return 1
    elif point == 1:
        return x
    p0 = 1
    p1 = x
    for i in range(point - 1):
        temp = ((2 * i + 3) * x * p1 - (i + 1) * p0) / (i + 2)
        p0 = p1
        p1 = temp
    return p1

你可能感兴趣的:(Python,数值计算)