Řešte soustavu lineárních rovnic:
                            2*x+y+3*z=1
                            x+4*y-2*z=-3

Řešení příslušné homogenní soustavy:

r1h:2*x+y+3*z=0; r2h:x+4*y-2*z=0;

linsolve_params:false$
Resenih:linsolve([r1h,r2h],[x,y,z]);

W:ev([x,y,z],Resenih);

load(draw)$

draw3d(surface_hide=true,xaxis=true,yaxis=true,zaxis=true,line_width=5,color=orange,parametric(W[1],W[2],W[3],z,-5,5));

draw3d(surface_hide=true,xaxis=true,yaxis=true,zaxis=true,line_width=5,parametric(W[1],W[2],W[3],z,-5,5),line_width=1,color=yellow,implicit(r1h,x,-10,10,y,-10,10,z,-10,10),color=pink,implicit(r2h,x,-10,10,y,-10,10,z,-10,10));

Řešení dané nehomogenní soustavy:

r1:2*x+y+3*z=1; r2:x+4*y-2*z=-3;

linsolve_params:false$
Reseni:linsolve([r1,r2],[x,y,z]);

K:ev([x,y,z],Reseni);

draw3d(surface_hide=true,xaxis=true,yaxis=true,zaxis=true,line_width=5,color=red,parametric(K[1],K[2],K[3],z,-5,5),color=orange,parametric(W[1],W[2],W[3],z,-5,5));

draw3d(surface_hide=true,xaxis=true,yaxis=true,zaxis=true,line_width=5,color=red,parametric(K[1],K[2],K[3],z,-5,5),line_width=1,color=blue,implicit(r1,x,-10,10,y,-10,10,z,-10,10),color=green,implicit(r2,x,-10,10,y,-10,10,z,-10,10),line_width=5,color=orange,parametric(W[1],W[2],W[3],z,-5,5),line_width=1,color=yellow,implicit(r1h,x,-10,10,y,-10,10,z,-10,10),color=pink,implicit(r2h,x,-10,10,y,-10,10,z,-10,10));