def hermite_basis_poly(x_values, index): lag_basis = langrange_basis_poly(x_values, index) value = lag_basis.deriv()(x_values[index]) poly = 1 - 2 * value * poly1d(x_values[index], True) herm_basis = (lag_basis ** 2) * poly return herm_basis
def hermite_basis_deriv_poly(x_values, index): lag_basis = langrange_basis_poly(x_values, index) poly = poly1d(x_values[index], True) herm_basis_deriv = (lag_basis ** 2) * poly return herm_basis_deriv