clear;
clc;
a = [1 2];
b = [-1 1];
c = sym('c',[2 1]);
eq1 = a*c==-1;
eq2 = b*c==3;
sol = vpasolve([eq1, eq2],c);
solc = [sol.c1 sol.c2];
disp(solc)