A=[1 2 3;-1 3 7;9 0 3];

b=[1,4,7];

#方法1
x1=A\b;

#方法2
x2=linsolve(A,b)

using TyOptimization

f=[-40,-30];
A=[3 2];
b=[1000];
lb=[0,0];
ub=[Inf,Inf];

x,fval=linprog(f,A,b,[],[],lb,ub)
fval=-fval