Как запрограммировать (C++, Delphi, VB) нахождение переменных в системе уравнений?

k11*x1 + k12*x2 + .+kn1*xn = s1
k12*x1 + k22*x2 + .+kn2*xn = s2
k13*x1 + k23*x2 + .+kn3*xn = s3
k14*x1 + k24*x2 + .+kn4*xn = s4
k15*x1 + k25*x2 + .+kn5*xn = S

Числа n, s1, s2, s3, s4 - задаются пользователем.
Коэффициенты kij (0

Числа n, s1, s2, s3, s4 - задаются пользователем.
Коэффициенты kij (0

Числа n, s1, s2, s3, s4 - задаются пользователем.
Коэффициенты kij занесены в систему и являются константами.
Необходимо найти такие x1, x2, ..xn (необязательно целочисленные) , чтобы сумма S была минимальной.