Jbears

lebebfkadal

Oct 12th, 2016
118
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Python 0.47 KB | None | 0 0
  1. def leb(num_points,x):
  2.    
  3.     data = get_x_i(num_points-1)
  4.     N=num_points-1
  5.     lg=0.0
  6.     lagrange_basis = numpy.ones(N + 1,dtype=float)
  7.     for i in xrange(N + 1):
  8.          for j in xrange(N + 1):
  9.             if i != j:
  10.                 print lagrange_basis
  11.                
  12.                 lagrange_basis[i] *= (x - data[j]) / (data[i] - data[j])
  13.        
  14.     for i in xrange(N+1):
  15.         lg+=numpy.abs(lagrange_basis[i])
  16.        
  17.    
  18.    
  19.     return lg
Advertisement
Add Comment
Please, Sign In to add comment