The function, |
|
|
|
1 import numpy as np 2 D=np.array([[ -0.736], [2.480], [3.000], [2.000]]) 3 print(np.transpose(D)) 4 # [[-0.736 2.48 3. 2. ]] 5 import sympy as sym 6 x = sym.Symbol('x') 7 P3=D[0]+(x-0.3)*D[1]+(x-0.3)*(x-1)*D[2]+(x-0.3)*(x-1)*(x-0.7)*D[3] 8 print(P3) 9 # [2.48*x + 2.0*(x - 1)*(x - 0.7)*(x - 0.3) + 3.0*(x - 1)*(x - 0.3) - 1.48] 10 print(sym.expand(2.48*x + 2.0*(x - 1)*(x - 0.7)*(x - 0.3) + 3.0*(x - 1)*(x - 0.3) - 1.48)) 11 # 2.0*x**3 - 1.0*x**2 + 1.0*x - 1.0 |