kill(all);

Van = facebook.com/LearningCalculus;

A:matrix([a,0,0],[-a,a,0],[a,0,-3*a]);

b:matrix([7],[5],[1]);

x:first(linsolve_by_lu(A,b));

I:matrix([1],[1],[1]);

x.I;