> #bellenot\
with(plots):
--------------------------------------------------------------------------------
> #1\
poly1:=[[1,-2],[3,-2],[3,2],[-3,2]];\
line1:=[[1,-2],[1,2]];\
PLOT(CURVES(line1, LINESTYLE(2)) ,POLYGONS(poly1,COLOUR(RGB,0.7,0.7,0.7)),\
TITLE(`By the Good Doctor`),AXESLABELS(`x`,`y`));

                  poly1 := [[1, -2], [3, -2], [3, 2], [-3, 2]]

                           line1 := [[1, -2], [1, 2]]
--------------------------------------------------------------------------------
> #2\
n:=20;\
poly2:=[seq([(-1+2*i/n)^2,-1+2*i/n],i=0..n)];\
PLOT(POLYGONS(poly2),TITLE(`By the Good Doctor`),AXESLABELS(`x`,`y`));

                                     n := 20

                     81            16            49
poly2 := [[1, -1], [---, -9/10], [----, -4/5], [---, -7/10], [9/25, -3/5],
                    100            25           100

    [1/4, -1/2], [4/25, -2/5], [9/100, -3/10], [1/25, -1/5], [1/100, -1/10],

    [0, 0], [1/100, 1/10], [1/25, 1/5], [9/100, 3/10], [4/25, 2/5], [1/4, 1/2],

                   49           16           81
    [9/25, 3/5], [---, 7/10], [----, 4/5], [---, 9/10], [1, 1]]
                  100           25          100
> #3\
poly3:=[[1,-2],[1,2],[-1,2],[-1,-2]];\
PLOT(POLYGONS(poly2,poly3,COLOUR(RGB,1,1,1,0.7,0.7,0.7)),TITLE(`By the Good Doctor`),AXESLABELS(`x`,`y`));

                  poly3 := [[1, -2], [1, 2], [-1, 2], [-1, -2]]
> #4\
poly3d:=[[2,0,0],[0,3,0],[0,0,6]];\
PLOT3D(POLYGONS(poly3d),STYLE(PATCH),AXES(NORMAL),\
TITLE(`By the Good Doctor`),AXESLABELS(`x`,`y`,`z`));\


                   poly3d := [[2, 0, 0], [0, 3, 0], [0, 0, 6]]
> #5\
n:=20;m:=20;\
f:=(x,y)->1 - ((y-1)^2)/3 - x/20;\
mesh:=[seq([seq([i/n,2*j/m,f(i/n,2*j/n)],i=0..n)],j=0..m)]:\
regionD:=[[0,0,0],[1,0,0],[1,2,0],[0,2,0]];\
PLOT3D(MESH(mesh),POLYGONS(regionD),STYLE(PATCH),AXES(NORMAL),\
TITLE(`By the Good Doctor`),AXESLABELS(`x`,`y`,`z`));

                                     n := 20

                                     m := 20

                                                  2
                     f := (x,y) -> 1 - 1/3 (y - 1)  - 1/20 x

             regionD := [[0, 0, 0], [1, 0, 0], [1, 2, 0], [0, 2, 0]]
--------------------------------------------------------------------------------
> 
--------------------------------------------------------------------------------