>>>> 4 * answers[0] + 5 * answers[1]
>20.0
>>>> 1 * answers[0] + 2 * answers[1]
>13.0
Так и есть. Для того чтобы напечатать меньше текста, вы также можете указать NumPy найти скалярное произведение массивов:
>>>> product = np.dot(coefficients, answers)
>>>> product
>array([20.,··13.])
Если решение верно, значения массива product должны быть близки к значениям массива dependents. Вы можете использовать функцию allclose(), чтобы проверить, являются ли массивы хотя бы приблизительно равными (они могут быть не полностью равными из-за округления чисел с плавающей точкой):
>>>> np.allclose(product, dependents)
>True
NumPy также имеет модули для работы с многочленами, преобразованиями Фурье, статистикой и распределением вероятностей.