import numpy as np
import matplotlib.pyplot as plt
from scipy.optimize import nnls
A = np.array([[60, 90, 120], [30, 120, 90]])
b = np.array([67.5, 60])
x, rnorm = nnls(A, b)
print(x)
print(rnorm)
решение по определенной системе линейных уравнений в питоне