Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- def leb(num_points,x):
- data = get_x_i(num_points-1)
- N=num_points-1
- lg=0.0
- lagrange_basis = numpy.ones(N + 1,dtype=float)
- for i in xrange(N + 1):
- for j in xrange(N + 1):
- if i != j:
- print lagrange_basis
- lagrange_basis[i] *= (x - data[j]) / (data[i] - data[j])
- for i in xrange(N+1):
- lg+=numpy.abs(lagrange_basis[i])
- return lg
Advertisement
Add Comment
Please, Sign In to add comment