restart: randomize(): with(Maplets[Tools]): with(Maplets[Elements]): with(plots):with(student):with(StringTools): StartEngine(); QSurfaces:= Maplet(onstartup = RunWindow(MAIN), Window[MAIN](title="Quadric Surfaces", [halign=none, [ Button("New Question", Evaluate(function = "askNQ")), HorizontalGlue(), Button("Quit", Shutdown()) ], [ border=true, caption="Consider the quadric equation:", MathMLViewer[EQN](width=500, height=60, background=white) ], [ border=true, caption="Write the equation in standard form:", TextField[stdform](width=60), Button("Check", Evaluate(function = "check_std")), Button("Show", Evaluate(function = "show_std")) ], [ border=true, caption="Identify the type of quadric surface:", RadioButton[RB1]("Sphere",group=BG1), RadioButton[RB2]("Ellipsoid",group=BG1), RadioButton[RB3]("Hyperboloid",group=BG1), RadioButton[RB4]("Cone",group=BG1), RadioButton[RB5]("Paraboloid",group=BG1), RadioButton[RB6]("Cylinder",group=BG1), HorizontalGlue(), Button("Check", Evaluate(function = "check_shape")), Button("Show", Evaluate(function = "show_shape")) ], BoxColumn[typecol]( border=true, halign=none, valign=none, inset=0, spacing=0, caption="Properties of the surface:", [ # Sphere Label[Lspcenter]("Center: [p,q,r] =", visible=false), TextField[CenterS]("[ ]",width=10, visible=false), HorizontalGlue(), Label[Lspradius]("Radius: R =", visible=false), TextField[RadiusSp](width=5, visible=false), HorizontalGlue(), Button[Bspcheck]("Check", Evaluate(function = "check_sphere"), visible=false), Button[Bspshow]("Show", Evaluate(function = "show_sphere"), visible=false) ], [ # Ellipsoid Label[Lelcenter]("Center: [p,q,r] =",visible=false), TextField[CenterE]("[ ]",width=10, visible=false), HorizontalGlue(), Label[Lelrada]("Radii: a =",visible=false), TextField[RadiusEa](width=5, visible=false), Label[Lelradb]("b =", visible=false), TextField[RadiusEb](width=5, visible=false), Label[Lelradc]("c =", visible=false), TextField[RadiusEc](width=5, visible=false), HorizontalGlue(), Button[Bellcheck]("Check", Evaluate(function = "check_ellipsoid"), visible=false), Button[Bellshow]("Show", Evaluate(function = "show_ellipsoid"), visible=false) ], [ # Hyperboloid Label[Lhycenter]("Center: [p,q,r] =",visible=false), TextField[CenterH]("[ ]",width=10, visible=false), HorizontalGlue(), Label[Lhyaxis]("Axis:", visible=false), RadioButton[HyAxisx]("x",group=Axis, visible=false), RadioButton[HyAxisy]("y",group=Axis, visible=false), RadioButton[HyAxisz]("z",group=Axis, visible=false), HorizontalGlue(), Label[Lhysheets]("Sheets:", visible=false), RadioButton[onesheet]("1",group=Sheet, visible=false), RadioButton[twosheet]("2",group=Sheet, visible=false), HorizontalGlue(), Button[Bhycheck]("Check", Evaluate(function = "check_hyperboloid"), visible=false), Button[Bhyshow]("Show", Evaluate(function = "show_hyperboloid"), visible=false) ], [ # Cone Label[Lcovert]("Vertex: [p,q,r] =",visible=false), TextField[VertexC]("[ ]",width=10, visible=false), HorizontalGlue(), Label[Lcoaxis]("Axis:",visible=false), RadioButton[CoAxisx]("x-axis",group=Axis, visible=false), RadioButton[CoAxisy]("y-axis",group=Axis, visible=false), RadioButton[CoAxisz]("z-axis",group=Axis, visible=false), HorizontalGlue(), Button[Bcocheck]("Check", Evaluate(function = "check_cone"), visible=false), Button[Bcoshow]("Show", Evaluate(function = "show_cone"), visible=false) ], BoxColumn( valign=none, inset=0, spacing=0, [ # Paraboloid Label[Lparaaxis]("Axis:", visible=false), RadioButton[Paraaxisx]("x",group=Axis, visible=false), RadioButton[Paraaxisy]("y",group=Axis, visible=false), RadioButton[Paraaxisz]("z",group=Axis, visible=false), HorizontalGlue(), Label[Lparatype]("Type:", visible=false), RadioButton[RBP1]("Elliptic ",group=BGP, visible=false), RadioButton[RBP2]("Hyperbolic",group=BGP, visible=false), HorizontalGlue(), HorizontalGlue() ], [ Label[Lvert]("Vertex: [p,q,r] =", visible=false), TextField[ParaVertex]("[ ]",width=10, visible=false), HorizontalGlue(), Button[Bparacheck]("Check", Evaluate(function = "check_paraboloid"), visible=false), Button[Bparashow]("Show", Evaluate(function = "show_paraboloid"), visible=false) ] ), BoxColumn( valign=none, inset=0, spacing=0, [ # Cylinder Label[LAxiscyl]("Axis:", visible=false), RadioButton[Cylaxisx]("x",group=Axis, visible=false), RadioButton[Cylaxisy]("y",group=Axis, visible=false), RadioButton[Cylaxisz]("z",group=Axis, visible=false), HorizontalGlue(), Label[Lcyltype]("Type:", visible=false), RadioButton[RBC1]("Circular",group=BG2, visible=false), RadioButton[RBC2]("Elliptic",group=BG2, visible=false), RadioButton[RBC3]("Hyperbolic",group=BG2, visible=false), RadioButton[RBC4]("Parabolic",group=BG2, visible=false), HorizontalGlue(), HorizontalGlue() ], [ Label[Lcylcenter]("Center or Vertex: (Axis coordinate is 0.) [p,q,r] =", visible=false), TextField[Cylcenter]("[ ]",width=10, visible=false), HorizontalGlue(), Button[Bcylcheck]("Check", Evaluate(function = "check_cylinder"), visible=false), Button[Bcylshow]("Show", Evaluate(function = "show_cylinder"), visible=false) ] ) ), [ TextField[reply](width=10, editable=false, background=white) ], [ "Programmers: Shu Ho & Amy Su.", HorizontalGlue(), "Copyright 2008 P. Yasskin." ] ] ), ButtonGroup[BG1](), ButtonGroup[BG2](), ButtonGroup[BGP](), ButtonGroup[Axis](), ButtonGroup[Sheet]() ): rand6:=rand(1..6): rand2:=rand(1..2): randm33:=rand(-3..3): rand3:=rand(1..3): rand3:=rand(1..3): askNQ:=proc() global probtype, subprobtype, stdfrm,axis,radlen, a,b,c,g,h,i; local rgt, givenfrm; probtype:=rand6();#rand6(); subprobtype:=rand2(); g:=randm33(); h:=randm33(); i:=randm33(); a:=rand3(); b:=rand3(); c:=rand3(); radlen:=rand6()^2;#radius length rgt:=lcm(a^2,b^2,c^2); axis:=rand3(); Set(stdform = ""); Set(RB1=false); Set(RB2=false); Set(RB3=false); Set(RB4=false); Set(RB5=false); Set(RB6=false); if probtype=1 then # sphere stdfrm:=((x-g)^2)/radlen+((y-h)^2)/radlen+((z-i)^2)/radlen=1; givenfrm:=expand(stdfrm); elif probtype=2 then # ellipse stdfrm:=(x-g)^2/a^2+(y-h)^2/b^2+(z-i)^2/c^2=1; givenfrm:=expand(stdfrm)*rgt; elif probtype=3 then # hyperboloid if subprobtype=1 then #++- (zaxis) 1 sheet if axis=1 then stdfrm:=-(x-g)^2+(y-h)^2+(z-i)^2=1; givenfrm:=expand(stdfrm); elif axis=2 then stdfrm:=(x-g)^2-(y-h)^2+(z-i)^2=1; givenfrm:=expand(stdfrm); elif axis=3 then stdfrm:=(x-g)^2+(y-h)^2-(z-i)^2=1; givenfrm:=expand(stdfrm); end if; elif subprobtype=2 then #--+ (zaxis) 2 sheet if axis=1 then stdfrm:=(x-g)^2-(y-h)^2-(z-i)^2=1; givenfrm:=expand(stdfrm); elif axis=2 then stdfrm:=-(x-g)^2+(y-h)^2-(z-i)^2=1; givenfrm:=expand(stdfrm); elif axis=3 then stdfrm:= -(x-g)^2-(y-h)^2+(z-i)^2=1; givenfrm:=expand(stdfrm); end if; end if; elif probtype=4 then # Cone (x-a)^2 + (y-b)^2 = (z-c)^2 if axis = 1 then stdfrm:=(z-i)^2 + (y-h)^2 = (x-g)^2; givenfrm:=expand(stdfrm); elif axis = 2 then stdfrm:=(x-g)^2 + (z-i)^2 = (y-h)^2; givenfrm:=expand(stdfrm); elif axis = 3 then stdfrm:=(x-g)^2 + (y-h)^2 = (z-i)^2; givenfrm:=expand(stdfrm); end if; elif probtype=5 then # Paraboloid if subprobtype=1 then #++/-- (Elliptic)z-i=c+(x-g)^2/a^2+(y-h)^2/b^2 if axis=1 then stdfrm:=x-g=(z-i)^2/c^2+(y-h)^2/b^2; givenfrm:=expand(stdfrm); elif axis=2 then stdfrm:=y-h=(x-g)^2/a^2+(z-i)^2/c^2; givenfrm:=expand(stdfrm); elif axis=3 then stdfrm:=z-i=(x-g)^2/a^2+(y-h)^2/b^2; givenfrm:=expand(stdfrm); end if; elif subprobtype=2 then #+-/-+ (Saddle/Hyperbolic) z-i=c-(x-g)^2/a^2+(y-h)^2/b^2 if axis=1 then stdfrm:=x-g=(z-i)^2/c^2+(y-h)^2/b^2; givenfrm:=expand(stdfrm); elif axis=2 then stdfrm:=y-h=(x-g)^2/a^2-(z-i)^2/c^2; givenfrm:=expand(stdfrm); elif axis=3 then stdfrm:=z-i=(x-g)^2/a^2+(y-h)^2/b^2; givenfrm:=expand(stdfrm); end if; end if; elif probtype=6 then # Cylinder if subprobtype=1 then #(Circle)(x-g)^2+(y-h)^2=1; if axis=1 then stdfrm:=(z-i)^2+(y-h)^2=1; givenfrm:=expand(stdfrm); g:=0; elif axis=2 then stdfrm:=(x-g)^2+(z-i)^2=1; givenfrm:=expand(stdfrm); h:=0; elif axis=3 then stdfrm:=(x-g)^2+(y-h)^2=1; givenfrm:=expand(stdfrm); i:=0; end if; elif subprobtype=2 then #(Ellipse)(x-g)^2/a^2+(y-h)^2/b^2=1; if axis=1 then stdfrm:=(z-i)^2/c^2+(y-h)^2/b^2=1; givenfrm:=expand(stdfrm); g:=0; elif axis=2 then stdfrm:=(x-g)^2/a^2+(z-i)^2/c^2=1; givenfrm:=expand(stdfrm); h:=0; elif axis=3 then stdfrm:=(x-g)^2/a^2+(y-h)^2/b^2=1; givenfrm:=expand(stdfrm); i:=0; end if; end if; end if; Set(EQN = MathML[Export](givenfrm)); Set(reply=""); #HERE COMES THE MASSIVE LOADS OF VARIABLE FALSES! #Spheres: Set(Lspcenter(visible)=false); Set(CenterS(visible)=false, CenterS=""); Set(Lspradius(visible)=false); Set(RadiusSp(visible)=false, RadiusSp=""); Set(Bspcheck(visible)=false); Set(Bspshow(visible)=false); #Ellipsoid: Set(CenterE(visible)=false, CenterE=""); Set(RadiusEa(visible)=false, RadiusEa=""); Set(RadiusEb(visible)=false, RadiusEb=""); Set(RadiusEc(visible)=false, RadiusEc=""); Set(Bellcheck(visible)=false); Set(Bellshow(visible)=false); Set(Lelcenter(visible)=false); Set(Lelrada(visible)=false); Set(Lelradb(visible)=false); Set(Lelradc(visible)=false); #Hyperboloid Set(Lhycenter(visible)=false); Set(CenterH(visible)=false, CenterH=""); Set(Lhyaxis(visible)=false); Set(HyAxisx(visible)=false, HyAxisx=false); Set(HyAxisy(visible)=false, HyAxisy=false); Set(HyAxisz(visible)=false, HyAxisz=false); Set(onesheet(visible)=false, onesheet=false); Set(twosheet(visible)=false, twosheet=false); Set(Bhycheck(visible)=false); Set(Bhyshow(visible)=false); Set(Lhysheets(visible)=false); #Cone Set(Lcovert (visible)=false); Set(VertexC (visible)=false, VertexC=""); Set(Lcoaxis (visible)=false); Set(CoAxisx (visible)=false, CoAxisx=false); Set(CoAxisy (visible)=false, CoAxisy=false); Set(CoAxisz (visible)=false, CoAxisz=false); Set(Bcocheck (visible)=false); Set(Bcoshow(visible)=false); #Paraboloid Set(RBP1 (visible)=false, RBP1=false); Set(RBP2 (visible)=false, RBP2=false); Set(Lvert (visible)=false); Set(ParaVertex (visible)=false, ParaVertex=""); Set(Lparaaxis(visible)=false); Set(Paraaxisx(visible)=false, Paraaxisx=false); Set(Paraaxisy(visible)=false, Paraaxisy=false); Set(Paraaxisz(visible)=false, Paraaxisz=false); Set(Bparacheck (visible)=false); Set(Bparashow(visible)=false); Set(Lparatype(visible)=false); #Cylinder Set(RBC1 (visible)=false, RBC1=false); Set(RBC2 (visible)=false, RBC2=false); Set(RBC3 (visible)=false, RBC3=false); Set(RBC4 (visible)=false, RBC4=false); Set(LAxiscyl (visible)=false); Set(Cylaxisx (visible)=false, Cylaxisx=false); Set(Cylaxisy (visible)=false, Cylaxisy=false); Set(Cylaxisz (visible)=false, Cylaxisz=false); Set(Lcylcenter (visible)=false); Set(Cylcenter (visible)=false, Cylcenter=""); Set(Bcylcheck (visible)=false); Set(Bcylshow (visible)=false); Set(Lcyltype (visible)=false); end proc: check_std:=proc() global stdfrm; local useranswer; useranswer:= Get(stdform::equation): if simplify(lhs(useranswer)-lhs(stdfrm))=0 and simplify(rhs(useranswer)-rhs(stdfrm))=0 then Set(reply ="Good Job.") else Set(reply = "Try again.") end if; end proc: show_std:=proc() global stdfrm; Set(stdform = stdfrm); Set(reply="Standard form shown."); end proc: check_shape:=proc() global probtype; local usernum; if Get(RB1::truefalse) then usernum:=1 end if: if Get(RB2::truefalse) then usernum:=2 end if: if Get(RB3::truefalse) then usernum:=3 end if: if Get(RB4::truefalse) then usernum:=4 end if: if Get(RB5::truefalse) then usernum:=5 end if: if Get(RB6::truefalse) then usernum:=6 end if: if usernum=probtype then Set(reply="Correct! Keep going."); #MASSIVE SEGMENT OF REPETITION OCCURS NOW!! if probtype=1 then Set(Lspcenter(visible)=true); Set(CenterS(visible)=true); Set(Lspradius(visible)=true); Set(RadiusSp(visible)=true); Set(Bspcheck(visible)=true); Set(Bspshow(visible)=true); elif probtype=2 then Set(CenterE(visible)=true); Set(RadiusEa(visible)=true); Set(RadiusEb(visible)=true); Set(RadiusEc(visible)=true); Set(Bellcheck(visible)=true); Set(Bellshow(visible)=true); Set(Lelcenter(visible)=true); Set(Lelrada(visible)=true); Set(Lelradb(visible)=true); Set(Lelradc(visible)=true); elif probtype=3 then Set(Lhycenter(visible)=true); Set(CenterH(visible)=true); Set(Lhyaxis(visible)=true); Set(HyAxisx(visible)=true); Set(HyAxisy(visible)=true); Set(HyAxisz(visible)=true); Set(onesheet(visible)=true); Set(twosheet(visible)=true); Set(Bhycheck(visible)=true); Set(Bhyshow(visible)=true); Set(Lhysheets(visible)=true); elif probtype=4 then Set(Lcovert(visible)=true); Set(VertexC(visible)=true); Set(Lcoaxis(visible)=true); Set(CoAxisx(visible)=true); Set(CoAxisy(visible)=true); Set(CoAxisz(visible)=true); Set(Bcocheck(visible)=true); Set(Bcoshow(visible)=true); elif probtype=5 then Set(RBP1 (visible)=true); Set(RBP2 (visible)=true); Set(Lvert (visible)=true); Set(ParaVertex (visible)=true); Set(Lparaaxis (visible)=true); Set(Paraaxisx (visible)=true); Set(Paraaxisy (visible)=true); Set(Paraaxisz (visible)=true); Set(Bparacheck (visible)=true); Set(Bparashow(visible)=true); Set(Lparatype(visible)=true); elif probtype=6 then Set(RBC1 (visible)=true); Set(RBC2 (visible)=true); Set(RBC3 (visible)=true); Set(RBC4 (visible)=true); Set(LAxiscyl (visible)=true); Set(Cylaxisx (visible)=true); Set(Cylaxisy (visible)=true); Set(Cylaxisz (visible)=true); Set(Lcylcenter (visible)=true); Set(Cylcenter (visible)=true); Set(Bcylcheck (visible)=true); Set(Bcylshow (visible)=true); Set(Lcyltype (visible)=true); end if; else Set(reply="Incorrect. Try again.") end if; end proc: show_shape:=proc() global probtype; local usernum; Set(reply="Try it yourself next time."); #MASSIVE SEGMENT OF REPETITION OCCURS NOW!! if probtype=1 then Set(RB1=true); Set(Lspcenter(visible)=true); Set(CenterS(visible)=true); Set(Lspradius(visible)=true); Set(RadiusSp(visible)=true); Set(Bspcheck(visible)=true); Set(Bspshow(visible)=true); elif probtype=2 then Set(RB2=true); Set(CenterE(visible)=true); Set(RadiusEa(visible)=true); Set(RadiusEb(visible)=true); Set(RadiusEc(visible)=true); Set(Bellcheck(visible)=true); Set(Bellshow(visible)=true); Set(Lelcenter(visible)=true); Set(Lelrada(visible)=true); Set(Lelradb(visible)=true); Set(Lelradc(visible)=true); elif probtype=3 then Set(RB3=true); Set(Lhycenter(visible)=true); Set(CenterH(visible)=true); Set(Lhyaxis(visible)=true); Set(HyAxisx(visible)=true); Set(HyAxisy(visible)=true); Set(HyAxisz(visible)=true); Set(onesheet(visible)=true); Set(twosheet(visible)=true); Set(Bhycheck(visible)=true); Set(Bhyshow(visible)=true); Set(Lhysheets(visible)=true); elif probtype=4 then Set(RB4=true); Set(Lcovert(visible)=true); Set(VertexC(visible)=true); Set(Lcoaxis(visible)=true); Set(CoAxisx(visible)=true); Set(CoAxisy(visible)=true); Set(CoAxisz(visible)=true); Set(Bcocheck(visible)=true); Set(Bcoshow(visible)=true); elif probtype=5 then Set(RB5=true); Set(RBP1 (visible)=true); Set(RBP2 (visible)=true); Set(Lvert (visible)=true); Set(ParaVertex (visible)=true); Set(Lparaaxis (visible)=true); Set(Paraaxisx (visible)=true); Set(Paraaxisy (visible)=true); Set(Paraaxisz (visible)=true); Set(Bparacheck (visible)=true); Set(Bparashow(visible)=true); Set(Lparatype(visible)=true); elif probtype=6 then Set(RB6=true); Set(RBC1 (visible)=true); Set(RBC2 (visible)=true); Set(RBC3 (visible)=true); Set(RBC4 (visible)=true); Set(LAxiscyl (visible)=true); Set(Cylaxisx (visible)=true); Set(Cylaxisy (visible)=true); Set(Cylaxisz (visible)=true); Set(Lcylcenter (visible)=true); Set(Cylcenter (visible)=true); Set(Bcylcheck (visible)=true); Set(Bcylshow (visible)=true); Set(Lcyltype (visible)=true); end if; end proc: check_sphere:=proc() global g, h, i, CenterS, RadiusSp, radlen; local centerSp, radiusSp; #ghi centerSp:= Get(CenterS::anything); radiusSp:= Get(RadiusSp::integer); if radiusSp - sqrt(radlen) = 0 and centerSp = [g,h,i] then Set('reply' = "Correct! Try another question.") else Set('reply' = "Try again.") end if; end proc: show_sphere:=proc() global g,h,i, radlen; Set(CenterS = [g,h,i]); Set(RadiusSp =sqrt(radlen)); Set(reply = "Answer has been unveiled."); end proc: check_ellipsoid:=proc() global g,h,i,a,b,c; local centerEll, ellrada, ellradb, ellradc; centerEll:=Get(CenterE::anything); ellrada:=Get(RadiusEa::integer); ellradb:=Get(RadiusEb::integer); ellradc:=Get(RadiusEc::integer); if centerEll = [g,h,i] and ellrada = a and ellradb = b and ellradc = c then Set(reply = "You got it right!"); else Set(reply = "Try again."); end if; end proc: show_ellipsoid:=proc() global g,h,i,a,b,c; Set(CenterE = [g,h,i]); Set(RadiusEa = a); Set(RadiusEb = b); Set(RadiusEc = c); Set(reply = "The answer has been shown."); end proc: check_hyperboloid:=proc() global g, h, i, subprobtype, axis; local centerHy, useraxis, usersheet; centerHy:= Get(CenterH::anything); if Get(HyAxisx ::truefalse) then useraxis:=1 end if: if Get(HyAxisy ::truefalse) then useraxis:=2 end if: if Get(HyAxisz ::truefalse) then useraxis:=3 end if: if Get(onesheet::truefalse) then usersheet:=1 end if: if Get(twosheet::truefalse) then usersheet:=2 end if: if centerHy = [g,h,i] and useraxis = axis and usersheet = subprobtype then Set(reply = "Good job"); else Set(reply = "Try again."); end if; end proc: show_hyperboloid:=proc() global g, h, i, subprobtype, axis; Set(CenterH = [g,h,i]); if axis = 1 then Set(HyAxisx = true); elif axis = 2 then Set(HyAxisy = true); elif axis = 3 then Set(HyAxisz = true); end if; if subprobtype = 1 then Set(onesheet = true) elif subprobtype = 2 then Set(twosheet = true) end if; end proc: check_cone:=proc() global g,h,i, axis; local uservertex, useraxis; uservertex:=Get(VertexC::anything); if Get(CoAxisx ::truefalse) then useraxis:=1 end if: if Get(CoAxisy ::truefalse) then useraxis:=2 end if: if Get(CoAxisz ::truefalse) then useraxis:=3 end if: if useraxis = axis and uservertex = [g,h,i] then Set(reply = "Good job!"); else Set(reply = "Try again."); end if; end proc: show_cone:=proc() global g,h,i, axis; Set(VertexC = [g,h,i]); if axis = 1 then Set(CoAxisx = true); elif axis = 2 then Set(CoAxisy = true); elif axis = 3 then Set(CoAxisz = true); end if; end proc: check_paraboloid:=proc() global g, h, i, subprobtype, axis; local userVertex, useraxis, usertype; userVertex:= Get(ParaVertex::anything); if Get(Paraaxisx ::truefalse) then useraxis:=1 end if: if Get(Paraaxisy ::truefalse) then useraxis:=2 end if: if Get(Paraaxisz ::truefalse) then useraxis:=3 end if: if Get(RBP1::truefalse) then usertype:=1 end if: if Get(RBP2::truefalse) then usertype:=2 end if: if userVertex = [g,h,i] and useraxis = axis and usertype = subprobtype then Set(reply = "Good job"); else Set(reply = "Try again."); end if; end proc: show_paraboloid:=proc() global g,h,i,axis; Set(ParaVertex = [g,h,i]); if axis = 1 then Set(Paraaxisx = true); elif axis = 2 then Set(Paraaxisy = true); elif axis = 3 then Set(Paraaxisz = true); end if; if subprobtype = 1 then Set(RBP1 = true) elif subprobtype = 2 then Set(RBP2 = true) end if; end proc: check_cylinder:=proc() global g,h,i,axis, subprobtype; local usercenter, useraxis, usertype; usercenter:=Get(Cylcenter::anything); if Get(RBC1::truefalse) then usertype:=1 end if: if Get(RBC2::truefalse) then usertype:=2 end if: if Get(RBC3::truefalse) then usertype:=3 end if: if Get(RBC4::truefalse) then usertype:=4 end if: if Get(Cylaxisx::truefalse) then useraxis:=1 end if: if Get(Cylaxisy::truefalse) then useraxis:=2 end if: if Get(Cylaxisz::truefalse) then useraxis:=3 end if: if useraxis = axis and usercenter = [g,h,i] and usertype = subprobtype then Set(reply = "Good job!"); else Set(reply = "Try again."); end if; end proc: show_cylinder:=proc() global g,h,i, axis; Set(Cylcenter = [g,h,i]); if subprobtype = 1 then Set(RBC1 = true) elif subprobtype = 2 then Set(RBC2 = true) end if; if axis = 1 then Set(Cylaxisx = true); elif axis = 2 then Set(Cylaxisy = true); elif axis = 3 then Set(Cylaxisz = true); end if; end proc: Maplets[Display]( QSurfaces );