1.引入库
import numpy as np
2.正文
import numpy as np enter = input('请输入若干个三维数组,分别表示x,y,以及一阶导数的值用逗号隔开:') data = enter.split(',') data = np.array(data,dtype=float) ji = 0 for i in data: ji+=1 n = int(ji/3) data = data.reshape(n,3) data_x = np.array([]) data_y = np.array([]) data_d = np.array([]) for i in range(0,n): xi = data[i,0] yi = data[i,1] zi = data[i,2] data_x = np.append(data_x,[xi],axis=0) data_y = np.append(data_y,[yi],axis=0) data_d = np.append(data_d,[zi],axis=0) def lagrange_basis_function(x,n): a1 = np.array([],dtype=float) for i in range(0,n): p,q=1,1 for j in range(0,i): q *= x-data_x[j] p *= data_x[i]-data_x[j] for j in range(i+1