use plots in
module()
local q,p1,p2,p3;
q:=9*x^2+z^2-18*x-6*z+9;
p1:=plots:-implicitplot3d(9*x^2+z^2-18*x-6*z+9=0,x=0..2,y=1..5,z=0..6, style=surface,grid=[25,25,25]);
p2:=spacecurve([[1,0,3],[1,6,3]],color=red,thickness=3,numpoints=2);
p3:=display(p1,p2,scaling=constrained,axes=frame,view=[0..2,0..6,0..6], tickmarks=[3,6,6],orientation=[-35,70,0]);
print(p3);
end module:
end use: