A=[1 2 3]; B=[4 5 6]; C=[7 8 9]; mat=[A;B;C]; if rank(mat)==3 fprintf('Uklad jest liniowo niezalezny'); return; end syms x y; f1=x*B+y*C-A; f2=x*A+y*C-B; f3=x*A+y*B-C; [x1,y1]=solve(f1); [x2,y2]=solve(f2); [x3,y3]=solve(f3); v1=[x1 y1] v2=[x2 y2] v3=[x3 y3]