! I've realized that there's a slight variation we should play on poles and vibes. ! ! I guess I should have added Poles and Vibes to the list I just sent you ! and Michael, because what we have is excellent and this new thing will take a ! little effort to get right. (Remember, I wanted to kill the a and b ! as axis names.) ! ! What is missing is phase. An important lesson is that adding two sinusoidal ! waves with the same frequency can give zero back - if they are perfectly out ! of phase. ! ! Phase doesn't much matter if there is only one damped sine around - it just ! shifts the picture one way or another. But if thre are two they have a relative ! phase, which shows up in the pictures. ! ! One way to deal would be to have a phase button somewhere in the "w" control ! area, which would bring up a phase slider giving phi between -pi and pi, ! and replacing Bcos(dt) with B cos(dt - phi). The B then becomes the ! amplitude. Whenever I am in this situation I also look for a little polar ! diagram with angle phi and magnitude B, but this may get too confusing so ! maybe we should skip it at first. ! ! The phase has nothing to do with the pole diagram, you can't see it from that ! picture. !Poles and Vibrations: Better! - one window, using cos instead of sin.... ! !The window seems to decide which poles to move by taking the one with the !nearest imaginary part. ! !When we take say B = 0, the red poles should disappear along with the sliders !for their position and their value should also vanish. Sometimes they do (but !the poles are still left on the window) - this part seems a little balky. ! !There is a mysterious pair of slider markers between b and d.... !! File: PolesVibrations.mit !! July 28, 2001 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 PUBLIC largefonts PUBLIC title$ LET toolHgt= 540 LET toolWid= 780 LET window$= "MIT Visual Math Project" LET colorscheme= 0 LET title$ = "Poles and Vibrations" SUB ThisProgram CALL PolesAndVibrations 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= SetTextFont(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+20,toolhdr-10,txt$,titleclr) END SUB SUB CopyRight(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 DATA 10001000010001000000010000000000 DATA 10001000010001000000010000000000 DATA 10001000010001001110010110010110 DATA 11111000011111010001011001011001 DATA 10001000010001010001010001010001 DATA 10001000010001010001010001010001 DATA 10001010010001001110010001010001 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 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,.2,.4 ! Back ELSE SET COLOR MIX(0) b4,b4,b4 ! Back SET COLOR MIX(0) .7,.7,.7 ! 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 SET COLOR MIX( 7) b5,b2,b0 ! 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) b2,b2,b5 ! blue SET COLOR MIX(12) b5,b0,b5 ! magenta SET COLOR MIX(13) b5,b3,b0 ! magenta ELSE SET COLOR MIX( 7) b3,b1,b0 ! red SET COLOR MIX( 8) b4,b3,b0 ! yellow SET COLOR MIX( 9) b0,b2,b0 ! green SET COLOR MIX(10) b0,b3,b3 ! cyan SET COLOR MIX(11) b2,b0,b4 ! blue SET COLOR MIX(12) b3,b0,b3 ! magenta SET COLOR MIX(13) b5,b3,b0 ! magenta END IF FOR i= 1 to 13 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 = drkmid LET axisclr = litgry LET axislabelclr= white LET titleclr = litgry LET rightsclr = litmif LET numberlineclr= litgry LET slotdrkclr= black LET slotlgtclr= litmid ELSE LET planeclr= white LET gridclr = litgry LET rimclr = drkmid LET axisclr = drkgry LET axislabelclr= black LET titleclr= drkgry LET rightsclr= drkgry LET numberlineclr= drkgry LET slotdrkclr= drkgry LET slotlgtclr= white 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 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 CALL CopyRight(iLft+20,iBas-3,rightsclr) 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 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= round(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= round(fHgt/ystp) FOR i= 0 to lines LET n = fBas + 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,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 PLOT wx0,wBas; wx0,wTop-3 DRAW arrow3 with rotate(-pi/2) * shift(wx0,wTop-3) END IF END SUB SUB EdgesVrt(wAxis,wBas,wTop,fBas,fTop,first,stp1,stp2,nstp) 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,-1) CALL VNumberLineLabels(wBas,wTop,wAxis,fBas,fTop,first,stp1,stp2,nstp,-1) END SUB SUB EdgesHrz(wLft,wRgt,wBas,fLft,fRgt,first,stp1,stp2,nstp) 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,1) CALL HNumberLineLabels(wLft,wRgt,wBas,fLft,fRgt,first,stp1,stp2,nstp,1) END SUB SUB EdgesHrzPi(wLft,wRgt,wBas,fLft,fRgt,first,stp1,stp2,nstp) 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,1) CALL HNumberLineLabelsPi(wLft,wRgt,wBas,fLft,fRgt,first,stp1,stp2,nstp,1) 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 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= 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= 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 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(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 ! short ticks LET ticks= int(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(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 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-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) !LET form$ = "-%.#" END IF CALL SetTextFont(1,9,"normal") LET ticks= round(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" ELSE LET n$= trim$(using$(form$,abs(n))) 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= round(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 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= 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 hwndx2(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 hwndx3(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 hwndx4(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 END MODULE MODULE GraphParts ! ---- visual parts PICTURE arrow2 PLOT 2,0; 0, 2; 0,-2; 2,0 PLOT 0,2; 0,-2 PLOT 1,1; 1,-1 END PICTURE PICTURE arrow3 PLOT 3, 0; 0,-3; 0,3; 3,0 PLOT 1,-2; 1, 2 PLOT 2,-1; 2, 1 END PICTURE PICTURE arrow4 PLOT -2,-2; -2,2; 0,0; -2,-2 PLOT -3,-3; -3,3; 0,0; -3,-3 PLOT -4,-4; -4,4; 0,0; -4,-4 END PICTURE PICTURE pi9 PLOT 0,-4; 1,-5; 5,-5; 6,-6 PLOT 1, 0; 2,-1; 2,-5 PLOT 4, 0; 5,-1; 5,-5 END PICTURE PICTURE Diamond PLOT 0,2; 0,-2 PLOT -1,1; -1,-1 PLOT -2,0 PLOT 1,1; 1,-1 PLOT 2,0 END PICTURE PICTURE diamond5(clr) SET COLOR clr PLOT -2, 0 PLOT -1,-1; -1,1 PLOT -0,-2; -0,2 PLOT 1,-1; 1,1 PLOT 2, 0 SET COLOR black PLOT -3,0; 0,3; 3,0; 0,-3; -3,0 END PICTURE PICTURE diamond7(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 ! SET COLOR black ! PLOT -4,0; 0,4; 4,0; 0,-4; -4,0 END PICTURE PICTURE Root(rWid) LET Lft= -rWid PLOT Lft+1,-12; 0,-12 PLOT Lft ,-11; 0,-11 PLOT Lft+1,-12; Lft-5,0; Lft-7,-4 PLOT Lft+2,-12; Lft-4,0; Lft-6,-4 END PICTURE PICTURE IntegralSign(clr) SET COLOR clr PLOT 0, -9; 0, 2 PLOT 1, -9; 1, 2 PLOT 1,-10; 2,-10 PLOT -1, 3; 0, 3 END PICTURE PICTURE Summation(clr) SET COLOR clr PLOT 0, 2; 7, 2 PLOT 0, 2; 6, -4; 0,-10 PLOT 1, 2; 7, -4; 1,-10 PLOT 0,-10; 7,-10 END PICTURE PICTURE PlusMinus(clr) SET COLOR clr PLOT -3, 0; 3, 0 PLOT -3,-6; 3,-6 PLOT 0,-9; 0,-3 END PICTURE END MODULE MODULE InterAction 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 SUB DrawButton(l,r,b,t,tBas,t$) SET COLOR white BOX LINES l,r-1,b-1,t SET COLOR black BOX LINES l+1,r,b,t+1 SET COLOR drkgry BOX AREA l+1,r-1,b-1,t+1 IF t$<>"" 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 ! *** SUB PolesAndVibrations 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,slideclr DECLARE PUBLIC workLft,workRgt,workBas,workTop,workMidx DECLARE PUBLIC qLft,qRgt,qBas,qTop DECLARE PUBLIC infLft,infRgt,infBas,infTop PUBLIC h1a,h2b,h1olda,h2oldb PUBLIC h3c,h4d,h3oldc,h4oldd PUBLIC h6j,h5k,h6oldj,h5oldk ! --- help screen array --- DIM info$(1:1) MAT READ info$ DATA "Information on poles and vibrations" ! ---------- Utility functions --- DEF clamp(n,lo,hi)= min(max(n,lo),hi) DEF roundn(n,step)= round(n/step)*step DEF e= exp(1) ! --- colors --- LET abclr= green LET cdclr= red ! ---------- Graphing plane parameters and methods ---------- ! --- 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 LET w1Flag= 1 LET w1Lft= workRgt-210 ! pixel bounds LET w1Rgt= w1Lft+160 LET w1Top= workBas-190 LET w1Bas= w1Top+160 LET w1fLft= -0.4 ! function bounds LET w1fRgt= 0.4 LET w1fTop= 8 LET w1fBas= -8 LET w1Xax$= "" ! axis labels LET w1Yax$= "" LET w1xGridstep= 0.1 ! horizontal grid intervals LET w1yGridstep= 1.0 ! vertical grid intervals LET w1SxTik = 0.1 ! horizontal axis Tik marks LET w1LxTik = 0.2 LET w1xLabel= 0.2 LET w1Firstx= w1fLft LET w1SyTik = 1 ! vertical axis Tik marks LET w1LyTik = 4 LET w1yLabel= 4 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 Plane1Redraw CALL Plane1Init CALL Plane1Update CALL t1Init CALL t5Set CALL h1DrawSlider(h1a) CALL h2DrawSlider(h2b) END SUB SUB Plane1Init CALL DrawPlane1(1,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$ END SUB SUB Plane1Refresh BOX SHOW w1gridLayer$ at w1Lft-5,w1Bas+5 END SUB ! ------------------------------------------ ! --- plane 2 data --- DECLARE PUBLIC w2Lft,w2Rgt,w2Bas,w2Top,w2Midx,w2Midy DECLARE PUBLIC w2fLft,w2fRgt,w2fBas,w2fTop,w2x0,w2y0 DECLARE PUBLIC w2Firstx, w2SxTik, w2LxTik, w2xLabel, w2xGridstep DECLARE PUBLIC w2Firsty, w2SyTik, w2LyTik, w2yLabel, w2yGridstep DECLARE PUBLIC w2xPiFlag LET w2Flag= 1 LET w2Lft = workMidx-336 ! pixel bounds LET w2Rgt = w2Lft+672 LET w2Top = workTop+35 !LET w2Bas = w2Top+192 LET w2Bas = w2Top+224 LET w2fLft= 0 ! function bounds * pi LET w2fRgt= 8 LET w2fTop= 8 LET w2fBas= -8 LET w2xGridstep= 0.5 ! horizontal grid intervals LET w2yGridstep= 1 ! vertical grid intervals LET w2xAx$= "t" ! axis labels LET w2yAx$= "x" LET w2SxTik = 0.25 ! horizontal axis Tik marks LET w2LxTik = 0.50 LET w2xLabel= 1 LET w2Firstx= w2fLft LET w2xPiFlag= 1 LET w2SyTik = 1 ! vertical axis Tik marks LET w2LyTik = 2 LET w2yLabel= 2 LET w2Firsty= w2fBas ! --- plane 2 methods --- DECLARE DEF w2fncx, w2fncy, w2wndx, w2wndy ! window/function transforms CALL Plane2Variables SUB Plane2Init CALL DrawPlane2(1,1,1) ! grid, axes, zeroaxes CALL Grid1 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 Grid1 LOCAL i,wx SET COLOR drkmid FOR i= 2 to 6 step 2 LET wx= w2wndx(i) PLOT wx,w2Top+1; wx,w2Bas-1 NEXT i END SUB ! ------------------------------------------ ! --- plane 3 data --- DECLARE PUBLIC w3Lft,w3Rgt,w3Bas,w3Top,w3Midx,w3Midy DECLARE PUBLIC w3fLft,w3fRgt,w3fBas,w3fTop,w3x0,w3y0 DECLARE PUBLIC w3Firstx, w3SxTik, w3LxTik, w3xLabel, w3xGridstep DECLARE PUBLIC w3Firsty, w3SyTik, w3LyTik, w3yLabel, w3yGridstep LET w3Flag= 1 ! plane is turned on or off LET w3Lft = workMidx+46 ! pixel bounds LET w3Rgt = w3Lft+128 LET w3Top = w1Top LET w3Bas = w1Bas LET w3fLft= -0.4 ! function bounds LET w3fRgt= 0.4 LET w3fTop= 8 LET w3fBas= -8 LET w3xAx$= "c" ! axis labels LET w3yAx$= "d" LET w3xGridstep= 0.1 ! grid line intervals LET w3yGridstep= 1 LET w3SxTik = 0.1 ! axis Tik marks LET w3LxTik = 0.2 LET w3xLabel= 0.2 LET w3Firstx= w3fLft LET w3SyTik = 1 LET w3LyTik = 4 LET w3yLabel= 4 LET w3Firsty= w3fBas ! --- plane 3 methods --- DECLARE DEF w3fncx,w3fncy,w3wndx,w3wndy ! window/function transforms CALL Plane3Variables SUB Plane3Clear BOX CLEAR w3Lft-20,w3Rgt+20,workBas-5,w3Top-20 END SUB SUB Plane3Redraw CALL Plane3Init CALL Plane3Update CALL t2init !CALL t3Init !CALL t4Init CALL t5Set CALL h3DrawSlider(h3c) CALL h4DrawSlider(h4d) END SUB SUB Plane3Init CALL DrawPlane3(1,1,1) ! grid, axes, zeroaxes CALL SetTextFont(1,12,"bold") CALL PlotTextLJ(w3Rgt+8,w3y0+3,w3xAx$,axislabelclr) ! axis labels CALL PlotTextCJ(w3x0,w3Top-10,w3yAx$,axislabelclr) BOX KEEP w3Lft-5,w3Rgt+5,w3Bas+5,w3Top-5 in w3gridLayer$ END SUB SUB Plane3Refresh BOX SHOW w1gridLayer$ at w1Lft-5,w1Bas+5 END SUB ! ---------- Slider parameters and methods ---------- ! ----------- vertical sliders ------------ LET slideHeight= 96 LET c1$= "A" LET c2$= "B" ! --- vertical slider 1 --- DECLARE DEF v1Fncy,v1wndy ! window/function transforms DECLARE PUBLIC v1axis,v1wLft,v1wRgt,v1wBas,v1wTop,v1sBas,v1sTop DECLARE PUBLIC v1fBas,v1fTop,v1First,v1STik,v1LTik,v1Label DECLARE PUBLIC v1name$,v1m,v1form$,v1clr LET v1clr = slideclr LET v1name$= "a" LET v1form$= "-%.##" LET v1m = 1 LET v1a = 0 LET v1olda = -999 LET v1axis = w1Lft - 36 LET v1wBas = w1Bas LET v1wTop = v1wBas-slideHeight LET v1fBas = -0.4 LET v1fTop = 0.4 LET v1STik = 0.1 LET v1LTik = 0.2 LET v1Label= 0.2 LET v1First= v1fBas CALL v1SliderVariables ! ----------- horizontal sliders ------------ LET slideWidth= 96 ! --- horizontal slider 1 --- DECLARE DEF h1Fncx ! window/function transforms DECLARE PUBLIC h1axis,h1wLft,h1wRgt,h1wBas,h1wTop,h1sLft,h1sRgt DECLARE PUBLIC h1fLft,h1fRgt,h1First,h1STik,h1LTik,h1Label DECLARE PUBLIC h1name$,h1m,h1form$,h1clr LET h1clr = abclr LET h1name$= "a" LET h1form$= "-%.##" LET h1m = 1 LET h1a = 0 LET h1olda = -999 LET h1axis = w2Bas + 180 LET h1wLft = w2Lft LET h1wRgt = h1wLft+slideWidth LET h1fLft = -0.4 LET h1fRgt = 0.4 LET h1STik = 0.1 ! short tick marks LET h1LTik = 0.2 ! long tick marks LET h1Label= 0.2 ! labels LET h1First= h1fLft ! first tick mark CALL h1SliderVariables ! --- Slider 2 --- DECLARE DEF h2Fncx ! window/function transforms DECLARE PUBLIC h2axis,h2wLft,h2wRgt,h2wBas,h2wTop,h2sLft,h2sRgt DECLARE PUBLIC h2fLft,h2fRgt,h2First,h2STik,h2LTik,h2Label DECLARE PUBLIC h2name$,h2m,h2form$,h2clr LET h2clr = abclr LET h2name$ = "b" LET h2form$ = "--%.##" LET h2m = 1 LET h2b = 1 LET h2oldb = -999 LET h2axis = h1axis+40 LET h2wLft = h1wLft LET h2wRgt = h2wLft+slidewidth LET h2fLft = -8 LET h2fRgt = 8 LET h2STik = 1 LET h2LTik = 4 LET h2Label= 4 LET h2First= h2fLft CALL HSlider2Variables ! --- slider 3 --- DECLARE DEF h3Fncx ! window/function transforms DECLARE PUBLIC h3axis,h3wLft,h3wRgt,h3wBas,h3wTop,h3sLft,h3sRgt DECLARE PUBLIC h3fLft,h3fRgt,h3First,h3STik,h3LTik,h3Label DECLARE PUBLIC h3name$,h3m,h3form$,h3clr LET h3clr = cdclr LET h3name$ = "c" LET h3form$ = "-%.##" LET h3m = 1 LET h3c = 0 LET h3oldc = -999 LET h3axis = h1axis LET h3wLft = h1wRgt + 100 LET h3wRgt = h3wLft+slideWidth LET h3fLft = -0.4 LET h3fRgt = 0.4 LET h3STik = 0.1 LET h3LTik = 0.2 LET h3Label= 0.2 LET h3First= h3fLft CALL HSlider3Variables ! --- slider 4 --- DECLARE DEF h4Fncx ! window/function transforms DECLARE PUBLIC h4axis,h4wLft,h4wRgt,h4wBas,h4wTop,h4sLft,h4sRgt DECLARE PUBLIC h4fLft,h4fRgt,h4First,h4STik,h4LTik,h4Label DECLARE PUBLIC h4name$,h4m,h4form$,h4clr LET h4clr = cdclr LET h4name$ = "d" LET h4form$ = "--%.##" LET h4m = 1 LET h4d = 1 LET h4oldd = -999 LET h4axis = h2axis LET h4wLft = h3wLft LET h4wRgt = h4wLft+slideWidth LET h4fLft = -8 LET h4fRgt = 8 LET h4STik = 1 LET h4LTik = 4 LET h4Label= 4 LET h4First= h4fLft CALL HSlider4Variables ! --- slider 5 --- DECLARE DEF h5Fncx ! window/function transforms DECLARE PUBLIC h5axis,h5wLft,h5wRgt,h5wBas,h5wTop,h5sLft,h5sRgt DECLARE PUBLIC h5fLft,h5fRgt,h5First,h5STik,h5LTik,h5Label DECLARE PUBLIC h5name$,h5m,h5form$,h5clr LET h5clr = slideclr LET h5name$ = c2$ ! B LET h5form$ = "--%.##" LET h5m = 1 LET h5k = 1 LET h5oldk = -999 LET h5axis = w1Top+30 LET h5wLft = h3wLft LET h5wRgt = h3wRgt LET h5fLft = 0 LET h5fRgt = 2 LET h5STik = 0 LET h5LTik = 0.5 LET h5Label= 0.5 LET h5First= h5fLft CALL HSlider5Variables ! --- slider 6 --- DECLARE DEF h6Fncx ! window/function transforms DECLARE PUBLIC h6axis,h6wLft,h6wRgt,h6wBas,h6wTop,h6sLft,h6sRgt DECLARE PUBLIC h6fLft,h6fRgt,h6First,h6STik,h6LTik,h6Label DECLARE PUBLIC h6name$,h6m,h6form$,h6clr LET h6clr = slideclr LET h6name$ = c1$ ! A LET h6form$ = "--%.##" LET h6m = 1 LET h6j = 1 LET h6oldj = -999 LET h6axis = h5axis LET h6wLft = h1wLft LET h6wRgt = h1WRgt LET h6fLft = 0 LET h6fRgt = 2 LET h6STik = 0 LET h6LTik = 0.5 LET h6Label= 0.5 LET h6First= h6fLft CALL HSlider6Variables ! ---------- Text Output Rects ---------- ! --- text rectangle 1 --- LET t1BasLn= h1wTop - 30 LET t1Lft = h1wLft + 23 LET t1Rgt = h1wRgt + 30 LET t1Bas = t1BasLn + 5 LET t1Top = t1BasLn - 15 SUB t1Label CALL SuperSubScriptRJ(t1Lft,t1BasLn,"z = ",abclr) END SUB SUB t1Set LOCAL r$,i$,z$ CALL t1Clear LET r$ = trim$(using$("--%.##",h1a)) LET i$ = trim$(using$("--%.##i",abs(h2b))) CALL StringWidth(r$,rw) LET pmx= t1lft + rw + 10 LET ix = pmx + 3 + 6 CALL PlotTextLJ(t1Lft,t1BasLn,r$,abclr) DRAW PlusMinus(abclr) with shift(pmx,t2BasLn) CALL PlotTextLJ(ix,t1BasLn,i$,abclr) ! LET z$ = r$ & " +- " & i$ ! CALL PlotTextLJ(t1Lft,t1BasLn,z$,abclr) END SUB SUB t1Clear BOX CLEAR t1Lft-2,t1Rgt,t1Bas,t1Top END SUB SUB t1Init CALL t1Label CALL t1Set END SUB ! --- text rectangle 2 --- LET t2BasLn= t1BasLn LET t2Lft = h3wLft + 23 LET t2Rgt = h3wRgt + 30 LET t2Bas = t2BasLn + 5 LET t2Top = t2BasLn - 15 SUB t2Label CALL SuperSubScriptRJ(t2Lft,t2BasLn,"w = ",cdclr) END SUB SUB t2Set LOCAL r$,i$,z$ CALL t2Clear LET r$ = trim$(using$("--%.##",h3c)) LET i$ = trim$(using$("--%.##i",abs(h4d))) CALL StringWidth(r$,rw) LET pmx= t2lft + rw + 10 LET ix = pmx + 6 + 3 CALL PlotTextLJ(t2Lft,t2BasLn,r$,cdclr) DRAW PlusMinus(cdclr) with shift(pmx,t2BasLn) CALL PlotTextLJ(ix,t2BasLn,i$,cdclr) END SUB SUB t2Clear BOX CLEAR t2Lft-2,t2Rgt,t2Bas,t2Top END SUB SUB t2Init CALL t2Label CALL t2Set END SUB ! --- text rectangle 3 --- LET t3BasLn= t1BasLn LET t3Lft = h3wLft + 23 LET t3Rgt = h3wRgt + 30 LET t3Bas = t3BasLn + 5 LET t3Top = t3BasLn - 15 SUB t3Label CALL SuperSubScriptRJ(t3Lft,t3BasLn,"z_[3] = ",cclr) END SUB SUB t3Set LOCAL r$,i$,z$ CALL t3Clear LET r$ = trim$(using$("--%.##",h3c)) LET i$ = trim$(using$("--%.##i",abs(h4d))) LET z$ = r$ & " + " & i$ CALL PlotTextLJ(t3Lft,t3BasLn,z$,cclr) END SUB SUB t3Clear BOX CLEAR t3Lft-2,t3Rgt,t3Bas,t3Top END SUB SUB t3Init CALL t3Label CALL t3Set END SUB ! ! ! ! --- text rectangle 4 --- ! ! LET t4BasLn= t3BasLn + 20 ! LET t4Lft = t3Lft ! LET t4Rgt = t3Rgt ! LET t4Bas = t4BasLn + 5 ! LET t4Top = t4BasLn - 15 ! ! SUB t4Label ! CALL SuperSubScriptRJ(t4Lft,t4BasLn,"z_[4] = ",cclr) ! END SUB ! ! SUB t4Set ! LOCAL r$,i$,z$ ! CALL t4Clear ! LET r$ = trim$(using$("--%.##",h3c)) ! LET i$ = trim$(using$("--%.##i",abs(h4d))) ! LET z$ = r$ & " + " & i$ ! CALL PlotTextLJ(t4Lft,t4BasLn,z$,cclr) ! END SUB ! ! SUB t4Clear ! BOX CLEAR t4Lft-2,t4Rgt,t4Bas,t4Top ! END SUB ! ! SUB t4Init ! CALL t4Label ! CALL t4Set ! END SUB ! --- reset pairs --- SUB ts12 CALL t1Set !CALL t2Set END SUB SUB ts34 CALL t2Set !CALL t4Set END SUB ! --- Function Names --- ! --- text rectangle 5 --- LET t5BasLn= w2Bas + 43 LET t5Lft = w2Lft LET t5Lfte = w2Lft + 240 LET t5Rgt = w2Rgt - 150 LET t5Bas = t5BasLn + 5 LET t5Top = t5BasLn - 15 SUB GetFunctionText IF h5k>0 and h6j>0 then LET f1$= "f(t) = " & c1$ & "e^[at]cos(bt) + " & c2$ & "e^[ct]cos(dt)" !LET f2$= "g(t) = " & c1$ & "e^[at] + " & c2$ & "e^[ct]" IF h1a>h3c then LET f2$= "g(t) = " & c1$ & "e^[at]" ELSE IF h1a0 then LET f1$= "f(t) = " & c1$ & "e^[at]cos(bt)" LET f2$= "g(t) = " & c1$ & "e^[at]" ELSE IF h5k>0 then LET f1$= "f(t) = " & c2$ & "e^[ct]cos(dt)" LET f2$= "g(t) = " & c2$ & "e^[ct]" ELSE LET f1$= "f(t) = 0" LET f2$= "g(t) = 0" END IF END SUB ! IF h1a>h3c then ! LET y= h6j*e^(h1a*x) ! ELSE IF h1ah6j then LET f2$= "g(t) = Be^[ct]" ELSE LET f2$= "g(t) = Ae^[ct]" END IF ! IF h5k>0 and h6j>0 then ! IF h1a>h3c then ! LET f2$= "g(t) = " & c1$ & "e^[at]" ! ELSE IF h1a0 then ! LET f2$= "g(t) = " & c1$ & "e^[at]" ! ELSE IF h5k>0 then ! LET f2$= "g(t) = " & c2$ & "e^[ct]" ! ELSE ! LET f2$= "g(t) = 0" ! END IF END SUB ! DEF envl(x)= h6j*e^(h1a*x) + h5k*e^(h3c*x) ! --- functions, parameters, equations, constants --- DEF envl(x)= h6j*e^(h1a*x) + h5k*e^(h3c*x) DEF wave(x)= h6j*e^(h1a*x)*cos(h2b*x) + h5k*e^(h3c*x)*cos(h4d*x) ! --- Draw the screen --- LET h1a= -0.1 LET h2b= 2 LET h3c= 0.1 LET h4d= 4 LET h5k= 0 LET h5oldk= 0 LET h6j,h6oldj= 1 CALL InitScreen SUB InitScreen BOX CLEAR worklft,workrgt,workbas,worktop IF colorscheme=1 then LET aclr= red ! envelope LET bclr= blue ! wave LET cclr= green ! complex points LET slideclr= drkgry ! slider names ELSE LET aclr= yellow ! envelope LET bclr= cyan ! wave LET cclr= green ! complex points LET slideclr= litgry ! slider names END IF LET h5clr,h6clr= slideclr LET abclr= green LET cdclr= red CALL Plane1Init CALL h6DrawSlider(h6j) CALL h1DrawSlider(h1a) CALL h2DrawSlider(h2b) CALL t1Init LET w1Flag=1 ! CALL t2Init ! LET h6j= 0 ! CALL Plane3Init ! CALL t4Init CALL h5DrawSlider(h5k) IF b>0 then CALL h3DrawSlider(h3c) CALL h4DrawSlider(h4d) CALL t2Init LET w3Flag= 1 ELSE BOX CLEAR h3wLft-20,h3wRgt+50,h4wBas+10,t2Top-10 LET w3Flag= 0 END IF CALL Plane2Init CALL t5Init CALL UpdateGraphs CALL Plane1Update ! CALL v1DrawSlider(v1a) END SUB ! ----------------- Event manager ----------------- DO CALL MouseDown(mx,my,ms) IF mx>w1Lft and mxw1Top and myw3Lft and mxw3Top and myv1wLft and mxv1sTop and myh1sLft and mxh1wTop and myh2sLft and mxh2wTop and myh3sLft and mxh3wTop and myh4sLft and mxh4wTop and myh5sLft and mxh5wTop and myh6sLft and mxh6wTop and myinfLft and mxinfTop and myqLft and mxqTop and myh1olda or h2b<>h2oldb then CALL UpdateGraphs CALL t1Set CALL h1Mark(h1a) CALL h2mark(h2b) LET h1olda= h1a LET h2oldb= h2b CALL t5SetG END IF GET MOUSE: mx,my,ms LET mx = clamp(mx,w1Lft,w1Rgt) LET my = clamp(my,w1Top,w1Bas) LOOP until ms=3 ELSE IF w1Flag=0 then ! a is off DO LET h3c= round(w1fncx(clamp(mx,w1Lft,w1Rgt)),2) LET h4d= round(w1fncy(clamp(my,w1Top,w1Bas)),1) IF h3c<>h3oldc or h4d<>h4oldd then CALL UpdateGraphs CALL t2Set CALL h3Mark(h3c) CALL h4mark(h4d) LET h3oldc= h3c LET h4oldd= h4d CALL t5SetG END IF GET MOUSE: mx,my,ms LET mx = clamp(mx,w1Lft,w1Rgt) LET my = clamp(my,w1Top,w1Bas) LOOP until ms=3 ELSE LET wa = w1Wndx( h1a) LET wb1= w1Wndy( h2b) LET wb2= w1Wndy(-h2b) LET wc = w1Wndx( h3c) LET wd1= w1Wndy( h4d) LET wd2= w1Wndy(-h4d) LET abd1= sqr( (mx-wa)^2 + (my-wb1)^2 ) LET abd2= sqr( (mx-wa)^2 + (my-wb2)^2 ) LET cdd1= sqr( (mx-wc)^2 + (my-wd1)^2 ) LET cdd2= sqr( (mx-wc)^2 + (my-wd2)^2 ) IF (abd1h1olda or h2b<>h2oldb then CALL UpdateGraphs CALL t1Set CALL h1Mark(h1a) CALL h2mark(h2b) LET h1olda= h1a LET h2oldb= h2b CALL t5SetG END IF GET MOUSE: mx,my,ms LET mx = clamp(mx,w1Lft,w1Rgt) LET my = clamp(my,w1Top,w1Bas) LOOP until ms=3 ELSE DO LET h3c= round(w1fncx(clamp(mx,w1Lft,w1Rgt)),2) LET h4d= round(w1fncy(clamp(my,w1Top,w1Bas)),1) IF h3c<>h3oldc or h4d<>h4oldd then CALL UpdateGraphs CALL t2Set CALL h3Mark(h3c) CALL h4mark(h4d) LET h3oldc= h3c LET h4oldd= h4d CALL t5SetG END IF GET MOUSE: mx,my,ms LET mx = clamp(mx,w1Lft,w1Rgt) LET my = clamp(my,w1Top,w1Bas) LOOP until ms=3 END IF END IF END SUB SUB Plane3MouseDrag(mx,my) ! DO ! GET MOUSE: mx,my,ms ! LET h3c= round(w3fncx(clamp(mx,w3Lft,w3Rgt)),2) ! LET h4d= round(w3fncy(clamp(my,w3Top,w3Bas)),1) ! IF h3c<>h3oldc or h4d<>h4oldd then ! CALL UpdateGraphs ! CALL ts34 ! CALL h3mark(h3c) ! CALL h4mark(h4d) ! LET h3oldc= h3c ! LET h4oldd= h4d ! CALL T5SetG ! END IF ! LOOP until ms=3 END SUB ! --- Read Vertical Sliders SUB v1MouseClick(mx,my) CALL MouseUp(mx,my,ms) LET v1a= v1Fncy(clamp(my,v1wTop,v1wBas)) LET v1a= roundn(v1a,v1STik) IF v1a<>v1olda then CALL v1Action LET v1olda= v1a END IF END SUB SUB v1MouseDrag(mx,my) DO GET MOUSE: mx,my,ms LET v1a= round(v1Fncy(clamp(my,v1wTop,v1wBas)),2) IF v1a<>v1olda then CALL v1Action LET v1olda= v1a END IF LOOP until ms=3 END SUB SUB v1Action CALL v1Mark(v1a) ! CALL UpdateGraphs ! CALL ts12 END SUB ! --- Read Slider 1 --- SUB h1MouseClick(mx,my) CALL MouseUp(mx,my,ms) LET h1a= h1Fncx(clamp(mx,h1wLft,h1wRgt)) LET h1a= roundn(h1a,h1STik) IF h1a<>h1olda then CALL h1Action LET h1olda= h1a END IF END SUB SUB h1MouseDrag(mx,my) DO GET MOUSE: mx,my,ms LET h1a= round(h1Fncx(clamp(mx,h1wLft,h1wRgt)),2) IF h1a<>h1olda then CALL h1Action LET h1olda= h1a END IF LOOP until ms=3 END SUB SUB h1Action CALL h1Mark(h1a) CALL UpdateGraphs CALL ts12 CALL T5SetG END SUB ! --- Read Slider 2 --- SUB h2MouseClick(mx,my) CALL MouseUp(mx,my,ms) LET h2b= h2Fncx(clamp(mx,h2wLft,h2wRgt)) LET h2b= roundn(h2b,h2STik) IF h2b<>h2oldb then CALL h2Action LET h2oldb= h2b END IF END SUB SUB h2MouseDrag(mx,my) CALL Plane2Refresh DO GET MOUSE: mx,my,ms LET h2b= round(h2Fncx(clamp(mx,h2wLft,h2wRgt)),2) IF h2b<>h2oldb then CALL h2Action LET h2oldb= h2b END IF LOOP until ms=3 END SUB SUB h2Action CALL h2mark(h2b) BOX SHOW Env$ at w2Lft,w2Bas CALL Plane2Wave CALL Plane1Refresh CALL Plane1Update CALL ts12 END SUB ! --- Read Slider 3 --- SUB h3MouseClick(mx,my) CALL MouseUp(mx,my,ms) LET h3c= h3Fncx(clamp(mx,h3wLft,h3wRgt)) LET h3c= roundn(h3c,h3STik) IF h3c<>h3oldc then CALL h3Action LET h3oldc= h3c END IF END SUB SUB h3MouseDrag(mx,my) DO GET MOUSE: mx,my,ms LET h3c= round(h3Fncx(clamp(mx,h3wLft,h3wRgt)),2) IF h3c<>h3oldc then CALL h3Action LET h3oldc= h3c END IF LOOP until ms=3 END SUB SUB h3Action CALL h3mark(h3c) CALL UpdateGraphs CALL ts34 CALL T5SetG END SUB ! --- Read Slider 4 --- SUB h4MouseClick(mx,my) CALL MouseUp(mx,my,ms) LET h4d= h4Fncx(clamp(mx,h4wLft,h4wRgt)) LET h4d= roundn(h4d,h4STik) IF h4d<>h4oldd then CALL h4Action LET h4oldd= h4d END IF END SUB SUB h4MouseDrag(mx,my) CALL Plane2Refresh DO GET MOUSE: mx,my,ms LET h4d= round(h4Fncx(clamp(mx,h4wLft,h4wRgt)),2) IF h4d<>h4oldd then CALL h4Action LET h4oldd= h4d END IF LOOP until ms=3 END SUB SUB h4Action CALL h4mark(h4d) BOX SHOW Env$ at w2Lft,w2Bas CALL Plane2Wave CALL Plane3Refresh CALL Plane3Update CALL ts34 END SUB ! --- Read Slider 5 --- SUB h5MouseClick(mx,my) CALL MouseUp(mx,my,ms) LET h5k= h5Fncx(clamp(mx,h5wLft,h5wRgt)) LET h5k= roundn(h5k,h5LTik) IF h5k<>h5oldk then CALL h5Action LET h5oldk= h5k END IF END SUB SUB h5MouseDrag(mx,my) CALL Plane2Refresh DO GET MOUSE: mx,my,ms LET h5k= round(h5Fncx(clamp(mx,h5wLft,h5wRgt)),2) IF h5k<>h5oldk then CALL h5Action LET h5oldk= h5k END IF LOOP until ms=3 END SUB SUB h5Action CALL h5mark(h5k) IF h5k=0 and h5oldk<>0 then BOX CLEAR h3wLft-20,h3wRgt+50,h4wBas+10,t2Top-10 LET w3Flag= 0 !CALL Plane1Update CALL UpdateGraphs CALL t5Set ELSE IF h5k<>0 and h5oldk=0 then !CALL Plane3Redraw CALL h3DrawSlider(h3c) CALL h4DrawSlider(h4d) CALL t2Init LET w3Flag= 1 !CALL Plane1Update CALL UpdateGraphs CALL t5Set END IF BOX SHOW Env$ at w2Lft,w2Bas CALL Plane2Refresh CALL Plane2Envelope CALL Plane2Wave !CALL UpdateGraphs END SUB ! --- Read Slider 6 --- SUB h6MouseClick(mx,my) CALL MouseUp(mx,my,ms) LET h6j= h6Fncx(clamp(mx,h6wLft,h6wRgt)) LET h6j= roundn(h6j,h6LTik) IF h6j<>h6oldj then CALL h6Action LET h6oldj= h6j END IF END SUB SUB h6MouseDrag(mx,my) CALL Plane2Refresh DO GET MOUSE: mx,my,ms LET h6j= round(h6Fncx(clamp(mx,h6wLft,h6wRgt)),2) IF h6j<>h6oldj then CALL h6Action LET h6oldj= h6j END IF LOOP until ms=3 END SUB SUB h6Action CALL h6mark(h6j) IF h6j=0 and h6oldj<>0 then !CALL Plane1Clear BOX CLEAR h1wLft-20,h1wRgt+50,h2wBas+10,t1Top-10 LET w1Flag= 0 !CALL Plane1Update CALL UpdateGraphs CALL t5Set ELSE IF h6j<>0 and h6oldj=0 then !CALL Plane1Redraw CALL h1DrawSlider(h1a) CALL h2DrawSlider(h2b) CALL t1Init LET w1Flag= 1 !CALL Plane1Update CALL UpdateGraphs CALL t5Set END IF BOX SHOW Env$ at w2Lft,w2Bas CALL Plane2Refresh CALL Plane2Envelope CALL Plane2Wave !CALL UpdateGraphs END SUB ! --- Graph drawing methods --- SUB UpdateGraphs !IF w1Flag=1 then CALL Plane1Refresh CALL Plane1Update !END IF ! IF w3Flag=1 then ! CALL Plane3Refresh ! CALL Plane3Update ! END IF !IF w2Flag=1 then CALL Plane2Refresh CALL Plane2Envelope BOX KEEP w2Lft,w2Rgt,w2Bas,w2Top in Env$ CALL Plane2Wave !END IF END SUB SUB Plane1Update ! --- Points Plane 1 CALL Plane1Refresh IF w1Flag=1 then LET wx = w1wndx( h1a) LET wy1= w1wndy( h2b) LET wy2= w1wndy(-h2b) DRAW diamond7(abclr) with shift(wx,wy1) DRAW diamond7(abclr) with shift(wx,wy2) END IF IF w3Flag=1 then LET wx = w1wndx( h3c) LET wy1= w1wndy( h4d) LET wy2= w1wndy(-h4d) DRAW diamond7(cdclr) with shift(wx,wy1) DRAW diamond7(cdclr) with shift(wx,wy2) END IF END SUB SUB Plane3Update ! --- Points Plane 3 CALL Plane1Update ! LOCAL wx,wy1,wy2 ! SET COLOR cclr ! LET wx = w3wndx( h3c) ! LET wy1= w3wndy( h4d) ! LET wy2= w3wndy(-h4d) END SUB SUB Plane2Wave ! --- Wave Plane2 LOCAL x,wx,oldwx,y,wy,oldwy,wx1,wy1,wx2,wy2 LET oldwx= w2Lft LET x = w2fncx(w2Lft) LET y = wave(x) LET oldwy= w2wndy(y) FOR wx= w2Lft to w2Rgt step 2 LET x = w2fncx(wx) * pi LET y = wave(x) LET wy= w2wndy(y) SET COLOR bclr IF wy>w2Top and wyw2Top and oldwyw2Top and wyw2Top and oldwyw2Bas and oldwyw2Top) then IF wy>oldwy then PLOT oldwx,w2top+1; wx,w2bas-1 ELSE PLOT wx,w2top+1; oldwx,w2bas-1 END IF END IF LET oldwx= wx LET oldwy= wy NEXT wx PLOT END SUB ! DEF envl(x)= h6j*e^(h1a*x) + h5k*e^(h3c*x) ! DEF wave(x)= h6j*e^(h1a*x)*cos(h2b*x) + h5k*e^(h3c*x)*cos(h4d*x) SUB Plane2Envelope ! --- Envelope Plane 2 LOCAL x,wx,y,wy SET COLOR aclr SET CURSOR 1,1 ! LET h6j= round(h6j,5) ! PRINT h1a,h3c,h6j FOR wx= w2Lft to w2Rgt step 4 LET x= w2fncx(wx) * pi IF (h5k<>0 and h6j<>0 and h1a>h3c) or h5k=0 then LET y= h6j*e^(h1a*x) ELSE IF (h6j<>0 and h5k<>0 and h1a0 and h6j<>0 then IF h6j>=h5k then LET y= h6j*e^(h1a*x) ELSE LET y= h5k*e^(h3c*x) END IF ELSE LET y= envl(x) END IF ! DEF envl(x)= h6j*e^(h1a*x) + h5k*e^(h3c*x) IF y>w2fBas and y