restart;
f:=x^2: g:=x^4: h:=x^3:
plot([f,g,h],x=0..2);
g:=9*x^2+25*y^2-18*x-100*y-116=0;
with(plots);
implicitplot(g, x=-10..10, y=-10..10, scaling=constrained);
x := 2 *cos(t); y := 2 *sin(t);
plot([x, y, t = 0..2*Pi], scaling = constrained);
restart;
solve([x^2-y-5=0, x-y+1=0], [x, y]);
h:=x^2+y^2;
plot3d(h, x=-1..5, y=-2..8);
with(plots);
implicitplot3d(h=z, x=-5..5, y=-5..5, z=0..50, grid = [5, 25, 25]);
x:=2*cos(a)*cos(b); y:=3*sin(a)*cos(b); z:=4*sin(b);
plot3d([x,y,z], a=0..2*Pi, b=-Pi/2..Pi/2, scaling=constrained);