One option:
syms x y
eqn1 = x+y == 250;
eqn2= 10*x+45*y ==24000;
eqn3 = x+3*y ==3360;
eqn4 = 10*x+15*y ==1440;
eqns = [eqn1; eqn2; eqn3; eqn4];
vars = symvar(eqns)
[A,b] = equationsToMatrix(eqns,vars);
B = linsolve(double(A),double(b))
x = B(1)
y = B(2)
.
