!! File: TargetPractice ! IDE - 1995 ! BDH - 2001 ! MIT - 02/01/2002 - © Hubert Hohn PUBLIC PCflag,Mac5flag,M68Kflag,Unixflag,xmax,ymax PUBLIC toolLft,toolRgt,toolBas,toolTop,toolhdr,toolHgt,toolWid ! tool boundaries PUBLIC winLft,winRgt,winBas,winTop,winHgt,winWid ! window PUBLIC workLft,workRgt,workBas,workTop,workMidx ! work area PUBLIC qLft,qRgt,qBas,qTop PUBLIC infLft,infRgt,infBas,infTop PUBLIC black,drkgry,drkmid,midgry,litmid,litgry,white PUBLIC red,yellow,green,cyan,blue,magenta,colorscheme PUBLIC planeclr,gridclr,rimclr,axisclr,axislabelclr,titleclr,rightsclr PUBLIC numberlineclr,slotdrkclr,slotlgtclr,slideclr,headerclr,btnclr PUBLIC largefonts PUBLIC title$,SLUMode LET toolHgt= 520 LET toolWid= 720 LET window$= "MIT Visual Math Project" LET colorscheme= 0 LET title$= "Target Practice" SUB ThisProgram CALL TargetPractice END SUB !! ---------------------------------------------------------- !! --- Start Unix Header and Subs --- !LET Unixflag= 1 !ASK PIXELS winWid,winHgt !LET winLft= 0 !LET winTop= 0 !LET winRgt= winWid-1 !LET winBas= winHgt-1 !SET WINDOW 0,winRgt,winBas,0 !CALL Palette ! !CALL ToolPanel !CALL ThisProgram !CLEAR ! !END !EXTERNAL ! !MODULE UnixParts ! SHARE CharWidth ! ! SUB SetTextFont(font,size,style$) ! LET font$= "-adobe-courier-" ! IF style$= "normal" then ! LET style$= "medium-r-normal--" ! ELSE ! LET style$= "bold-r-normal--" ! END IF ! IF size=9 then ! LET size$= str$(10) ! ELSE ! LET size$= str$(size) ! END IF ! LET test= SetFont(font$&style$&size$&"*") ! ! IF size=9 then ! LET CharWidth= 6 ! ELSE IF size=12 then ! numeric output - axis labels ! LET CharWidth= 7 ! ELSE IF size=14 then ! rare ! LET CharWidth= 8 ! ELSE IF size=18 then ! rare ! LET CharWidth= 10 ! END IF ! END SUB ! ! SUB StringWidth(sw$,sl) ! string width in pixels ! ! LET sl= StrWidth(sw$) ! LET chars= len(sw$) ! LET sl = chars*CharWidth ! END SUB ! ! SUB SetLineWeight(wgt) ! ! CALL PenSize(wgt,wgt) ! END SUB ! ! SUB BoxDisk(Lft,Rgt,Bas,Top) ! CALL Fill_Circle(Lft,Rgt,Bas,Top) ! END SUB !END MODULE ! --------------------------------------------------------- !! --- Start Mac TB4 Header and Subs --- !LET M68Kflag = 1 !LIBRARY "MacTools*" !ASK PIXELS winWid,winHgt !LET winLft= 0 !LET winTop= 0 !LET winRgt= winWid-1 !LET winBas= winHgt-1 !SET WINDOW 0,winRgt,winBas,0 !CALL Palette !CLEAR ! !CALL ToolPanel !CALL ThisProgram !CLEAR ! !END !EXTERNAL ! !MODULE Mac4Parts ! SUB SetTextFont(font,size,style$) ! CALL MacTextFont(font) ! CALL MacTextSize(size) ! CALL MacTextFace(style$) ! END SUB ! ! SUB StringWidth(sw$,sl) ! DECLARE DEF MacStringWidth ! LET sl= MacStringWidth(sw$) ! END SUB ! ! SUB SetLineWeight(wgt) ! CALL MacPenSize(wgt,wgt) ! END SUB ! ! SUB BoxDisk(Lft,Rgt,Bas,Top) ! CALL MacPaintOval(Lft,Rgt,Bas,Top) ! END SUB !END MODULE ! ! --- End Mac4 Header and Subs --- ! --------------------------------------------------------- ! ! --- Start Cross-Platform TB5 header and subs --- ! LIBRARY "c:\TB Silver\TBLibs\TrueCtrl.trc" ! windows !LIBRARY "c:\TBGold\TBLibs\TrueCtrl.trc" ! windows !LIBRARY "c:\program files\TB Gold\TBLIbs\TrueCtrl.trc" ! windows !LIBRARY ":TBLibs:TrueCtrl.trc" ! macintosh PUBLIC old_priority LET new_priority= 19200 ! 9600 is default CALL Priority(new_priority,old_priority) PUBLIC WinID DECLARE PUBLIC OBJM_SET,OBJM_SYSINFO LET winHgt= toolHgt LET winWid= toolWid LET WinID = 0 DIM values(1) CALL TC_Init CALL Object(OBJM_SYSINFO,WinID,"MACHINE",system$,values()) IF system$="MAC" then LET Mac5flag= 1 ELSE IF system$="WIN32" then LET PCflag = 1 END IF CALL TC_SetUnitsToPixels ! 5.1 and up needs this CALL TC_GetScreenSize(scrnLft,scrnRgt,scrnBas,scrnTop) LET winLft= int((scrnRgt-scrnLft-winWid)/2) LET winRgt= winLft+winWid-1 LET winTop= int((scrnBas-scrnTop-winHgt)/2) + 10 LET winBas= winTop+winHgt-1 CALL TC_Win_Create (WinID,"TITLE",winLft,winRgt,winBas,winTop) LET values(1)= 2 CALL Object(OBJM_SET, WinID, "TYPE", "", values()) IF PCflag=1 then LET values(1)= 1 CALL Object(OBJM_SET, WinID, "SOLID MIX", "", values()) END IF LET values(1)= 0 CALL TC_SetRect(WinID,winLft,winRgt,winBas,winTop) CALL TC_Win_SetTitle (WinID,window$) CALL TC_Show(WinID) SET MODE "COLORSTANDARD" ASK PIXELS winWid,winHgt LET winLft= 0 LET winTop= 0 LET winRgt= winWid-1 LET winBas= winHgt-1 SET WINDOW 0,winRgt,winBas,0 CALL Palette IF PCflag=1 then LET values(1)= 0 CALL Object(OBJM_SET, WinID, "SOLID MIX", "", values()) CALL TC_Win_RealizePalette(WinID) ! PC needs this in 5.1 CALL TC_Win_SetFont(WinID,"arial",9,"plain") CALL StringWidth("0",sw) IF sw>7 then LET largefonts=1 else LET largefonts=0 END IF CALL TC_Win_Switch(WinID) CALL ToolPanel CALL ThisProgram CLEAR CALL TC_CleanUp END EXTERNAL MODULE TB5Parts SUB StringWidth(sw$,sl) DECLARE PUBLIC WinID LET sl= StrWidth(WinID,sw$) END SUB SUB SetLineWeight(wgt) DECLARE PUBLIC OBJM_SET DECLARE PUBLIC WinID DIM values(1) LET values(1)= wgt CALL Object(OBJM_SET,WinID, "WIDTH", "", values()) END SUB SUB SetTextFont(font,size,style$) DECLARE PUBLIC WinID,Mac5flag,PCflag,largefonts IF Mac5flag=1 then SELECT CASE font CASE 4 LET font$= "Courier" CASE 16 LET font$= "Times" CASE else LET font$= "Geneva" END SELECT ELSE IF PCflag=1 then IF largefonts=1 then IF size<12 then LET size= 7 ELSE IF size=14 then LET size= 10 ELSE IF size=18 then LET size= 12 ELSE IF size=24 then LET size= 14 ELSE IF size=12 then LET size= 8 ELSE LET size= round(72/96 * size * .8) END IF ELSE IF size<12 then LET size= 8 ELSE IF size=14 then LET size= 12 ELSE IF size=12 then LET size= 9 ELSE IF size=18 then LET size= 14 ELSE IF size=24 then LET size= 18 ELSE LET size= round(72/96 * size) END IF END IF SELECT CASE font CASE 4 LET font$= "Courier New" CASE 16 LET font$= "Times New Roman" CASE else LET font$= "Arial" END SELECT END IF IF style$= "normal" then LET style$= "plain" CALL TC_Win_SetFont(WinID,font$,size,style$) END SUB SUB BoxDisk(Lft,Rgt,Bas,Top) BOX DISK Lft,Rgt,Bas,Top END SUB END MODULE ! --- End TB5 header and subs --- !! --------------------------------------------------------- MODULE Interface DECLARE PUBLIC PCflag,M68kflag,Mac5flag,Unixflag,xmax,ymax DECLARE PUBLIC toolLft,toolRgt,toolBas,toolTop,toolhdr,toolHgt,toolWid ! tool boundaries DECLARE PUBLIC winLft,winRgt,winBas,winTop,winHgt,winWid ! window DECLARE PUBLIC workLft,workRgt,workBas,workTop,workMidx ! work area DECLARE PUBLIC qLft,qRgt,qBas,qTop DECLARE PUBLIC infLft,infRgt,infBas,infTop DECLARE PUBLIC iLft,iRgt,iBas,iTop DECLARE PUBLIC black,drkgry,drkmid,midgry,litmid,litgry,white DECLARE PUBLIC red,yellow,green,cyan,blue,magenta,colorscheme DECLARE PUBLIC planeclr,gridclr,rimclr,axisclr,axislabelclr,titleclr,rightsclr DECLARE PUBLIC numberlineclr,slotdrkclr,slotlgtclr,slideclr,headerclr,btnclr DECLARE PUBLIC largefonts DECLARE PUBLIC title$ SUB ToolPanel LET toolLft= int((winWid-toolWid)/2) LET toolRgt= toolLft+toolWid-1 LET toolTop= int((winHgt-toolHgt)/2) LET toolBas= toolTop+toolHgt-1 LET toolhdr= toolTop+27 LET toolmid= int((toolLft+toolRgt)/2) IF M68kFlag=1 or Unixflag=1 then ASK PIXELS xpix,ypix LET xmax= xpix-1 LET ymax= ypix-1 SET COLOR drkgry BOX AREA 0,xmax,ymax,0 CALL Rim(toolLft,toolRgt,toolBas,toolTop,0) ELSE SET COLOR white BOX LINES toollft,toolrgt,toolbas,tooltop END IF CALL HRule(toollft+4,toolrgt-4,toolhdr,3,0) LET workLft = toolLft+ 5 LET workRgt = toolRgt- 5 LET workTop = toolhdr+ 2 LET workBas = toolBas- 5 LET workMidx= int((workLft+workRgt)/2) CALL SetTextFont(1,9,"bold") LET qLft= toolRgt-40 LET qTop= toolTop+5 LET qRgt= qLft+35 LET qBas= qTop+17 CALL DrawButton(qLft,qRgt,qBas,qTop,5,"Quit") LET infLft= qLft - 40 ! info button LET infTop= qTop LET infRgt= infLft + 35 LET infBas= qBas CALL DrawButton(infLft,infRgt,infBas,infTop,5,"Help") CALL SetToolTitle(title$) ! LET cl= toolLft + 30 ! LET cr= toolRgt - 30 ! LET cy= workBas+2 ! SET COLOR white ! PLOT cl,cy-1; cr,cy-1 ! SET COLOR black ! PLOT cl+1,cy; cr+1,cy !CALL CopyRight(workLft+30,cy+9,rightsclr) END SUB SUB SetToolTitle(txt$) !BOX CLEAR qRgt+5,infLft-5,toolhdr-5,toolTop+5 CALL SetTextFont(1,12,"bold") CALL PlotTextLJ(worklft+30,toolhdr-10,txt$,titleclr) END SUB SUB CopyRightHC(hp,yp,c) ! 142 pixels LET x= hp LET y= yp-7 LET x= x+35 SET COLOR c FOR v= 1 to 7 READ data$ FOR h= 1 to 107 LET n$= data$(h:h) IF n$= "1" then PLOT x+h,y+v-1 END IF NEXT h NEXT v DATA 100001000010001000000010000000000000000000000000000001000000001000001110000000000000000100000000000000010010 DATA 100001000010001000000010000000000000000000000000000001000000001000010001000000000000000100000000000000010010 DATA 100001000010001001110010110010110000000111010110001111000000001000010000000111010110011111010001001110010010 DATA 111111000011111010001011001011001000001001011001010001000000001000010000001001011001000100010001010001010010 DATA 100001000010001010001010001010001000010001010001010001000000001000010000010001010001000100010101011110010010 DATA 100001000010001010001010001010001000010011010001010001000010001000010001010011010001000100011011010000010010 DATA 100001010010001001110010001010001000001101010001001111000001110010001110001101010001000010010001001111010010 END SUB SUB CopyRightH(hp,yp,c) ! Width = 87 pixels ! Width = 131 pixels LET x= hp LET y= yp-7 SET COLOR c FOR v= 0 to 7 READ data$ FOR h= 0 to 30 LET n$= data$(h:h) IF n$= "1" then PLOT x+h,y+v END IF NEXT h NEXT v DATA 00111100000000000000000000000 DATA 01000010000110001100011000100 DATA 10011001001001010010100101100 DATA 10100001000001010010100100100 DATA 10100001000010010010100100100 DATA 10011001000100010010100100100 DATA 01000010001111001100011001110 DATA 00111100000000000000000000000 LET x= x+35 FOR v= 1 to 7 READ data$ FOR h= 1 to 33 LET n$= data$(h:h) IF n$= "1" then PLOT x+h,y+v-1 END IF NEXT h NEXT v ! midpoint - 34 DATA 10001000010001000000010000000000 DATA 10001000010001000000010000000000 DATA 10001000010001001110010110010110 DATA 11111000011111010001011001011001 DATA 10001000010001010001010001010001 DATA 10001000010001010001010001010001 DATA 10001010010001001110010001010001 END SUB SUB CopyRightHM(hp,yp,c) ! Width = 87 pixels ! Width = 131 pixels LET x= hp LET y= yp-7 LET x= x+36 FOR v= 1 to 7 READ data$ FOR h= 1 to 68 LET n$= data$(h:h) IF n$= "1" then PLOT x+h,y+v-1 END IF NEXT h NEXT v DATA 00000000000000001000010001000010001001001001000000000000 DATA 00000000000000001000010001000011011000001001000000000000 DATA 00111010110001111000010001000010101001001001000111001011 DATA 01001011001010001000011111000010101001001001001000101100 DATA 10001010001010001000010001000010001001001001001111101000 DATA 10011010001010001000010001000010001001001001001000001000 DATA 01101010001001111000010001010010001001001001000111001000 END SUB SUB Rim(l,r,b,t,c) SET COLOR 0 BOX AREA l,r,b,t SET COLOR black PLOT l+3,b-3; l+3,t+3; r-3,t+3 PLOT l+1,b; r,b; r,t+1 SET COLOR white PLOT l+4,b-3; r-3,b-3; r-3,t+4 PLOT l,b-1; l,t; r-1,t SET COLOR c BOX AREA l+5,r-5,b-5,t+5 END SUB SUB HRule(l,r,b,h,c) LET t= b-h SET COLOR c BOX AREA l,r,b,t SET COLOR black PLOT l,b; r,b SET COLOR white PLOT l,t; r,t END SUB ! SUB Panel(l,r,b,t,c) ! SET COLOR black ! BOX AREA l,r,b,t ! SET COLOR litmid ! BOX LINES l,r,b,t ! END SUB SUB Panel(l,r,b,t,c) SET COLOR 0 BOX AREA l,r,b,t SET COLOR black PLOT l,b; l,t; r,t SET COLOR litgry PLOT l,b; r,b; r,t SET COLOR c BOX AREA l+1,r-1,b-1,t+1 END SUB ! SUB Panel(l,r,b,t,c) ! SET COLOR 0 ! BOX AREA l,r,b,t ! SET COLOR white ! PLOT l+1,b; r,b; r,t+1 ! PLOT l+3,b-3; l+3,t+3; r-3,t+3 ! SET COLOR black ! PLOT l+4,b-3; r-3,b-3; r-3,t+4 ! PLOT l,b-1; l,t; r-1,t ! SET COLOR c ! BOX AREA l+5,r-5,b-5,t+5 ! END SUB ! ---- Color Subs ----- SUB Palette ! LET b0= 0 ! LET b1= 0.30 ! LET b2= 0.40 ! LET b3= 0.60 ! LET b4= 0.70 ! LET b5= 1.00 LET b0= 0 LET b1= 0.20 LET b2= 0.40 LET b3= 0.60 LET b4= 0.80 LET b5= 1.00 IF colorscheme=0 then !SET COLOR MIX(0) b0,.3,.5 ! Back SET COLOR MIX(0) 0,.2,.4 ! Back ELSE SET COLOR MIX(0) b4,b4,b4 ! Back ! SET COLOR MIX(0) .7,.7,.7 ! Back ! SET COLOR MIX(0) .7,.7,.8 ! Back END IF SET BACK 0 SET COLOR 0 CLEAR SET COLOR MIX( 1) b0,b0,b0 ! black SET COLOR MIX( 2) b1,b1,b1 ! dark gray SET COLOR MIX( 3) b2,b2,b2 ! dark mid SET COLOR MIX( 4) b3,b3,b3 ! lite mid SET COLOR MIX( 5) b4,b4,b4 ! lite gray SET COLOR MIX( 6) b5,b5,b5 ! white IF colorscheme=0 then ! blue tool SET COLOR MIX( 7) b5,b1,b1 ! red SET COLOR MIX( 8) b5,b5,b0 ! yellow SET COLOR MIX( 9) b1,b5,b0 ! green SET COLOR MIX(10) b0,b5,b5 ! cyan SET COLOR MIX(11) b4,b3,b5 ! blue SET COLOR MIX(12) b5,b0,b5 ! magenta ELSE ! white tool SET COLOR MIX( 7) b4,b1,b0 ! red SET COLOR MIX( 8) b4,b3,b0 ! yellow SET COLOR MIX( 9) b0,b3,b0 ! green SET COLOR MIX(10) b0,b3,b3 ! cyan SET COLOR MIX(11) b1,b0,b4 ! blue SET COLOR MIX(12) b3,b0,b3 ! magenta END IF FOR i= 1 to 12 SET COLOR i NEXT i LET black = 1 LET drkgry = 2 LET drkmid = 3 LET midgry = 3 LET litmid = 4 LET litgry = 5 LET white = 6 LET red = 7 LET yellow = 8 LET green = 9 LET cyan = 10 LET blue = 11 LET magenta= 12 IF colorscheme=0 then LET planeclr= black LET gridclr = drkgry LET rimclr = litmid LET axisclr = litgry LET axislabelclr = white LET titleclr = white LET rightsclr = litgry LET numberlineclr= litgry LET slotdrkclr= black LET slotlgtclr= litmid LET headerclr = blue LET btnclr= drkgry ELSE LET planeclr= white LET gridclr = litgry LET rimclr = drkmid LET axisclr = drkgry LET axislabelclr= black LET titleclr= black LET rightsclr= drkgry LET numberlineclr= drkgry LET slotdrkclr= drkgry LET slotlgtclr= white LET headerclr= blue LET btnclr= drkmid END IF END SUB SUB PaletteTest ASK PIXELS xpix,ypix LET xmax= xpix-1 LET ymax= ypix-1 FOR i= 1 to 12 SET COLOR i LET x= workLft + 100 + i*40 BOX AREA x,x+19,workBas,workTop NEXT i GET POINT: mx,my DATA 0,.35,.40,.54,.67,.80,1 END SUB END MODULE MODULE Help DECLARE PUBLIC black,drkgry,drkmid,midgry,litmid,litgry,white DECLARE PUBLIC red,yellow,green,cyan,blue,magenta,colorscheme DECLARE PUBLIC planeclr,gridclr,rimclr,axisclr,axislabelclr,titleclr DECLARE PUBLIC workLft,workRgt,workBas,workTop,workMidx ! work area DECLARE PUBLIC M68kFlag DECLARE PUBLIC toolLft,toolRgt,toolBas,toolTop,toolhdr DECLARE PUBLIC SLUmode SUB InfoPage(Info$()) LET iLft= workLft LET iRgt= workRgt-1 LET iBas= workBas-2 LET iTop= workTop SET COLOR white BOX AREA iLft,iRgt,iBas,iTop IF SLUmode=0 then CALL CopyRightH(iLft+30,iBas-3,c) ELSE CALL CopyRightHC(iLft+30,iBas-3,c) END IF LET marg1 = iLft+30 LET marg2 = marg1+24 LET txtLft= marg1 LET txtRgt= iRgt-30 LET Basln = iTop+50 LET lnspc = 15 LET hlfspc= 8 LET clr = black SET COLOR clr CALL SetTextFont(1,12,"bold") LET linecount= ubound(Info$) FOR i= 1 to linecount LET txt$ = Info$(i) LET first$= txt$(1:1) IF first$= "x" or first$= "y" then LET txtLft= marg2 CALL SuperSubScriptLJ(txtLft,Basln,txt$,clr) LET Basln = Basln+lnspc ELSE IF first$= "-" then LET txtLft= marg2 LET txt$(1:1)= "" ELSE LET txtLft= marg1 END IF CALL Paragraph(txtLft,txtRgt,Basln,lnspc,txt$) END IF LET Basln= Basln+hlfspc NEXT i LET midpage= int((iLft+iRgt)/2) CALL SetTextFont(1,9,"bold") LET cbtnLft= midpage-35 LET cbtnRgt= midpage+35 LET cbtnBas= iBas-10 LET cbtnTop= cbtnBas-17 CALL DrawButton(cbtnLft,cbtnRgt,cbtnBas,cbtnTop,5,"Continue") CALL MouseDown(mx,my,ms) IF mx>cbtnLft and mxcbtnTop and myWidth or Rgtpnt=endofpar LET line$ = trim$(para$(Lftpnt:Rgtpnt)) CALL PlotTextLJ(txtLft,Basln,line$,clr) LET Basln = Basln + lnspc LET Lftpnt= min(Rgtpnt+1,endofpar) LOOP until Lftpnt=endofpar END SUB END SUB SUB Equation(Bas) LET Midx= (iLft+iRgt)/2 CALL SetTextFont(1,12,"bold") LET length$= "A = k ert" CALL StringWidth(length$,sl) LET left= txtLft+30 LET nt$ = "A = k e^[rt]" CALL SuperSubScriptLJ(left,Bas,nt$,clr) END SUB END MODULE ! --- external support libraries --- ! ---- graphing methods SUB ClipVert(Bas,Top,x1,y1,x2,y2) ! this handles vertical clipping for pixel values LOCAL xd,yd,xm LET yd= y2-y1 LET xd= x2-x1 LET m = xd/yd IF y1>Bas then ! Clip Bas LET x1= x2 + m*(Bas-y2) LET y1= Bas ELSE IF y1Bas then ! Clip Bas LET x2= x1 + m*(Bas-y1) LET y2= Bas ELSE IF y20 then LET lines= abs(int(fWid/xstp)) FOR i= 0 to lines LET n = fLft + i*xstp LET wx= wndx(n) PLOT wx,wTop+1; wx,wBas-1 NEXT i END IF IF ystp>0 then LET lines= abs(int(fHgt/ystp)) IF fHgt<0 then LET first= fTop else LET first= fBas FOR i= 0 to lines LET n = first + i*ystp LET wy= wndy(n) PLOT wLft+1,wy; wRgt-1,wy NEXT i END IF SET COLOR rimclr BOX LINES wLft,wRgt,wBas,wTop END SUB SUB ZeroAxes(wx0,wy0,wLft,wRgt,wBas,wTop,dir,aclr) IF wy0>=wTop and wy0<=wBas then SET COLOR aclr PLOT wLft,wy0; wRgt+3,wy0 DRAW arrow3 with rotate(0) * shift(wRgt+3,wy0) END IF IF wx0>=wLft and wx0<=wRgt then SET COLOR aclr IF dir=1 then PLOT wx0,wBas; wx0,wTop-3 DRAW arrow3 with rotate(-pi/2) * shift(wx0,wTop-3) ELSE PLOT wx0,wBas+3; wx0,wTop DRAW arrow3 with rotate(pi/2) * shift(wx0,wBas+3) END IF END IF END SUB SUB EdgesVrt(wAxis,wBas,wTop,fBas,fTop,first,stp1,stp2,nstp,dir) BOX CLEAR wAxis-35,wAxis-5,wBas+4,wTop-4 BOX CLEAR wAxis- 4,wAxis ,wBas ,wTop CALL VNumberLineTiks(wBas,wTop,wAxis,fBas,fTop,first,stp1,stp2,nstp,dir) CALL VNumberLineLabels(wBas,wTop,wAxis,fBas,fTop,first,stp1,stp2,nstp,dir) END SUB SUB EdgesVrtPi(wAxis,wBas,wTop,fBas,fTop,first,stp1,stp2,nstp,dir) BOX CLEAR wAxis-35,wAxis-5,wBas+4,wTop-4 BOX CLEAR wAxis- 4,wAxis ,wBas ,wTop CALL VNumberLineTiks(wBas,wTop,wAxis,fBas,fTop,first,stp1,stp2,nstp,dir) CALL VNumberLineLabelsPi(wBas,wTop,wAxis,fBas,fTop,first,stp1,stp2,nstp,dir) END SUB SUB EdgesHrz(wLft,wRgt,wBas,fLft,fRgt,first,stp1,stp2,nstp,dir) IF dir=1 then BOX CLEAR wLft,wRgt,wBas+1,wBas+4 BOX CLEAR wLft-15,wRgt+20,wBas+4,wBas+15 END IF CALL HNumberLineTiks(wLft,wRgt,wBas,fLft,fRgt,first,stp1,stp2,nstp,dir) CALL HNumberLineLabels(wLft,wRgt,wBas,fLft,fRgt,first,stp1,stp2,nstp,dir) END SUB SUB EdgesHrzPi(wLft,wRgt,wBas,fLft,fRgt,first,stp1,stp2,nstp,dir) BOX CLEAR wLft,wRgt,wBas+1,wBas+4 BOX CLEAR wLft-15,wRgt+20,wBas+4,wBas+15 CALL HNumberLineTiks(wLft,wRgt,wBas,fLft,fRgt,first,stp1,stp2,nstp,dir) CALL HNumberLineLabelsPi(wLft,wRgt,wBas,fLft,fRgt,first,stp1,stp2,nstp,dir) END SUB END MODULE MODULE MathGraphs DECLARE PUBLIC black,drkgry,drkmid,midgry,litmid,litgry,white DECLARE PUBLIC red,yellow,green,cyan,blue,magenta,colorscheme DECLARE PUBLIC planeclr,gridclr,rimclr,axisclr,axislabelclr,titleclr DECLARE PUBLIC numberlineclr,slotdrkclr,slotlgtclr,slideclr SUB VNumberLineTiks(wBas,wTop,wAxis,fBas,fTop,first,stp1,stp2,nstp,ndir) LOCAL fy,wHgt,fHgt,ticks,i,n,wy,s,n$,sl,form$,places DEF wndy(fy)= wBas - wHgt/fHgt*(fy-fBas) LET wHgt= wBas-wTop LET fHgt= fTop-fBas SET COLOR numberlineclr IF stp1>0 or stp2>0 then PLOT wAxis,wBas; wAxis,wTop ! Y boundary axis END IF IF stp1>0 then ! short ticks LET ticks= int(abs(fHgt/stp1)) FOR i= 0 to ticks LET n = first + i*stp1 LET wy= wndy(n) IF ndir=-1 then PLOT wAxis,wy; wAxis-2,wy ELSE IF ndir=1 then PLOT wAxis,wy; wAxis+2,wy END IF NEXT i END IF IF stp2>0 then ! long ticks LET ticks= int(abs(fHgt/stp2)) FOR i= 0 to ticks LET n = first + i*stp2 LET wy= wndy(n) IF ndir=-1 then PLOT wAxis,wy; wAxis-4,wy ELSE IF ndir=1 then PLOT wAxis,wy; wAxis+4,wy END IF NEXT i END IF END SUB SUB VNumberLineLabels(wBas,wTop,wAxis,fBas,fTop,first,stp1,stp2,nstp,ndir) LOCAL fx,wHgt,fHgt,ticks,i,n,wy,s,n$,sl,form$,places DEF wndy(fy)= wBas - wHgt/fHgt*(fy-fBas) LET wHgt= wBas-wTop LET fHgt= fTop-fBas IF ndir=-1 then IF stp1=0 and stp2=0 then LET rgtln= wAxis-3 ELSE LET rgtln= wAxis-6 END IF ELSE IF ndir=1 then IF stp1=0 and stp2=0 then LET rgtln= wAxis+7 ELSE LET rgtln= wAxis+11 END IF END IF SET COLOR numberlineclr IF stp1>0 or stp2>0 then PLOT wAxis,wBas; wAxis,wTop ! Y boundary axis END IF IF nstp>0 then ! labels IF int(nstp)= nstp then ! build format LET form$= "------#" ELSE LET places= abs(int(log10(nstp))) LET form$ = "--%." & repeat$("#",places) END IF CALL SetTextFont(1,9,"normal") LET ticks= int(abs(fHgt/nstp)) FOR i= 0 to ticks LET n = first + i*nstp LET n = round(n,7) LET wy= wndy(n) LET s = sgn(n) IF s=0 then LET n$= "0" !CALL SetTextFont(1,9,"bold") ELSE LET n$= trim$(using$(form$,n)) !CALL SetTextFont(1,9,"normal") END IF CALL PlotTextRJ(rgtln,wy+3,n$,numberlineclr) NEXT i END IF END SUB SUB VNumberLineLabelsPi(wBas,wTop,wAxis,fBas,fTop,first,stp1,stp2,nstp,ndir) LOCAL fx,wHgt,fHgt,ticks,i,n,wy,s,n$,sl,form$,places DEF pi$= "w" DEF wndy(fy)= wBas - wHgt/fHgt*(fy-fBas) LET wHgt= wBas-wTop LET fHgt= fTop-fBas IF ndir=-1 then IF stp1=0 and stp2=0 then LET rgtln= wAxis-3 ELSE LET rgtln= wAxis-6 END IF ELSE IF ndir=1 then IF stp1=0 and stp2=0 then LET rgtln= wAxis+7 ELSE LET rgtln= wAxis+11 END IF END IF SET COLOR numberlineclr IF stp1>0 or stp2>0 then PLOT wAxis,wBas; wAxis,wTop ! Y boundary axis END IF IF nstp>0 then ! labels IF int(nstp)= nstp then ! build format LET form$= "------#" ELSE LET places= abs(int(log10(nstp))) LET form$ = "--%." & repeat$("#",places) END IF CALL SetTextFont(1,9,"normal") LET ticks= abs(int(fHgt/nstp)) FOR i= 0 to ticks LET n = first + i*nstp LET n = round(n,5) LET wy = wndy(n) LET abn= abs(n) LET s = sgn(n) CALL SetTextFont(1,9,"normal") SELECT CASE abn CASE 0 LET n$= "0" !CALL SetTextFont(1,9,"bold") CASE 0.5 LET n$= pi$ & "/2" CASE 1 LET n$= pi$ CASE 1.5 LET n$= "3" & pi$ & "/2" CASE 2.5 LET n$= "5" & pi$ & "/2" CASE 3.5 LET n$= "7" & pi$ & "/2" CASE 2,3,4,5,6,7,8,9,10,11,12 LET n$= trim$(str$(abn)) & pi$ CASE else END SELECT CALL StringWidth(n$,sl) LET left= rgtln - sl LET p= pos(n$,pi$) LET l= len(n$) LET basln= wy+3 IF p=0 then CALL PlotTextRJ(rgtln,basln,n$,numberlineclr) ELSE IF p=1 then LET n$= n$(2:l) DRAW pi9 with shift(left,basln) CALL PlotTextRJ(rgtln,basln,n$,numberlineclr) ELSE IF p=l then LET n$= n$(1:l-1) CALL PlotTextLJ(left,basln,n$,numberlineclr) CALL StringWidth(n$,sw) DRAW pi9 with shift(left+sw,basln) ELSE LET Lft$= n$(1:p-1) LET Rgt$= n$(p+1:l) CALL PlotTextLJ(left,basln,Lft$,numberlineclr) CALL StringWidth(Lft$,sw) DRAW pi9 with shift(left+sw,basln) CALL PlotTextLJ(left+sw+8,basln,Rgt$,numberlineclr) END IF IF s=-1 then CALL PlotTextRJ(left-1,basln,"-",numberlineclr) NEXT i END IF END SUB SUB HNumberLineTiks(wLft,wRgt,wBas,fLft,fRgt,first,stp1,stp2,nstp,ndir) LOCAL wWid,fWid,ticks,i,n,wx DEF wndx(fx)= wLft + wWid/fWid*(fx-fLft) ! function to window LET wWid= wRgt-wLft LET fWid= fRgt-fLft SET COLOR numberlineclr IF stp1<>0 or stp2<>0 then PLOT wLft,wBas; wRgt,wBas ! axis line if ticks are on END IF IF stp1>0 then ! short ticks LET ticks= abs(int(fWid/stp1)) FOR i= 0 to ticks LET n = first + i*stp1 LET wx= wndx(n) IF ndir=-1 then PLOT wx,wBas; wx,wBas-2 ELSE IF ndir=1 then PLOT wx,wBas; wx,wBas+2 END IF NEXT i END IF IF stp2>0 then ! long ticks LET ticks= abs(int(fWid/stp2)) FOR i= 0 to ticks LET n = first + i*stp2 LET wx= wndx(n) IF ndir=-1 then PLOT wx,wBas; wx,wBas-4 ELSE IF ndir=1 then PLOT wx,wBas; wx,wBas+4 END IF NEXT i END IF END SUB SUB HNumberLineLabels(wLft,wRgt,wBas,fLft,fRgt,first,stp1,stp2,nstp,ndir) LOCAL wWid,fWid,places,form$,ticks,i,n,wx,s,n$,sl,left,basln DEF wndx(fx)= wLft + wWid/fWid*(fx-fLft) ! function to window LET wWid= wRgt-wLft LET fWid= fRgt-fLft IF ndir=-1 then IF stp1=0 and stp2=0 then LET basln= wBas-5 ELSE LET basln= wBas-7 END IF ELSE IF ndir=1 then IF stp1=0 and stp2=0 then LET basln= wBas+11 ELSE LET basln= wBas+14 END IF END IF SET COLOR numberlineclr IF stp1<>0 or stp2<>0 then PLOT wLft,wBas; wRgt,wBas ! axis line if ticks are on END IF IF nstp>0 then ! labels IF int(nstp)= nstp then ! build format LET form$ = "----#" ELSE LET places= abs(int(log10(nstp))) LET form$ = "--%." & repeat$("#",places) !LET form$ = "-%.#" END IF CALL SetTextFont(1,9,"normal") LET ticks= abs(int(fWid/nstp)) FOR i= 0 to ticks LET n = first + i*nstp LET n = round(n,5) LET wx= wndx(n) LET s = sgn(n) IF s=0 then LET n$= "0" !CALL SetTextFont(1,9,"bold") ELSE LET n$= trim$(using$(form$,abs(n))) !CALL SetTextFont(1,9,"normal") END IF CALL StringWidth(n$,sl) LET left= round(wx - sl/2) + 1 CALL PlotTextLJ(left,basln,n$,numberlineclr) IF s=-1 then CALL PlotTextRJ(left-1,basln,"-",numberlineclr) NEXT i END IF END SUB SUB HNumberLineLabelsPi(wLft,wRgt,wBas,fLft,fRgt,first,stp1,stp2,nstp,ndir) LOCAL wWid,fWid,places,form$,ticks,i,n,wx,abn,s,n$,sl,lft$,rgt$,left,basln DEF pi$= "w" DEF wndx(fx)= wLft + wWid/fWid*(fx-fLft) ! function to window LET wWid= wRgt-wLft LET fWid= fRgt-fLft IF ndir=-1 then IF stp1=0 and stp2=0 then LET basln= wBas-3 ELSE LET basln= wBas-7 END IF ELSE IF ndir=1 then IF stp1=0 and stp2=0 then LET basln= wBas+10 ELSE LET basln= wBas+14 END IF END IF SET COLOR numberlineclr IF stp1<>0 or stp2<>0 then PLOT wLft,wBas; wRgt,wBas ! axis line if ticks are on END IF IF nstp>0 then ! labels IF int(nstp)= nstp then ! build format LET form$ = "----#" ELSE LET places= abs(int(log10(nstp))) LET form$ = "--%." & repeat$("#",places) END IF CALL SetTextFont(1,9,"normal") LET ticks= abs(int(fWid/nstp)) FOR n= fLft to fRgt step nstp LET n = round(n,5) LET wx = wndx(n) LET abn= abs(n) LET s = sgn(n) CALL SetTextFont(1,9,"normal") SELECT CASE abn CASE 0 LET n$= "0" !CALL SetTextFont(1,9,"bold") CASE 0.5 LET n$= pi$ & "/2" CASE 1 LET n$= pi$ CASE 1.5 LET n$= "3" & pi$ & "/2" CASE 2.5 LET n$= "5" & pi$ & "/2" CASE 3.5 LET n$= "7" & pi$ & "/2" CASE else LET n$= trim$(str$(abn)) & pi$ !CASE else END SELECT CALL StringWidth(n$,sl) LET left= wx - sl/2 + 1 LET p= pos(n$,pi$) LET l= len(n$) IF p=0 then CALL PlotTextLJ(left,basln,n$,numberlineclr) ELSE IF p=1 then LET n$= n$(2:l) DRAW pi9 with shift(left,basln) CALL PlotTextLJ(left+8,basln,n$,numberlineclr) ELSE IF p=l then LET n$= n$(1:l-1) CALL PlotTextLJ(left,basln,n$,numberlineclr) CALL StringWidth(n$,sw) DRAW pi9 with shift(left+sw,basln) ELSE LET Lft$= n$(1:p-1) LET Rgt$= n$(p+1:l) CALL PlotTextLJ(left,basln,Lft$,numberlineclr) CALL StringWidth(Lft$,sw) DRAW pi9 with shift(left+sw,basln) CALL PlotTextLJ(left+sw+8,basln,Rgt$,numberlineclr) END IF IF s=-1 then CALL PlotTextRJ(left-1,basln,"-",numberlineclr) NEXT n END IF END SUB END MODULE ! ----- library slider subs MODULE v1slider PUBLIC v1axis,v1wLft,v1wRgt,v1wBas,v1wTop,v1wHgt,v1sBas,v1sTop PUBLIC v1fBas,v1fTop,v1First,v1STik,v1LTik,v1Label,v1fHgt,v1yf,v1y2,v1y4,v1yl PUBLIC v1fRatio,v1wRatio,v1name$,v1m,v1form$,v1clr,v1slot$ DEF v1fncy(wy)= v1fBas + v1fRatio*(v1wBas-wy) ! window to function DEF v1wndy(fy)= v1wBas + v1wRatio*(fy-v1fBas) ! function to window SUB v1SliderVariables CALL vWnd(v1axis,v1wLft,v1wRgt,v1wBas,v1wTop,v1wHgt,v1sBas,v1sTop) CALL vFnc(v1fBas,v1fTop,v1fHgt) CALL vRatios(v1fHgt,v1wHgt,v1fRatio,v1wRatio) END SUB SUB v1DrawSlider(n) BOX CLEAR v1wLft,v1wRgt,v1wBas+15,v1wTop-20 CALL SliderSlotV(v1wBas,v1wTop,v1axis) CALL SliderAxisVrt(v1wBas,v1wTop,v1axis,v1fBas,v1fTop,v1first,v1STik,v1LTik,v1Label) BOX KEEP v1axis,v1wRgt,v1sBas,v1sTop in v1slot$ CALL PlotVSliderName(v1axis,v1wBas,v1name$,v1clr) CALL v1Mark(n) END SUB SUB v1Mark(n) LET wy= v1wBas - v1wRatio*(n-v1fBas) BOX SHOW v1slot$ at v1axis,v1sBas CALL SliderKnobV(v1wRgt-5,wy) CALL PlotVSliderValue(v1axis,v1wTop-5,using$(v1form$,n*v1m),v1clr) END SUB END MODULE MODULE h1slider PUBLIC h1axis,h1wLft,h1wRgt,h1wBas,h1wTop,h1wWid,h1sLft,h1sRgt PUBLIC h1fLft,h1fRgt,h1fWid,h1First,h1STik,h1LTik,h1Label ! h1xf,h1x2,h1x4,h1xl PUBLIC h1fratio,h1wratio,h1name$,h1m,h1form$,h1clr,h1slot$ DEF h1fncx(wx)= h1fLft + h1fratio*(wx-h1wLft) ! window to function DEF h1wndx(fx)= h1wLft + h1wratio*(fx-h1fLft) ! function to window SUB h1SliderVariables CALL hWnd(h1axis,h1wLft,h1wRgt,h1wBas,h1wTop,h1wWid,h1sLft,h1sRgt) CALL hFnc(h1fLft,h1fRgt,h1fWid) CALL hRatios(h1fWid,h1wWid,h1fratio,h1wratio) END SUB SUB h1DrawSlider(n) BOX CLEAR h1wLft-20,h1wRgt+35,h1wBas,h1wTop CALL SliderSlotH(h1wLft,h1wRgt,h1axis) CALL SliderAxisHrz(h1wLft,h1wRgt,h1axis,h1fLft,h1fRgt,h1First,h1STik,h1LTik,h1Label) BOX KEEP h1sLft,h1sRgt,h1wBas,h1axis in h1slot$ CALL PlotSliderName(h1wLft,h1wBas,h1name$,h1clr) CALL h1Mark(n) END SUB SUB h1Mark(n) LET wx= h1wLft + h1wratio*(n-h1fLft) BOX SHOW h1slot$ at h1sLft,h1wBas CALL SliderKnob(wx,h1wBas-5) CALL PlotSliderValue(h1wRgt,h1wBas,using$(h1form$,n*h1m),h1clr) END SUB END MODULE MODULE hSlider2 PUBLIC h2axis,h2wLft,h2wRgt,h2wBas,h2wTop,h2wWid,h2sLft,h2sRgt PUBLIC h2fLft,h2fRgt,h2First,h2STik,h2LTik,h2Label,h2fWid,h2xf,h2x2,h2x4,h2xl PUBLIC h2fratio,h2wratio,h2name$,h2m,h2form$,h2clr PUBLIC hslot2$ DEF h2Fncx(wx)= h2fLft + h2fratio*(wx-h2wLft) ! window to function DEF h2Wndx(fx)= h2wLft + h2wratio*(fx-h2fLft) ! function to window SUB hSlider2Variables CALL hWnd(h2axis,h2wLft,h2wRgt,h2wBas,h2wTop,h2wWid,h2sLft,h2sRgt) CALL hFnc(h2fLft,h2fRgt,h2fWid) CALL hRatios(h2fWid,h2wWid,h2fratio,h2wratio) END SUB SUB h2DrawSlider(n) BOX CLEAR h2wLft-20,h2wRgt+35,h2wBas,h2wTop CALL SliderSlotH(h2wLft,h2wRgt,h2axis) CALL SliderAxisHrz(h2wLft,h2wRgt,h2axis,h2fLft,h2fRgt,h2First,h2STik,h2LTik,h2Label) BOX KEEP h2sLft,h2sRgt,h2wBas,h2axis in hSlot2$ CALL PlotSliderName(h2wLft,h2wBas,h2name$,h2clr) CALL h2Mark(n) END SUB SUB h2Mark(n) LET wx= h2wLft + h2wratio*(n-h2fLft) BOX SHOW hSlot2$ at h2sLft,h2wBas CALL SliderKnob(wx,h2wBas-5) CALL PlotSliderValue(h2wRgt,h2wBas,using$(h2form$,n*h2m),h2clr) END SUB END MODULE MODULE hSlider3 PUBLIC h3axis,h3wLft,h3wRgt,h3wBas,h3wTop,h3wWid,h3sLft,h3sRgt PUBLIC h3fLft,h3fRgt,h3First,h3STik,h3LTik,h3Label,h3fWid,h3xf,h3x2,h3x4,h3xl PUBLIC h3fratio,h3wratio,h3name$,h3m,h3form$,h3clr PUBLIC hslot3$ DEF h3Fncx(wx)= h3fLft + h3fratio*(wx-h3wLft) ! window to function DEF h3Wndx(fx)= h3wLft + h3wratio*(fx-h3fLft) ! function to window SUB hSlider3Variables CALL hWnd(h3axis,h3wLft,h3wRgt,h3wBas,h3wTop,h3wWid,h3sLft,h3sRgt) CALL hFnc(h3fLft,h3fRgt,h3fWid) CALL hRatios(h3fWid,h3wWid,h3fratio,h3wratio) END SUB SUB h3DrawSlider(n) BOX CLEAR h3wLft-20,h3wRgt+35,h3wBas,h3wTop CALL SliderSlotH(h3wLft,h3wRgt,h3axis) CALL SliderAxisHrz(h3wLft,h3wRgt,h3axis,h3fLft,h3fRgt,h3First,h3STik,h3LTik,h3Label) BOX KEEP h3sLft,h3sRgt,h3wBas,h3axis in hSlot3$ CALL PlotSliderName(h3wLft,h3wBas,h3name$,h3clr) CALL h3Mark(n) END SUB SUB h3Mark(n) LET wx= h3wLft + h3wratio*(n-h3fLft) BOX SHOW hSlot3$ at h3sLft,h3wBas CALL SliderKnob(wx,h3wBas-5) CALL PlotSliderValue(h3wRgt,h3wBas,using$(h3form$,n*h3m),h3clr) END SUB END MODULE MODULE hSlider4 PUBLIC h4axis,h4wLft,h4wRgt,h4wBas,h4wTop,h4wWid,h4sLft,h4sRgt PUBLIC h4fLft,h4fRgt,h4First,h4STik,h4LTik,h4Label,h4fWid,h4xf,h4x2,h4x4,h4xl PUBLIC h4fratio,h4wratio,h4name$,h4m,h4form$,h4clr PUBLIC hSlot4$ DEF h4Fncx(wx)= h4fLft + h4fratio*(wx-h4wLft) ! window to function DEF h4Wndx(fx)= h4wLft + h4wratio*(fx-h4fLft) ! function to window SUB hSlider4Variables CALL hWnd(h4axis,h4wLft,h4wRgt,h4wBas,h4wTop,h4wWid,h4sLft,h4sRgt) CALL hFnc(h4fLft,h4fRgt,h4fWid) CALL hRatios(h4fWid,h4wWid,h4fratio,h4wratio) END SUB SUB h4DrawSlider(n) BOX CLEAR h4wLft-20,h4wRgt+45,h4wBas,h4wTop CALL SliderSlotH(h4wLft,h4wRgt,h4axis) CALL SliderAxisHrz(h4wLft,h4wRgt,h4axis,h4fLft,h4fRgt,h4First,h4STik,h4LTik,h4Label) BOX KEEP h4sLft,h4sRgt,h4wBas,h4axis in hSlot4$ CALL PlotSliderName(h4wLft,h4wBas,h4name$,h4clr) CALL h4Mark(n) END SUB SUB h4Mark(n) LET wx= h4wLft + h4wratio*(n-h4fLft) BOX SHOW hSlot4$ at h4sLft,h4wBas CALL SliderKnob(wx,h4wBas-5) CALL PlotSliderValue(h4wRgt,h4wBas,using$(h4form$,n*h4m),h4clr) END SUB END MODULE MODULE hSlider5 PUBLIC h5axis,h5wLft,h5wRgt,h5wBas,h5wTop,h5wWid,h5sLft,h5sRgt PUBLIC h5fLft,h5fRgt,h5First,h5STik,h5LTik,h5Label,h5fWid,h5xf,h5x2,h5x5,h5xl PUBLIC h5fratio,h5wratio,h5name$,h5m,h5form$,h5clr PUBLIC hSlot5$ DEF h5Fncx(wx)= h5fLft + h5fratio*(wx-h5wLft) ! window to function DEF hwndx5(fx)= h5wLft + h5wratio*(fx-h5fLft) ! function to window SUB hSlider5Variables CALL hWnd(h5axis,h5wLft,h5wRgt,h5wBas,h5wTop,h5wWid,h5sLft,h5sRgt) CALL hFnc(h5fLft,h5fRgt,h5fWid) CALL hRatios(h5fWid,h5wWid,h5fratio,h5wratio) END SUB SUB h5DrawSlider(n) BOX CLEAR h5wLft-20,h5wRgt+45,h5wBas,h5wTop CALL SliderSlotH(h5wLft,h5wRgt,h5axis) CALL SliderAxisHrz(h5wLft,h5wRgt,h5axis,h5fLft,h5fRgt,h5First,h5STik,h5LTik,h5Label) BOX KEEP h5sLft,h5sRgt,h5wBas,h5axis in hSlot5$ CALL PlotSliderName(h5wLft,h5wBas,h5name$,h5clr) CALL h5Mark(n) END SUB SUB h5Mark(n) LET wx= h5wLft + h5wratio*(n-h5fLft) BOX SHOW hSlot5$ at h5sLft,h5wBas CALL SliderKnob(wx,h5wBas-5) CALL PlotSliderValue(h5wRgt,h5wBas,using$(h5form$,n*h5m),h5clr) END SUB END MODULE MODULE hSlider6 PUBLIC h6axis,h6wLft,h6wRgt,h6wBas,h6wTop,h6wWid,h6sLft,h6sRgt PUBLIC h6fLft,h6fRgt,h6First,h6STik,h6LTik,h6Label,h6fWid,h6xf,h6x2,h6x6,h6xl PUBLIC h6fratio,h6wratio,h6name$,h6m,h6form$,h6clr PUBLIC hSlot6$ DEF h6Fncx(wx)= h6fLft + h6fratio*(wx-h6wLft) ! window to function DEF hwndx6(fx)= h6wLft + h6wratio*(fx-h6fLft) ! function to window SUB hSlider6Variables CALL hWnd(h6axis,h6wLft,h6wRgt,h6wBas,h6wTop,h6wWid,h6sLft,h6sRgt) CALL hFnc(h6fLft,h6fRgt,h6fWid) CALL hRatios(h6fWid,h6wWid,h6fratio,h6wratio) END SUB SUB h6DrawSlider(n) BOX CLEAR h6wLft-20,h6wRgt+46,h6wBas,h6wTop CALL SliderSlotH(h6wLft,h6wRgt,h6axis) CALL SliderAxisHrz(h6wLft,h6wRgt,h6axis,h6fLft,h6fRgt,h6First,h6STik,h6LTik,h6Label) BOX KEEP h6sLft,h6sRgt,h6wBas,h6axis in hSlot6$ CALL PlotSliderName(h6wLft,h6wBas,h6name$,h6clr) CALL h6Mark(n) END SUB SUB h6Mark(n) LET wx= h6wLft + h6wratio*(n-h6fLft) BOX SHOW hSlot6$ at h6sLft,h6wBas CALL SliderKnob(wx,h6wBas-5) CALL PlotSliderValue(h6wRgt,h6wBas,using$(h6form$,n*h6m),h6clr) END SUB END MODULE MODULE SliderParts DECLARE PUBLIC black,drkgry,drkmid,midgry,litmid,litgry,white DECLARE PUBLIC red,yellow,green,cyan,blue,magenta,colorscheme DECLARE PUBLIC planeclr,gridclr,rimclr,axisclr,axislabelclr,titleclr DECLARE PUBLIC numberlineclr,slotdrkclr,slotlgtclr,slideclr ! ------------- SUB vWnd(axx,wLft,wRgt,wBas,wTop,wHgt,sBas,sTop) LET wHgt = wBas-wTop LET wLft = axx-13 LET wRgt = axx+13 LET sBas = wBas+5 LET sTop = wTop-5 END SUB SUB vFnc(fBas,fTop,fHgt) LET fHgt = fTop-fBas END SUB SUB vRatios(vfHgt,vwHgt,vfratio,vwratio) LET vfratio= vfHgt/vwHgt LET vwratio= vwHgt/vfHgt END SUB SUB PlotVSliderName(Mid,Bas,name$,clr) CALL SetTextFont(1,12,"bold") CALL PlotTextCJ(Mid,Bas+13,name$,clr) END SUB SUB PlotVSliderValue(Mid,Bas,value$,clr) CALL SetTextFont(1,12,"bold") BOX CLEAR mid-20,mid+20,Bas,Bas-13 CALL PlotTextCJ(Mid,Bas-3,trim$(value$),clr) END SUB SUB SliderslotV(wb,wt,axis) LET sx= axis+7 BOX CLEAR axis+1,axis+14,wb+4,wt-4 SET COLOR slotdrkclr BOX AREA sx,sx-1,wb+1,wt-1 SET COLOR slotlgtclr PLOT sx+1,wb+1; sx+1,wt-1 END SUB SUB PhaseSlotV(sLft,sRgt,syb,syt,sx,sliderV$) BOX CLEAR sLft,sRgt,syb,syt SET COLOR white PLOT sx+1,syb-4; sx+1,syt+4 BOX KEEP sLft,sRgt,syb,syt in sliderV$ END SUB SUB SliderAxisVrt(wBas,wTop,wAxis,fBas,fTop,first,stp1,stp2,nstp) CALL VNumberLineTiks(wBas,wTop,wAxis,fBas,fTop,first,stp1,stp2,nstp,-1) CALL VNumberLineLabels(wBas,wTop,wAxis,fBas,fTop,first,stp1,stp2,nstp,-1) END SUB SUB SliderKnobV(kx,ky) SET COLOR black BOX AREA kx-3,kx+4,ky+3,ky-2 BOX AREA kx-4,kx+5,ky+2,ky-1 SET COLOR litmid BOX AREA kx-3,kx+3,ky+2,ky-2 BOX AREA kx-4,kx+4,ky+1,ky-1 SET COLOR white PLOT kx-4,ky; kx-6,ky SET COLOR litgry PLOT kx-3,ky-2; kx+2,ky-2 PLOT kx-3,ky+2; kx+2,ky+2 PLOT kx-4,ky-1; kx-4,ky+1 SET COLOR white PLOT kx-2,ky; kx+2,ky END SUB ! ------------- SUB hWnd(axy,wLft,wRgt,wBas,wTop,wWid,sLft,sRgt) LET wWid = wRgt-wLft LET wTop = axy-13 LET wBas = axy+13 LET sLft = wLft-5 LET sRgt = wRgt+5 END SUB SUB hFnc(fLft,fRgt,fWid) LET fWid = fRgt-fLft END SUB SUB hRatios(hfWid,hwWid,hfratio,hwratio) LET hfratio= hfWid/hwWid LET hwratio= hwWid/hfWid END SUB SUB PlotSliderName(Lft,Bas,name$,clr) CALL SetTextFont(1,12,"bold") CALL PlotTextRJ(Lft-5,Bas-3,name$,clr) END SUB SUB PlotSliderValue(Rgt,Bas,value$,clr) CALL SetTextFont(1,12,"bold") BOX CLEAR Rgt+3,Rgt+50,Bas,Bas-13 CALL PlotTextLJ(Rgt+7,Bas-3,trim$(value$),clr) END SUB SUB SliderSlotH(wl,wr,axis) LET sy= axis+7 BOX CLEAR wl-4,wr+4,axis+14,axis+1 SET COLOR slotdrkclr BOX AREA wl-1,wr+1,sy,sy-1 SET COLOR slotlgtclr PLOT wl-1,sy+1; wr+1,sy+1 END SUB SUB SliderAxisHrz(wLft,wRgt,wBas,fLft,fRgt,first,stp1,stp2,nstp) CALL HNumberLineTiks(wLft,wRgt,wBas,fLft,fRgt,first,stp1,stp2,nstp,-1) CALL HNumberLineLabels(wLft,wRgt,wBas,fLft,fRgt,first,stp1,stp2,nstp,-1) END SUB SUB SliderAxisHrzPi(wLft,wRgt,axy,fLft,fRgt,first,stp1,stp2,nstp) LOCAL wWid,fWid,n,n$,left,sl,places,form$ DEF pi$= "w" DEF fncx(wx)= fLft + fWid*(wx-wLft)/wWid ! win to func DEF wndx(fx)= round(wLft + wWid*(fx-fLft)/fWid) ! func to win LET wWid= wRgt-wLft LET fWid= fRgt-fLft CALL HNumberLineTiks(wLft,wRgt,wBas,fLft,fRgt,first,stp1,stp2,nstp,-1) CALL HNumberLineLabelsPi(wLft,wRgt,wBas,fLft,fRgt,first,stp1,stp2,nstp,-1) END SUB SUB SliderKnob(kx,ky) IF colorscheme=1 then SET COLOR black BOX AREA kx-2,kx+3,ky+4,ky-3 BOX AREA kx-1,kx+2,ky+5,ky-3 SET COLOR drkmid BOX AREA kx-2,kx+2,ky+2,ky-3 BOX AREA kx-1,kx+1,ky+4,ky-4 SET COLOR litgry PLOT kx-2,ky-3; kx-2,ky+2 PLOT kx+2,ky-3; kx+2,ky+2 PLOT kx-1,ky-4; kx+1,ky-4 SET COLOR black PLOT kx,ky-4; kx,ky-6 SET COLOR white PLOT kx,ky+3; kx,ky-2 ELSE SET COLOR black BOX AREA kx-2,kx+3,ky+4,ky-3 BOX AREA kx-1,kx+2,ky+5,ky-4 SET COLOR litmid BOX AREA kx-2,kx+2,ky+2,ky-3 BOX AREA kx-1,kx+1,ky+4,ky-4 SET COLOR litgry PLOT kx-2,ky-3; kx-2,ky+2 PLOT kx+2,ky-3; kx+2,ky+2 PLOT kx-1,ky-4; kx+1,ky-4 SET COLOR white PLOT kx,ky-4; kx,ky-6 PLOT kx,ky+3; kx,ky-2 END IF END SUB ! SUB hLine(name$,axy,wLft,wRgt,fLft,fRgt,xf,x2,x4,xl,pclr,slot$) ! CALL ClearSlider(wLft,wRgt,axy+15,axy-15) ! CALL SliderSlotH(wLft-3,wRgt+3,axy) ! CALL SliderAxisHrz(wLft,wRgt,axy,fLft,fRgt,xf,x2,x4,xl) ! CALL SetTextFont(1,12,"bold") ! CALL PlotTextRJ(wLft-5,axy+9,name$,pclr) ! END SUB ! ! SUB hLinePi(hname$,axy,wLft,wRgt,fLft,fRgt,xf,x2,x4,xl,pclr,slot$) ! CALL ClearSlider(wLft,wRgt,axy+15,axy-15) ! CALL SliderSlotH(wLft-3,wRgt+3,axy) ! CALL SliderAxisHrzPi(wLft,wRgt,axy,fLft,fRgt,xf,x2,x4,xl) ! CALL SetTextFont(1,12,"bold") ! CALL PlotTextRJ(wLft-5,axy+9,name$,pclr) ! END SUB ! SUB hMark(n,m,wLft,wRgt,wBas,fLft,ratio,form$,slot$,clr) ! LET wx= wLft + ratio*(n-fLft) ! BOX SHOW slot$ at wLft-3,wBas ! CALL SliderKnob(wx,wBas-5) ! BOX CLEAR wRgt+3,wRgt+50,wBas,wBas-12 ! CALL SetTextFont(1,12,"bold") ! CALL PlotTextLJ(wRgt+7,wBas-4,trim$(using$(form$,n*m)),clr) ! END SUB ! SUB ClearSlider(wLft,wRgt,wBas,wTop) ! BOX CLEAR wLft-20,wRgt+35,wBas,wTop ! END SUB ! PICTURE KnobH(clr) ! SET COLOR drkgry ! PLOT -2,-2; -2,3 ! PLOT -1,-3; -1,4 ! PLOT 0,-5; 0,4 ! PLOT 1,-3; 1,4 ! PLOT 2,-2; 2,3 ! SET COLOR white ! PLOT 0,-5; 0,2 ! END PICTURE ! PICTURE PhaseMark(clr) ! SET COLOR clr ! PLOT -2,-2; -2,2 ! PLOT -1,-3; -1,3 ! PLOT 0,-5; 0,5 ! PLOT 1,-3; 1,3 ! PLOT 2,-2; 2,2 ! END PICTURE END MODULE MODULE TypeSetting DECLARE PUBLIC black,drkgry,drkmid,midgry,litmid,litgry,white DECLARE PUBLIC red,yellow,green,cyan,blue,magenta,colorscheme DECLARE PUBLIC planeclr,gridclr,rimclr,axisclr,axislabelclr,titleclr ! ---- type subs ---- SUB PlotTextLJ(tx,ty,t$,clr) DECLARE PUBLIC M68kFlag SET COLOR clr IF M68kflag=1 then PLOT TEXT, AT tx,ty: t$ ! 4.04 ELSE PLOT TEXT, AT tx,ty+1: t$ ! TB5 END IF END SUB SUB PlotTextCJ(midb,ty,t$,clr) CALL StringWidth(t$,sl) LET Lft= midb - int(sl/2) CALL PlotTextLJ(Lft,ty,t$,clr) END SUB SUB PlotTextRJ(Rgt,ty,t$,clr) CALL StringWidth(t$,sl) LET Lft= Rgt - sl CALL PlotTextLJ(Lft,ty,t$,clr) END SUB SUB DropTextLJ(tx,ty,t$,clr) CALL PlotTextLJ(tx+1,ty+1,t$,black) CALL PlotTextLJ(tx,ty,t$,clr) END SUB SUB DropTextCJ(midb,ty,t$,clr) CALL StringWidth(t$,sl) LET Lft= midb - int(sl/2) CALL DropTextLJ(Lft,ty,t$,clr) END SUB SUB DropTextRJ(Rgt,ty,t$,clr) CALL StringWidth(t$,sl) LET Lft= Rgt - sl CALL DropTextLJ(Lft,ty,t$,clr) END SUB SUB AlignEqual(ex,ey,eq$,ec) LET e = pos(eq$,"=") LET Lft$= eq$(1:e) CALL StringWidth(Lft$,sl) CALL PlotTextLJ(ex-sl,ey,eq$,ec) END SUB SUB AlignDot(dx,dy,n$,dc) LET d = pos(n$,".") LET Lft$= n$(1:d) CALL StringWidth(Lft$,sl) CALL PlotTextLJ(dx-sl,dy,n$,dc) END SUB SUB MeasureSuperSub(in$,pix) LET out$= in$ LET l = len(out$) LET extra$= "" FOR i= 1 to l LET chr$= out$(i:i) IF chr$="^" or chr$="_" or chr$="[" or chr$="]" then LET extra$= extra$ & chr$ !LET out$(i:i)= "" END IF NEXT i CALL SetTextFont(1,12,"bold") CALL StringWidth(in$,sl) CALL StringWidth(extra$,el) LET pix= sl-el END SUB SUB SuperSubScriptLJ(sx,sy,super$,sclr) CALL SuperSubScript(sx,sy,super$,sclr) END SUB SUB SuperSubScriptRJ(sx,sy,super$,sclr) CALL MeasureSuperSub(super$,sl) LET Lft= sx - sl CALL SuperSubScript(Lft,sy,super$,sclr) END SUB SUB SuperSubScriptCJ(sx,sy,super$,sclr) CALL MeasureSuperSub(super$,sl) LET Lft= sx - int(sl/2) CALL SuperSubScript(Lft,sy,super$,sclr) END SUB SUB SuperSubScript(sx,sy,super$,sclr) LOCAL p,last,p0,p1,p2,ll,el,rl LOCAL Lft$,exp$,Rgt$ LET Lft$,exp$,Rgt$= "" LET p0 = pos(super$,"^") LET p1 = pos(super$,"_") LET last= len(super$) CALL SetTextFont(1,12,"bold") IF p0=0 and p1=0 then CALL PlotTextLJ(sx,sy,super$,sclr) ELSE IF p0=0 then LET p= p1 ELSE IF p1=0 then LET p= p0 ELSE LET p= min(p1,p0) END IF LET Lft$= super$(1:p-1) LET p2 = pos(super$,"]") IF p2=0 then EXIT SUB LET exp$= super$(p+2:p2-1) LET Rgt$= super$(p2+1:last) CALL StringWidth(Lft$,ll) CALL PlotTextLJ(sx,sy,Lft$,sclr) CALL SetTextFont(1,9,"bold") CALL StringWidth(exp$,el) IF p0=p then CALL PlotTextLJ(sx+ll,sy-5,exp$,sclr) ELSE IF p1=p then CALL PlotTextLJ(sx+ll,sy+3,exp$,sclr) END IF CALL SetTextFont(1,12,"bold") IF Rgt$<>"" then CALL SuperSubScriptLJ(sx+ll+el,sy,Rgt$,sclr) END IF END IF END SUB SUB CmplxStr(cmplx$,real,imag) ! Numeric to text LET x$= using$("-%.##",round(real,2)) LET y$= using$("%.##",abs(round(imag,2))) IF imag<0 then LET sign$= " - " else LET sign$= " + " LET cmplx$= "( " & x$ & sign$ & y$ & " i )" END SUB SUB CmplxText(Lft,Bas,rl,im,clr) LET rldot = Lft+35 LET imdot = rldot+50 LET ix = imdot+20 LET txtrl$= trim$(using$("--%.##",rl)) LET txtim$= trim$(using$("++%.##",im)) CALL PlotTextLJ(Lft,Bas,"(",clr) CALL AlignDot(rldot,Bas,txtrl$,clr) CALL AlignDot(imdot,Bas,txtim$,clr) CALL PlotTextLJ(ix,Bas,"i )",clr) END SUB SUB VarDot(lft,bas,v$,clr) SET COLOR clr CALL PlotTextLJ(lft,bas,v$,clr) CALL StringWidth(".",sl) IF v$="t" then CALL PlotTextLJ(lft+3-sl/2,bas-10,".",clr) ELSE CALL PlotTextLJ(lft+4-sl/2,bas-9, ".",clr) END IF END SUB SUB AlignSuper(asx,ty,pnts,clr,txt$) LET chx= asx SET COLOR clr LET lft$,rgt$,exp$= "" LET l = len(txt$) LET p = pos(txt$,"^") IF p=0 then SET COLOR clr !PLOT TEXT, AT chx,ty: txt$ CALL PlotTextLJ(chx,ty,txt$,clr) ELSE LET lft$= txt$(1:p-1) LET rgt$= txt$(p+1:l) LET l = len(rgt$) IF rgt$(1:1)= "(" then LET p = pos(rgt$,")") LET exp$= rgt$(1:p) LET rgt$= rgt$(p+1:l) ELSE IF len(rgt$)=1 then LET exp$= rgt$ LET rgt$="" ELSE FOR i= 2 to l LET c$= rgt$(i:i) IF c$= " " or c$= "+" or c$= "-" or c$= "*" or c$= ")" then LET exp$= rgt$(1:i-1) LET rgt$= rgt$(i:l) EXIT FOR ELSE IF i=l then LET exp$= rgt$(1:l) LET rgt$= "" EXIT FOR END IF NEXT i END IF CALL SetTextFont(1,12,"bold") CALL StringWidth(lft$,ll) CALL PlotTextLJ(chx,ty,lft$,clr) CALL SetTextFont(1, 9,"bold") CALL StringWidth(exp$,el) CALL PlotTextLJ(chx+ll,ty-pnts,exp$,clr) LET chx= chx+ll+el LET txt$= rgt$ LET lft$,rgt$,exp$= "" LET l = len(txt$) LET p = pos(txt$,"^") IF p<>0 then LET lft$= txt$(1:p-1) LET rgt$= txt$(p+1:l) LET l = len(rgt$) IF rgt$(1:1)= "(" then LET p = pos(rgt$,")") LET exp$= rgt$(1:p) LET rgt$= rgt$(p+1:l) ELSE IF len(rgt$)=1 then LET exp$= rgt$ LET rgt$="" ELSE FOR i= 2 to l LET c$= rgt$(i:i) IF c$= " " or c$= "+" or c$= "-" or c$= "*" or c$= ")" then LET exp$= rgt$(1:i-1) LET rgt$= rgt$(i:l) EXIT FOR ELSE IF i=l then LET exp$= rgt$(1:l) LET rgt$= "" EXIT FOR END IF NEXT i END IF CALL SetTextFont(1,12,"bold") CALL StringWidth(lft$,ll) CALL PlotTextLJ(chx,ty,lft$,clr) CALL SetTextFont(1, 9,"bold") CALL StringWidth(exp$,el) CALL PlotTextLJ(chx+ll,ty-pnts,exp$,clr) CALL SetTextFont(1,12,"bold") !PLOT TEXT, AT chx+ll+el,ty: rgt$ CALL PlotTextLJ(chx+ll+el,ty, rgt$,clr) ELSE CALL SetTextFont(1,12,"bold") !PLOT TEXT, AT chx,ty: txt$ CALL PlotTextLJ(chx,ty, txt$,clr) END IF END IF END SUB END MODULE MODULE PopMenu PUBLIC mlft,mrgt,mbas,mtop PUBLIC poplft,poprgt,popbas,poptop,pophgt,PopUp$ PUBLIC listlft,listrgt,listbas,listtop,eqspc,eqlft DECLARE PUBLIC workLft,workRgt,workBas,workTop,workMidx ! work area DECLARE PUBLIC black,drkgry,drkmid,midgry,litmid,litgry,white DECLARE PUBLIC red,yellow,green,cyan,blue,magenta,colorscheme DECLARE PUBLIC planeclr,gridclr,rimclr,axisclr,axislabelclr,titleclr,rightsclr DECLARE PUBLIC numberlineclr,slotdrkclr,slotlgtclr,slideclr,headerclr,btnclr SUB InitMenu(mlft,mtop,eq,eqcnt,names$(,)) LET mrgt= mlft+16 LET mbas= mtop+16 CALL ShowEquation(eq,names$) CALL PopUpTab(mlft,mrgt,mbas,mtop) CALL DrawPopUp(names$(,),eqcnt) END SUB SUB PopUpTab(mlft,mrgt,mbas,mtop) CALL DrawButton(mlft,mrgt,mbas,mtop,0,"") LET sx = mlft+8 LET sy = mtop+9 LET syu= sy + 2 SET COLOR white FOR i= 0 to 5 PLOT sx-i,syu-i; sx+i,syu-i NEXT i END SUB SUB DrawPopUp(names$(,),eqcnt) LOCAL temp$,oldenum,mx,my,ms,et,eb,clearflag,ty LET maxsw= 0 FOR i= 1 to eqcnt LET eq$= names$(i,0) CALL StringWidth(eq$,sw) IF sw>maxsw then LET maxsw= sw NEXT i CALL SetTextFont(1,12,"bold") LET eqspc = 21 LET pophgt= eqcnt*eqspc + 26 IF pophgt>workbas-mtop then LET popbas= workbas ELSE LET popbas= mtop+pophgt END IF LET poptop = popbas-pophgt LET poplft = mrgt+2 LET poprgt = poplft+maxsw+100 LET listlft= poplft+ 5 LET listrgt= poprgt- 5 LET listtop= poptop+22 LET listbas= listtop + eqcnt*eqspc LET eqlft = poplft+80 BOX KEEP poplft,poprgt,popbas,poptop in temp$ CALL Rim(poplft,poprgt,popbas,poptop,0) CALL DrawButton(poplft+5,poplft+20,poptop+20,poptop+5,0,"") FOR i= 1 to eqcnt LET et= listtop + (i-1)*eqspc LET eb= et+eqspc-1 CALL DrawButton(poplft+5,poprgt-5,eb,et,0,"") CALL DiffEq(eqlft,eb-6,white,"y",i,names$) NEXT i BOX KEEP poplft,poprgt,popbas,poptop in PopUp$ CALL Delay(1) BOX SHOW temp$ at poplft,popbas LET temp$= "" END SUB SUB ShowEquation(eq,names$(,)) CALL SetTextFont(1,12,"bold") BOX CLEAR mrgt+2,mrgt+150,mbas,mtop CALL DiffEq(mrgt+75,mbas-4,white,"y",eq,names$) END SUB SUB DiffEq(dex,dey,clr,x$,eq,name$(,)) CALL SetTextFont(1,12,"bold") IF eq<5 then CALL PlotTextRJ(dex,dey,"dy/dt =",clr) ELSE CALL PlotTextRJ(dex,dey,"t dy/dt =",clr) END IF LET eq$= name$(eq,0) ! CALL VarDot(dex,dey,x$,clr) !LET e$= x$&" = " & eq$ !let lft$= x$ & " =" ! CALL PlotTextRJ(dex,dey,lft$,clr) ! LET p = pos(e$,"^") ! CALL SetTextFont(1,12,"bold") ! IF p<>0 then ! LET e$= " = " & eq$ ! CALL AlignSuper(dex+8,dey,4,clr,eq$) ! ELSE CALL SuperSubScriptLJ(dex+8,dey,eq$,clr) ! END IF END SUB SUB ReadPanel(eq) LOCAL temp$,oldenum,mx,my,ms,et,eb,clearflag,ty BOX KEEP poplft,poprgt,popbas,poptop in temp$ BOX SHOW PopUp$ at poplft,popbas LET oldenum= -99999 CALL MouseDown(mx,my,ms) DO GET MOUSE: mx,my,ms IF mx>listlft and mxlisttop and myoldenum then LET oldenum= enum LET et = listtop + eqspc*(enum-1) LET eb = et+eqspc-1 BOX SHOW PopUp$ at poplft,popbas CALL ButtonDown(listlft,listrgt,eb,et) LET clearflag= 1 END IF ELSE IF clearflag=1 then LET clearflag= 0 LET oldenum = -99999 BOX SHOW PopUp$ at poplft,popbas END IF END IF LOOP until ms=3 IF mx>listlft and mxlisttop and my"" then LET midb= int((l+r)/2) + 1 CALL PlotTextCJ(midb,b-tBas,t$,white) END IF END SUB ! ---- Mouse Button Action ---- SUB ButtonDown(l,r,b,t) SET COLOR black PLOT l,b; l,t; r,t SET COLOR white PLOT l+1,b; r,b; r,t+1 END SUB SUB ButtonUp(l,r,b,t) SET COLOR white PLOT l,b; l,t; r,t SET COLOR 1 PLOT l+1,b; r,b; r,t+1 END SUB SUB MouseButtonUp(Lft,Rgt,Bas,Top,ms) CALL ButtonDown(Lft,Rgt,Bas,Top) DO GET MOUSE: mx,my,ms LOOP until ms=3 CALL ButtonUp(Lft,Rgt,Bas,Top) END SUB SUB MouseUp(mx,my,ms) DO GET MOUSE: mx,my,ms LOOP until ms=3 END SUB SUB MouseDown(mx,my,ms) DO GET MOUSE: mx,my,ms LOOP until ms=2 END SUB END MODULE MODULE Clipping SUB Clip(w1fLft,w1fRgt,w1fBas,w1fTop,x1,y1,x2,y2) ! (x1,y1) is within LET yd= y2-y1 LET xd= x2-x1 IF xd<>0 and yd<>0 then LET ym= yd/xd LET xm= xd/yd IF x2w1fRgt then ! clip right LET y2= y1 + ym*(w1fRgt-x1) LET x2= w1fRgt END IF IF y2w1fTop then ! clip top LET x2= x1 + xm*(w1fTop-y1) LET y2= w1fTop END IF ELSE IF yd=0 then IF x2w1fRgt then ! clip right LET y2= y1 + ym*(w1fRgt-x1) LET x2= w1fRgt END IF END IF END SUB END MODULE MODULE Timer SHARE unitTicks SUB SetTimer LET ticks1= 0 LET t1 = time DO LET ticks1= ticks1+1 LOOP until time-t1>=1 LET ticks2= 0 LET t1 = time DO LET ticks2= ticks2+1 LOOP until time-t1>=1 LET unitTicks= (ticks1+ticks2)/2 END SUB SUB Delay(ticksn) LOCAL count,tick,killTime LET count= round(ticksn*unitTicks) FOR tick= 1 to count LET killTime= time-t1 LET killTime= time-t1 NEXT tick END SUB END MODULE ! *** SUB TargetPractice DECLARE PUBLIC PCflag,M68kflag,Mac5flag,xmax,ymax DECLARE PUBLIC toolLft,toolRgt,toolBas,toolTop,toolhdr,toolHgt,toolWid ! tool boundaries DECLARE PUBLIC winLft,winRgt,winBas,winTop,winHgt,winWid ! window DECLARE PUBLIC workLft,workRgt,workBas,workTop,workMidx ! work area DECLARE PUBLIC qLft,qRgt,qBas,qTop DECLARE PUBLIC infLft,infRgt,infBas,infTop DECLARE PUBLIC iLft,iRgt,iBas,iTop DECLARE PUBLIC black,drkgry,drkmid,midgry,litmid,litgry,white DECLARE PUBLIC red,yellow,green,cyan,blue,magenta,colorscheme DECLARE PUBLIC planeclr,gridclr,rimclr,axisclr,axislabelclr,titleclr,rightsclr DECLARE PUBLIC numberlineclr,slotdrkclr,slotlgtclr,slideclr,headerclr,btnclr DECLARE PUBLIC mlft,mrgt,mbas,mtop DECLARE PUBLIC poplft,poprgt,popbas,poptop,pophgt,PopUp$ DECLARE PUBLIC listlft,listrgt,listbas,listtop,eqspc,eqlft ! --- info page --- DIM Info$(1:9) MAT READ Info$ DATA "Target Practice" DATA "" DATA "Click the pop-up menu button [v] to open the list of equations." DATA "Click an equation to select it." DATA "" DATA "Click the (t,y) plane to set a new target or draw solutions." DATA "Click [Clear Curves] to remove solutions from the plane." DATA "Click [Clear Target] to remove the current target from the plane." DATA "Click [Draw Field] to toggle the field on or off." ! --- utility functions --- DEF clamp(n,lo,hi)= min(max(n,lo),hi) DEF SetColor(n) = red + n ! --- Functions --- ! Equation list panel LET eqcnt= 4 LET eq = 1 DEF dxdt1(t,x) = 1 ! unstable node DEF dydt1(t,x) = x*x-t DEF dxdt2(t,x) = 1 ! stable node DEF dydt2(t,x) = sin(x) DEF dxdt3(t,x) = 1 ! saddle DEF dydt3(t,x) IF cos(t)<>0 then LET dydt3= x/cos(t) ELSE LET dydt3= 99999 END IF END DEF DEF dxdt4(t,x) = 1 ! unstable star DEF dydt4(t,x) = x*(1-x) DEF dxdt5(t,x) = 1 ! stable star DEF dydt5(t,x) = 1+x*x ! DEF dxdt6(t,x) = 1 ! unstable improper node ! DEF dydt6(t,x) ! IF x*x-t*t>=0 then ! LET dydt6= sqr(x*x-t*t) ! ELSE ! LET dydt6= 99999 ! END IF ! END DEF DEF dxdt6(t,x) = 1 ! stable improper node DEF dydt6(t,x) = tan(t) DEF dxdt7(t,x) = 1 ! unstable focus DEF dydt7(t,x) IF t*x>=0 then LET dydt7= sqr(t*x) ELSE LET dydt7= 99999 END IF END DEF DEF dxdt8(t,x) = 1 ! stable focus DEF dydt8(t,x) IF x<>0 then IF t/x>=0 then LET dydt8= sqr(t/x) ELSE LET dydt8= 99999 END IF ELSE LET dydt8= 99999 END IF END DEF DEF dxdt9(x,y)= 1 ! unstable node DEF dydt9(x,y) IF x<>0 then LET dydt9= y/x END DEF DEF dxdt10(x,y)= 1 ! unstable node DEF dydt10(x,y) IF x<>0 then LET dydt10= y/x + x*cos(x) END DEF ! SUB FetchEq(eq,x,y,dx,dy) ! SELECT CASE eq ! CASE 1 ! LET dx= dxdt1(x,y) ! LET dy= dydt1(x,y) ! CASE 2 ! LET dx= dxdt2(x,y) ! LET dy= dydt2(x,y) ! END SELECT ! END SUB SUB FetchEq(eq,t,x,dx,dy) SELECT CASE eq CASE 1 LET dx= dxdt1(t,x) LET dy= dydt1(t,x) CASE 2 LET dx= dxdt2(t,x) LET dy= dydt2(t,x) ! CASE 3 ! LET dx= dxdt3(t,x) ! LET dy= dydt3(t,x) CASE 3 LET dx= dxdt4(t,x) LET dy= dydt4(t,x) CASE 4 LET dx= dxdt5(t,x) LET dy= dydt5(t,x) ! CASE 6 ! LET dx= dxdt6(t,x) ! LET dy= dydt6(t,x) ! CASE 7 ! LET dx= dxdt7(t,x) ! LET dy= dydt7(t,x) ! CASE 8 ! LET dx= dxdt8(t,x) ! LET dy= dydt8(t,x) CASE 5 LET dx= dxdt9(t,x) LET dy= dydt9(t,x) CASE 6 LET dx= dxdt10(t,x) LET dy= dydt10(t,x) ! CASE 11 ! LET dx= dxdt11(t,x) ! LET dy= dydt11(t,x) END SELECT IF dy=99999 then LET err= 1 else LET err= 0 END SUB ! Graphing window LET wsize = 360 LET pisize= 176 ! --- plane 1 data --- DECLARE PUBLIC w1Lft, w1Rgt, w1Bas, w1Top, w1Midx, w1Midy DECLARE PUBLIC w1fLft, w1fRgt, w1fBas, w1fTop, w1x0, w1y0 DECLARE PUBLIC w1Firstx, w1SxTik, w1LxTik, w1xLabel, w1xGridstep DECLARE PUBLIC w1Firsty, w1SyTik, w1LyTik, w1yLabel, w1yGridstep DECLARE PUBLIC w1xPiFlag, w1wHgt, w1fHgt LET w1Flag = 1 LET w1xPiFlag= 0 LET w1Lft= workLft + 60 ! pixel bounds LET w1Rgt= w1Lft + wsize LET w1Top= workTop + 35 LET w1Bas= w1Top + wsize LET sz= 3 LET w1fLft= -sz ! function bounds LET w1fRgt= sz LET w1fTop= sz LET w1fBas= -sz LET w1Xax$= "t" ! axis labels LET w1Yax$= "y" LET w1xGridstep= 0 ! horizontal grid intervals LET w1yGridstep= 0 ! vertical grid intervals LET w1SxTik = 0.5 ! horizontal axis Tik marks LET w1LxTik = 1 LET w1xLabel= 1 LET w1Firstx= w1fLft LET w1SyTik = 0.5 ! vertical axis Tik marks LET w1LyTik = 1 LET w1yLabel= 1 LET w1Firsty= w1fBas ! --- Plane 1 methods --- DECLARE DEF w1fncx,w1fncy,w1wndx,w1wndy ! window/function transforms CALL Plane1Variables SUB Plane1Clear BOX CLEAR w1Lft-20,w1Rgt+20,workBas-5,w1Top-20 END SUB SUB Plane1Init CALL DrawPlane1(0,1,1) ! grid, axes, zeroaxes CALL SetTextFont(1,12,"bold") CALL PlotTextLJ(w1Rgt+8,w1y0+3,w1Xax$,axislabelclr) ! axis labels CALL PlotTextCJ(w1x0,w1Top-10,w1Yax$,axislabelclr) BOX KEEP w1Lft-5,w1Rgt+5,w1Bas+5,w1Top-5 in w1gridLayer$ LET w1graphLayer$= w1gridLayer$ END SUB SUB Plane1Refresh BOX SHOW w1gridLayer$ at w1Lft-5,w1Bas+5 END SUB SUB Plane1Graph BOX SHOW w1graphLayer$ at w1Lft-5,w1Bas+5 END SUB ! ------------------------------------------ ! --- plane 2 data --- DECLARE PUBLIC w2Lft,w2Rgt,w2Bas,w2Top,w2Midx,w2Midy DECLARE PUBLIC w2fLft,w2fRgt,w2fBas,w2fTop,w2fHgt,w2fWid,w2x0,w2y0 DECLARE PUBLIC w2Firstx, w2SxTik, w2LxTik, w2xLabel, w2xGridstep DECLARE PUBLIC w2Firsty, w2SyTik, w2LyTik, w2yLabel, w2yGridstep DECLARE PUBLIC w2xPiFlag LET w2Flag= 1 LET w2Lft = w1Rgt + 70 ! pixel bounds LET w2Rgt = w2Lft + 200 LET w2Top = w1Midy -100 LET w2Bas = w2Top + 200 LET sz2 = 4 LET w2fLft= -sz2 ! function bounds * pi LET w2fRgt= sz2 LET w2fTop= sz2 LET w2fBas= -sz2 LET w2xGridstep= 0 ! horizontal grid intervals LET w2yGridstep= 0 ! vertical grid intervals LET w2xAx$= "t" ! axis labels LET w2yAx$= "y" LET w2SxTik = .5 ! horizontal axis Tik marks LET w2LxTik = 1 LET w2xLabel= 1 LET w2Firstx= w2fLft LET w2xPiFlag= 0 LET w2SyTik = .5 ! vertical axis Tik marks LET w2LyTik = 1 LET w2yLabel= 1 LET w2Firsty= w2fBas ! --- plane 2 methods --- DECLARE DEF w2fncx, w2fncy, w2wndx, w2wndy ! window/function transforms CALL Plane2Variables SUB Plane2Init CALL DrawPlane2(0,0,0) ! grid, axes, zeroaxes ! CALL Grid2 ! SET COLOR black ! BOX AREA w2Lft,w2Rgt,w2Bas,w2Top CALL SetTextFont(1,12,"bold") ! CALL PlotTextLJ(w2Rgt+8,w2y0+3,w2xAx$,axislabelclr) ! axis labels ! CALL PlotTextCJ(w2x0,w2Top-10,w2yAx$,axislabelclr) BOX KEEP w2Lft-5,w2Rgt+5,w2Bas+5,w2Top-5 in w2gridLayer$ END SUB SUB Plane2Refresh BOX SHOW w2gridLayer$ at w2Lft-5,w2Bas+5 END SUB SUB Plane2Graph BOX SHOW w2graphLayer$ at w2Lft-5,w2Bas+5 END SUB SUB Grid2 LOCAL i,wx SET COLOR drkmid FOR i= 1 to 2 step 1 LET wx= w2wndx(i) PLOT wx,w2Top+1; wx,w2Bas-1 NEXT i END SUB ! ---------- Text Output Rects ---------- ! --- text rectangle 1 --- LET t1BasLn = w1Midy + 3 LET t1Lft = w1Rgt + 60 LET t1Rgt = workrgt LET t1Bas = t1BasLn + 5 LET t1Top = t1BasLn - 15 LET t1Label$= "" SUB t1Label CALL SuperSubScriptRJ(t1Lft,t1BasLn,t1Label$,white) END SUB SUB t1Set(t1x,t1y) CALL SetTextFont(1,12,"bold") LOCAL r$,i$,z$ LET x$= trim$(using$(f$,t1x)) LET y$= trim$(using$(f$,t1y)) LET txt$= "target: (" & x$ & ", " & y$ & ")" CALL t1Clear CALL PlotTextLJ(t1Lft,t1BasLn,txt$,cyan) END SUB SUB t1Clear BOX CLEAR t1Lft-2,t1Rgt,t1Bas,t1Top END SUB SUB t1Init CALL t1Set(0,0) END SUB ! ! --- text rectangle 2 --- ! ! LET t2BasLn= t1BasLn + 30 ! LET t2Lft = t1Lft ! LET t2Rgt = t1Rgt ! LET t2Bas = t2BasLn + 5 ! LET t2Top = t2BasLn - 15 ! LET t2Label$ = "" ! ! SUB t2Label ! CALL SuperSubScriptRJ(t2Lft,t2BasLn,t2Label$,cclr) ! END SUB ! ! SUB t2Set ! CALL SetTextFont(1,14,"bold") ! LOCAL r$,i$,z$ ! CALL t2Clear ! ! LET t$ = "dy/dt = v" ! CALL PlotTextLJ(t2Lft,t2BasLn,t$,white) ! END SUB ! ! SUB t2Clear ! BOX CLEAR t2Lft-2,t2Rgt,t2Bas,t2Top ! END SUB ! ! SUB t2Init ! CALL t2Label ! CALL t2Set ! END SUB ! --- Dynamic output --- LET lnspc= 14 LET ivlft= w1Rgt + 60 LET ivrgt= ivlft+80 LET ivtop= w1Midy+15 LET ivbas= ivtop + 6*lnspc LET iveqx= ivlft+27 LET ivdot= ivlft+55 LET ivf$ = "--%.##" LET ivtbas1= ivtop + 1*lnspc LET ivtbas2= ivtop + 2*lnspc LET ivtbas4= ivtop + 4*lnspc LET ivtbas5= ivtop + 5*lnspc ! ! SUB ClearInit ! !BOX CLEAR ivlft-15,ivrgt+15,ivbas,ivtop-20 ! ! LET initx,inity = 0 ! ! LET initmx,initmy= 0 ! END SUB ! ! SUB ShowInit(x0,y0) ! CALL ClearInit ! SET COLOR white ! PLOT ivlft,ivtop-2; ivrgt,ivtop-2 ! CALL SetTextFont(1,9,"bold") ! CALL PlotTextLJ(ivlft,ivtop-7,"Initial Values",white) ! ! CALL SetTextFont(1,12,"bold") ! CALL PlotTextRJ(iveqx,ivtbas1,"t = ",white) ! CALL PlotTextRJ(iveqx,ivtbas2,"y = ",white) ! CALL AlignDot(ivdot,ivtbas1,using$(ivf$,x0),white) ! CALL AlignDot(ivdot,ivtbas2,using$(ivf$,y0),white) ! ! CALL FetchEq(eq,x0,y0,dx,dy) ! ! ! CALL VarDot(ivlft,ivtbas4,"t",white) ! ! CALL PlotTextRJ(iveqx,ivtbas4,"= ",white) ! ! CALL AlignDot(ivdot,ivtbas4,using$(ivf$,dx),white) ! ! ! ! CALL VarDot(ivlft,ivtbas5,"y",white) ! ! CALL PlotTextRJ(iveqx,ivtbas5,"= ",white) ! ! CALL AlignDot(ivdot,ivtbas5,using$(ivf$,dy),white) ! END SUB ! --- LET tlft = ivLft LET trgt = tlft+80 LET ttop = w1Top-5 LET tbas = ttop + 6*lnspc LET teqx = tlft+27 LET tdot = tlft+60 LET tbas1= ttop + 1*lnspc LET tbas2= ttop + 2*lnspc LET tbas4= ttop + 4*lnspc LET tbas5= ttop + 5*lnspc SUB ClearOutput BOX CLEAR teqx,trgt,tbas,ttop+5 END SUB SUB OutputLabels CALL SetTextFont(1,12,"bold") CALL PlotTextRJ(teqx,tbas1,"t = ",white) CALL PlotTextRJ(teqx,tbas2,"y = ",white) !CALL VarDot(tlft,tbas4,"t",white) !CALL PlotTextRJ(teqx,tbas4,"t' = ",white) !CALL VarDot(tlft,tbas5,"y",white) !CALL PlotTextRJ(teqx,tbas5,"y' = ",white) !CALL CursorValues(x1,y1,dx,dy,dt) END SUB SUB CursorValues(x1,y1,dx,dy,dt) CALL ClearOutput CALL SetTextFont(1,12,"bold") LET f$= "---%.##" CALL AlignDot(tdot,tbas1,using$(f$,x1),white) CALL AlignDot(tdot,tbas2,using$(f$,y1),white) ! IF dx<1000 then ! CALL AlignDot(tdot,tbas4,using$(f$,dx),white) ! END IF ! IF abs(dy)<1000 then ! CALL AlignDot(tdot,tbas5,using$(f$,dy),white) ! END IF END SUB SUB ClearPrompt BOX CLEAR plft,prgt,pbas,ptop END SUB SUB Prompt CALL SetTextFont(1,12,"bold") LET plft = w1Lft LET prgt = w1Rgt LET pbasln= workbas-10 LET ptop = pbasln -15 LET pbas = pbasln + 5 CALL ClearPrompt IF targstate=0 and eq<8 then CALL PlotTextLJ(plft,pbasln,"Click the plane to set a target point...",white) ! ELSE ! CALL PlotTextLJ(plft,pbasln,"Click another point on the same trajectory...",white) END IF END SUB ! --- popup menu --- LET mlft = w1Rgt + 60 LET mtop = w1Midy + 30 LET mrgt = mlft + 16 LET mbas = mtop + 16 DIM eqnames$(1:4,0:1) MAT READ eqnames$ DATA "y^[2]-t","" DATA "sin(y)","" ! DATA "y/cos(t)","" DATA "y(1-y)","" DATA "1+y^[2]","" ! DATA "sqrt(y^[2]-t^[2])","" ! DATA "tan(t)","" ! DATA "sqrt(ty)","" ! DATA "sqrt(t/y)","" ! DATA "y","" ! DATA "y + t^[2] cos(t)","" SUB ClearEquation BOX CLEAR mrgt+8,workrgt,mbas,mtop END SUB ! --- init --- LET fieldclr = drkmid LET fieldState= 1 LET menuState = 0 LET dt,tstp = 1/16 LET targstate = 0 LET f$ = "-%.##" LET elft = tlft+32 CALL SetTimer CALL InitScreen CALL Prompt CALL InitMenu(mlft,mtop,eq,eqcnt,eqnames$) SUB InitScreen CALL Plane1Init !CALL t1Init CALL Buttons CALL OutputLabels CALL VectorField LET mrgt= mlft+16 LET mbas= mtop+16 CALL ShowEquation(eq,eqnames$) CALL PopUpTab(mlft,mrgt,mbas,mtop) IF targstate=1 then DRAW target(cyan) with shift(wtx1,wty1) CALL t1Set(ftx1,fty1) BOX KEEP w1Lft-5,w1Rgt+5,w1Bas+5,w1Top-5 in w1GraphLayer$ END IF END SUB ! --- event loop --- DO CALL SetTextFont(1,12,"bold") LET clearflag= 0 IF ms<>2 then DO GET MOUSE: mx,my,ms IF mx>w1Lft-3 and mxw1Top-3 and myoldx or my<>oldy then LET oldx= mx LET oldy= my LET x = w1Fncx(mx) LET y = w1Fncy(my) IF targstate=0 then CALL Plane1Graph !CALL Plane2Refresh IF abs(mx-w1x0)<3 and abs(my-w1y0)<3 then LET x= 0 LET y= 0 LET mx= w1x0 LET my= w1y0 END IF DRAW target(cyan) with shift( mx,my ) CALL t1Set(x,y) ELSE LET x= round(x,5) LET y= round(y,5) IF round(y,2)<>0 then CALL FetchEq(eq,x,y,dx,dy) LET ang= angle(dx,-dy) LET sa = 5*sin(ang) LET wy1= my+sa LET wy2= my-sa LET ca = 5*cos(ang) ! CALL VectorClip(t,x,dt,6) CALL CursorValues(x,y,dx,dy,dt) CALL Plane1Graph SET COLOR yellow PLOT mx-ca,wy2; mx+ca,wy1 END IF END IF LET clearflag= 1 END IF ELSE IF clearflag=1 then CALL CleanUp END IF LOOP until ms=2 IF clearflag=1 then CALL CleanUp END IF IF mx>w1Lft-3 and mxw1Top-3 and myw1Lft-3 and mxw1Top-3 and mycclft and mxcctop and myctlft and mxcttop and mysllft and mxsltop and mygLft and mxgTop and mymlft and mxmtop and my8 then ! CALL ClearPrompt ! CALL SetTarget(0,0,w1x0,w1y0) ! ELSE CALL Prompt ! END IF ELSE IF mx>inflft and mxinftop then CALL MouseButtonUp(inflft,infrgt,infbas,inftop,ms) CALL InfoPage(Info$()) CALL InitScreen ELSE IF mx>qlft and mxqtop then CALL MouseButtonUp(qlft,qrgt,qbas,qtop,ms) EXIT SUB ELSE CALL MouseUp(mx,my,ms) END IF LOOP SUB SetTarget(ftx,fty,wx,wy) LET targstate= 1 LET ftx1= w1Fncx(wx) LET fty1= w1Fncy(wy) LET wtx1= wx LET wty1= wy CALL t1Set(ftx1,fty1) DRAW target(cyan) with shift(wtx1,wty1) CALL ClearAndReset BOX KEEP w1Lft-5,w1Rgt+5,w1Bas+5,w1Top-5 in w1GraphLayer$ END SUB SUB CleanUp CALL Plane1Graph CALL ClearOutput IF targstate=0 then CALL t1Clear END IF LET clearflag= 0 END SUB ! --- targets --- SUB Trajectory(x0,y0,clr) CALL CursorValues(x,y,dx,dy,dt) LET hitflag = 0 LET exitflag= 0 LET near = 2 ! was .6 LET e = .005 LET adt = abs(dt) LET tclr= yellow LET dt = adt CALL Flow(x0,y0) LET dt = -adt ! now negative time CALL Flow(x0,y0) IF hitflag=1 then LET tclr= green LET dt = adt CALL Flow(x0,y0) LET dt = -adt CALL Flow(x0,y0) DRAW target(cyan) with shift(wtx1,wty1) END IF END SUB SUB Flow(x0,y0) LET t = x0 LET x = y0 LET wx1 = w1Wndx(t) LET wy1 = w1Wndy(x) SET COLOR tclr FOR i= 1 to 5000 ! draw trajectory CALL PlotPoints(t,x) IF abs(wx1-wtx1)4 or abs(t)>4 then EXIT FOR NEXT i PLOT END SUB SUB FlowCheck(x0,y0) LET t = x0 LET x = y0 LET wx1 = w1Wndx(t) LET wy1 = w1Wndy(x) FOR i= 1 to 5000 ! draw trajectory IF abs(wx1-wtx1)4 or abs(t)>4 then EXIT FOR END IF NEXT i END SUB SUB HitFlow(x0,y0,flowclr) LET t = x0 LET x = y0 LET wx1 = w1Wndx(t) LET wy1 = w1Wndy(x) SET COLOR flowclr FOR i= 1 to 5000 ! draw trajectory CALL PlotPoints(t,x) LET err = 0 LET oldx = x LET oldt = t LET oldwx= wx1 LET oldwy= wy1 CALL AdjustStep(oldwx,oldwy) CALL RungeKutta4(t,x) IF err=1 then EXIT FOR LET wx1= round(w1Wndx(t)) LET wy1= round(w1Wndy(x)) GET MOUSE: mx,my,s IF (abs(dx)4 or abs(t)>4 or s=2 then EXIT FOR END IF NEXT i PLOT END SUB SUB RungeKutta4(x,y) LET tstp= dt CALL FetchEq(eq,x,y,dx,dy) LET dx1= dx LET dy1= dy LET x1 = x + .5*dx1*tstp LET y1 = y + .5*dy1*tstp CALL FetchEq(eq,x1,y1,dx,dy) LET dx2= dx LET dy2= dy LET x2 = x + .5*dx2*tstp LET y2 = y + .5*dy2*tstp CALL FetchEq(eq,x2,y2,dx,dy) LET dx3= dx LET dy3= dy LET x3 = x + dx3*tstp LET y3 = y + dy3*tstp CALL FetchEq(eq,x3,y3,dx,dy) LET dx4= dx LET dy4= dy LET dy = (dy1 + 2*dy2 + 2*dy3 + dy4) / 6 LET dx = (dx1 + 2*dx2 + 2*dx3 + dx4) / 6 LET y = y + tstp*dy LET x = x + tstp*dx END SUB SUB AdjustStep(oldwx,oldwy) FOR i= 1 to 6 LET tx = x LET tt = t CALL RungeKutta4(tt,tx) LET wx1= w1Wndx(tt) LET wy1= w1Wndy(tx) LET pixels= max(abs(abs(wx1)-abs(oldwx)),abs(abs(wy1)-abs(oldwy))) IF pixels>=1 and pixels<2 then EXIT FOR ELSE IF pixels<1 then LET dt= dt*2 ELSE IF pixels>3 then LET dt= dt/2 END IF END IF NEXT i END SUB SUB FlowT(x0,y0) LET t= x0 LET x= y0 LET oldwy= 999 SET COLOR yellow FOR i= 1 to 5000 ! draw trajectory LET wx1= round(w1Wndx(t)) LET wy1= round(w1Wndy(x)) IF eq>4 and abs(wx1-w1x0)<2 and abs(wy1-w1y0)<2 then LET wx1= round(w1Wndx(0)) LET wy1= round(w1Wndy(0)) CALL PlotPoints(0,0) EXIT FOR END IF CALL PlotPoints(t,x) LET olddy= dy LET err= 0 CALL RungeKutta4(t,x) IF err=1 then EXIT FOR GET MOUSE: mx,my,s IF (abs(dx)4 or abs(t)>4 or s=2 then EXIT FOR END IF NEXT i PLOT END SUB SUB PlotPoints(t,x) IF wx1>=w1Lft and wx1<=w1Rgt and wy1>=w1Top and wy1<=w1Bas then PLOT wx1,wy1; ELSE PLOT END IF END SUB SUB VectorField CALL Plane1Refresh IF fieldState=1 then SET COLOR fieldclr ! vector field !LET dx= 1 ! aspect LET vstp= 20 LET hstp= vstp/2 FOR wy= w1top+hstp to w1bas-hstp step vstp FOR wx= w1Lft+hstp to w1Rgt-hstp step vstp LET x= w1Fncx(wx) LET y= w1Fncy(wy) IF round(y,2)<>0 then CALL FetchEq(eq,x,y,dx,dy) LET ang= angle(dx,-dy) LET sa = 4*sin(ang) LET wy1= wy+sa LET wy2= wy-sa LET ca = 4*cos(ang) PLOT wx-ca,wy2; wx+ca,wy1 !PLOT wx,wy; wx+ca,wy+sa END IF NEXT wx NEXT wy IF targstate=1 then DRAW target(cyan) with shift(wtx1,wty1) END IF BOX KEEP w1lft-5,w1rgt+5,w1bas+5,w1top-5 in w1GraphLayer$ ELSE IF targstate=1 then DRAW target(cyan) with shift(wtx1,wty1) END IF BOX KEEP w1lft-5,w1rgt+5,w1bas+5,w1top-5 in w1GraphLayer$ END IF END SUB SUB ClearAndReset CALL Plane1Refresh CALL VectorField IF targstate=1 then DRAW target(cyan) with shift(wtx1,wty1) BOX KEEP w1Lft-5,w1Rgt+5,w1Bas+5,w1Top-5 in w1GraphLayer$ END IF END SUB ! --- SUB Hit BOX KEEP wtx1-5,wtx1+5,wty1+5,wty1-5 in temp$ SET COLOR yellow BOX AREA wtx1-3,wtx1+3,wty1+3,wty1-3 CALL Delay(1/8) BOX SHOW temp$ at wtx1-5,wty1+5 PLOT LET hitflag= 1 END SUB ! --- PICTURE target(clr) SET COLOR clr PLOT -3, 0 PLOT -2,-1; -2,1 PLOT -1,-2; -1,2 PLOT 0,-3; 0,3 PLOT 1,-2; 1,2 PLOT 2,-1; 2,1 PLOT 3, 0 PLOT -4,0; 4,0 PLOT 0,-4; 0,4 END PICTURE PICTURE bigtarget(clr) !CALL SetLineWeight(3) SET COLOR clr PLOT -4,-1; 2,-1 PLOT -1,-4; -1,2 !CALL SetLineWeight(1) END PICTURE PICTURE vpointer PLOT 0, 0; 8,0 END PICTURE PICTURE VectPnt PLOT -3,-3; 0,0; -3,3 END PICTURE ! --- SUB Buttons CALL SetTextFont(1,9,"bold") LET bwid= 88 LET bhgt= 17 LET cclft= w1Lft LET ccrgt= cclft + bwid LET cctop= w1Bas + 30 LET ccbas= cctop + bhgt CALL DrawButton(cclft,ccrgt,ccbas,cctop,5,"Clear Curves") LET ctlft= ccrgt+2 LET ctrgt= ctlft+bwid LET cttop= cctop LET ctbas= cttop+bhgt CALL DrawButton(ctlft,ctrgt,ctbas,cttop,5,"Set Target") ! LET sllft= ctrgt+2 ! LET slrgt= sllft+bwid ! LET sltop= cctop ! LET slbas= sltop + bhgt ! CALL DrawButton(sllft,slrgt,slbas,sltop,5,"Draw Field") ! LET glft= sllft ! LET grgt= slrgt ! LET gtop= slBas+2 ! LET gbas= round(gTop+bhgt/2) ! LET gcnt= 4 ! LET gwid= round((slRgt-slLft)/4) ! FOR btn= 0 to gcnt-1 ! LET l= glft + btn*gwid ! LET r= l + gwid-2 ! CALL DrawButton(l,r,gbas,gtop,0,"") ! SET COLOR btn+2 ! BOX AREA l+1,r-1,gbas-1,gtop+1 ! NEXT btn END SUB END SUB