def basis_polys():
""" Returns the basis polynomials and their derivatives and antiderivatives
"""
nodes, _, _ = quad()
? = [lagrange(nodes,eye(N+1)[i]) for i in range(N+1)]
?Der = [[polyder(?_p, m=a) for ?_p in ?] for a in range(N+1)]
?Int = [polyint(?_p) for ?_p in ?]
return ?, ?Der, ?Int
评论列表
文章目录