col1:=proc(x1,y1,x2,y2,x3,y3)
local solution;
solution:= solve({y1=m*x1+b,y2=m*x2+b},{m,b});
 print(solution);
if (y3=subs(solution,m*x3+b)) then print (colinear);
else print(not colinear);
fi;
end;