!When you click in the window while the amination is going !you enter a state which computes this partial sum incorrectly. ! !In the case of the spike there is a clipping problem with !this rule when all the terms are turned on and t = 0. !Heat: Looking good. ! !- the chosen initial condition button should be highlighed somehow. ! !- the vertical axis should be u , not y. ! !- the "spike" should be marked with a vertical line segment extending from u = 0 to u = 1.5, !and ending with an arrow head. The segment might be drawn bolder than normal. ! !- Let's say "ramp" rather than "taper." ! !- Let's draw the white curve AFTER the terms have been drawn in, not before. ! !- Can we have the animation start from the current value of t, rather than always !from t = 0? ! !all of the above are easy ! ! !- We were going to speed this up by using symmetry when possible (ie in the spike and uniform cases). ! !I tried this but giving the nature of rolling averages it seems necessary to compute all the !way from left to right to get the right curve. ! !- I have to admit that the wire doesn't add very much for me. This color representation !is effective in two dimensions but somehow you don't get enough color, or you only get to !compare this color to the color at left and at right, and the effect isn't very striking. ! !Easy to turn off the wire. ! !- The initial animation is fun but it takes a while. Will it annoy the student? !I think at least we should abort it at a mouse button. ! !I am calculating and storing the values here for fast animation - we could put a !message up "computing curves" or something - it seemed a bit boring to wait in !front of a blank screen when you could get introduced to the three states of the wire. !I could put the data in a file, but it takes time to move it from a file into arrays too, !and files tend to separated from the programs they support. ! !- A vertical rule on the screen would encourage students to think of the white curve !as the sum of the colored ones. ! !I was thinking about a color coded vertical so you can see the pile of values. I think !I did that in the old rolling average sketch. ! !- There is some weird bug having to do with the finite element thing. You have changed !the pixel size for time from .02 to .01, but you don't recompute the finite element !solution every time. One effect of this is: If you start from uniform and step !forward two times, the two white curves are identical (and they shouldn't be). If you !now step back, so t = .01, you get a different curve than was shown at first for t = .01. !If you now switch to "spike," you get a broadened spike, exiting the window (and hence !not erased when you go on). ! !-Presumably this can be fixed just by computing 25 time steps instead of 50, to get a time !change of .01 instead of .02. You'll have to clip the smeared spike. ! !This may be just a rounding issue - I was rounding to the nearest .02 but might not !be doing it for the click steps. ! !-In fact, it might be a good idea to allow a "zoom in" on the very small values of t; !but perhaps that can wait. ! !-I'm not sure how you step back in the finite element method! It might be better to eliminate !the left arrow key. ! !With all the curves precomputed, forward or backward is just a matter of looking !up the list of values for a particular t. !Fix k = 1. ! !s = 4/(pi^2). ! !Delta t = .0004 ! !so that there are 50 steps of the finite element computation for !each time step of .02 in the Fourier tool. !! File: Heat Equation !! December 8, 2001 Hubert Hohn for Prof. Haynes Miller 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$ = "Fourier Series: heat dissipation from the ends of an insulated wire" SUB ThisProgram CALL heat 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+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,b1,b0 ! red SET COLOR MIX( 8) b5,b5,b0 ! yellow SET COLOR MIX( 9) b1,b5,b0 ! green SET COLOR MIX(10) b0,b5,b3 ! cyan SET COLOR MIX(11) b0,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 FOR i= 20 to 40 LET f = (i-20)/20 LET f0= 1-f SET COLOR MIX(i) f,0,f0 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 ! --- 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(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= abs(round(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(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(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= abs(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 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= 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 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= 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= 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 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+45,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 h2Slider 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 h2SliderVariables 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 h3Slider 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 h3SliderVariables 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-28,Bas-3,name$,clr) END SUB SUB PlotSliderValue(Rgt,Bas,value$,clr) CALL SetTextFont(1,12,"bold") BOX CLEAR Rgt+30,Rgt+80,Bas,Bas-13 CALL PlotTextLJ(Rgt+33,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 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 AlignChar(chx,ty,clr,txt$,char$) ! LET p = pos(txt$,char$) ! LET lft$= txt$(1:p) ! CALL StringWidth(lft$,sl) ! SET COLOR 10 ! PLOT TEXT, AT chx-sl+1,ty+1: txt$ ! SET COLOR clr ! PLOT TEXT, AT chx-sl,ty: txt$ ! END SUB ! SUB AlignDot(ec,ex,ey,eq$) ! LET e = pos(eq$,".") ! LET lft$= eq$(1:e) ! CALL StringWidth(lft$,sl) ! SET COLOR ec ! PLOT TEXT, AT ex-sl,ey: eq$ ! END SUB ! ! SUB AlignEqualMono(ec,ex,ey,eq$) ! LET e = pos(eq$,"=") ! LET lft$= eq$(1:e) ! LET rgt$= eq$(e+1:len(eq$)) ! CALL StringWidth(lft$,sl) ! SET COLOR 10 ! PLOT TEXT, AT ex-sl+1,ey+1: lft$ ! SET COLOR ec ! PLOT TEXT, AT ex-sl,ey: lft$ ! CALL MonoSpace(9,0,ec,ex,ey,rgt$) ! 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 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 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 !BOX KEEP iLft,iRgt,iBas,iTop in help$ SET COLOR white BOX AREA iLft,iRgt,iBas,iTop CALL CopyRight(iLft+20,ibas-3,2) !CALL Rim(iLft,iRgt,iBas,iTop,litgry) 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,size,"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 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 ticks3= 0 LET t1 = time DO LET ticks3= ticks3+1 LOOP until time-t1>=1 LET unitTicks= (ticks1+ticks2+ticks3)/3 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 ! --- End of Library --- ! *** heat Program *** SUB heat 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 Fourier Series and heat dissipation" ! ---------- Utility functions --- DEF clamp(n,lo,hi)= min(max(n,lo),hi) DEF roundn(n,step)= round(n/step)*step DEF e= exp(1) ! ------- functions ----- ! ------ square wave ---- DEF fn(t,x,term) LET trmval= term*2 - 1 LET fn= (4/pi) * (1/trmval) * e^( -trmval*trmval*k*t) * sin(trmval*x) END DEF DEF f(t,x,n) LET out= 0 LET c = k*t FOR term= 1 to n LET trmval= term*2 - 1 LET out = out + (1/trmval) * e^( -trmval*trmval*c) * sin(trmval*x) NEXT term LET f= 4/pi * out END DEF LET f$= "{4/pi} {e^[-kt]sin x + (e^[-9kt]/3) sin 3x + (e^[-25kt]/5) sin 5x + (e^[-49kt]/7) sin 7x + (e^[-81kt]/9) sin 9x}" ! ---------- 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 DECLARE PUBLIC w1xPiFlag, w1wHgt, w1fHgt, w1wWid, w1fWid LET w1Flag= 1 LET w1xPiFlag= 1 LET w1Lft = workLft+90 ! pixel bounds LET w1Rgt = w1Lft+600 LET w1Top = workTop+35 LET w1Bas = w1Top+280 LET w1fLft= 0 ! function bounds LET w1fRgt= 1 LET w1fTop= 1.5 LET w1fBas= -0.5 LET w1mlt = pi LET w1Xax$= "x" ! axis labels LET w1Yax$= "u" LET w1xGridstep= 0 ! horizontal grid intervals LET w1yGridstep= 0 ! vertical grid intervals LET w1SxTik = 0 ! horizontal axis Tik marks LET w1LxTik = 0.25 LET w1xLabel= 0.5 LET w1Firstx= w1fLft LET w1SyTik = 0.1 ! vertical axis Tik marks LET w1LyTik = 0.5 LET w1yLabel= 0.5 LET w1Firsty= w1fBas ! --- Plane 1 methods --- DECLARE DEF w1fncx,w1fncy,w1wndx,w1wndy ! window/function transforms CALL Plane1Variables DIM lftlist(0:1,0:1),rgtlist(0:1,0:1),sumlist(0:1,0:1) MAT redim lftlist(0:w1wWid,0:1), rgtlist(0:w1wWid,0:1), sumlist(0:w1wWid,0:1) 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$ 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, w2x0, w2y0 DECLARE PUBLIC w2Firstx, w2SxTik, w2LxTik, w2xLabel, w2xGridstep DECLARE PUBLIC w2Firsty, w2SyTik, w2LyTik, w2yLabel, w2yGridstep DECLARE PUBLIC w2xPiFlag, w2wHgt, w2fHgt LET w2Flag = 1 LET w2xPiFlag= 1 LET w2Lft= w1Lft ! pixel bounds LET w2Rgt= w1Rgt LET w2Top= w1Bas + 30 LET w2Bas= w2Top + 10 LET w2fLft= 0 ! function bounds LET w2fRgt= 1 LET w2fTop= 1.25 LET w2fBas= 0 LET w2mlt = pi LET w2Xax$= "x" ! axis labels LET w2Yax$= "y" LET w2xGridstep= 0 ! horizontal grid intervals LET w2yGridstep= 0 ! vertical grid intervals LET w2SxTik = 0 ! horizontal axis Tik marks LET w2LxTik = .25 LET w2xLabel= .5 LET w2Firstx= w2fLft LET w2SyTik = 0 ! vertical axis Tik marks LET w2LyTik = .5 LET w2yLabel= .5 LET w2Firsty= w2fBas ! --- Plane 2 methods --- DECLARE DEF w2fncx,w2fncy,w2wndx,w2wndy ! window/function transforms CALL Plane2Variables SUB Plane2Clear BOX CLEAR w2Lft-20,w2Rgt+20,workBas-5,w2Top-20 END SUB SUB Plane2Init SET COLOR 4 BOX LINES w2Lft,w2Rgt,w2Bas+2,w2Bas+1 BOX LINES w2Lft,w2Rgt,w2Top-1,w2Top-2 END SUB SUB Plane2Refresh BOX SHOW w2gridLayer$ at w2Lft-5,w2Bas+5 END SUB ! ----------- horizontal sliders ------------ LET slideWidth= 120 ! --- 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 = white LET h1name$= "t" LET h1form$= "--%.##" LET h1m = 1 LET h1a = 0 LET h1olda = -999 LET h1axis = w1Bas + 120 LET h1wLft = w1Lft LET h1wRgt = h1wLft + 100 LET h1fLft = 0 LET h1fRgt = 1 LET h1mlt = 1 LET h1STik = 0.1 ! short tick marks LET h1LTik = 0.5 ! long tick marks LET h1Label= 0.5 ! labels LET h1First= 0 ! first tick mark !LET h1xPix = (h1fRgt-h1fLft)/(h1wRgt-h1wLft) LET h1xPix = .02 CALL h1SliderVariables SUB h1Init(t) CALL h1DrawSlider(t) LET h1anmLft= h1wLft - 24 LET h1anmRgt= h1anmLft+18 LET h1anmTop= h1axis LET h1anmBas= h1axis+12 CALL DrawButton(h1anmLft,h1anmRgt,h1anmBas,h1anmTop,3,">>") ! LET midbtn= int((h1anmLft+h1anmRgt)/2) ! SET COLOR white ! DRAW StepArrow with shift(midbtn+2,h1anmBas-6) ! DRAW StepArrow with shift(midbtn-2,h1anmBas-6) LET h1lstLft= h1wRgt+ 6 LET h1lstRgt= h1lstLft+10 LET h1lstTop= h1anmtop LET h1lstBas= h1anmbas CALL DrawButton(h1lstLft,h1lstRgt,h1lstBas,h1lstTop,3,"<") LET h1rstLft= h1lstRgt+ 2 LET h1rstRgt= h1rstLft+10 LET h1rstTop= h1anmtop LET h1rstBas= h1anmbas CALL DrawButton(h1rstLft,h1rstRgt,h1rstBas,h1rstTop,3,">") END SUB PICTURE StepArrow PLOT -2,-4; 2,0; -2,4; -2,-4 END PICTURE ! --- horizontal 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 = white LET h2name$= "k" LET h2form$= "--%.##" LET h2m = 1 LET h2a = 0 LET h2olda = -999 LET h2axis = h1axis + 45 LET h2wLft = h1wLft LET h2wRgt = h1wRgt LET h2fLft = 0 LET h2fRgt = 4 LET h2mlt = 1 LET h2STik = 0.1 ! short tick marks LET h2LTik = 0.5 ! long tick marks LET h2Label= 0.5 ! labels LET h2First= 0 ! first tick mark LET h2xPix = (h2fRgt-h2fLft)/(h2wRgt-h2wLft) CALL h2SliderVariables SUB h2Init(k) CALL h2DrawSlider(k) LET h2anmLft= h2wLft - 24 LET h2anmRgt= h2anmLft+18 LET h2anmTop= h2axis LET h2anmBas= h2axis+12 CALL DrawButton(h2anmLft,h2anmRgt,h2anmBas,h2anmTop,3,">>") LET h2lstLft= h2wRgt+ 6 LET h2lstRgt= h2lstLft+10 LET h2lstTop= h2anmtop LET h2lstBas= h2anmbas CALL DrawButton(h2lstLft,h2lstRgt,h2lstBas,h2lstTop,3,"<") LET h2rstLft= h2lstRgt+ 2 LET h2rstRgt= h2rstLft+10 LET h2rstTop= h2anmtop LET h2rstBas= h2anmbas CALL DrawButton(h2rstLft,h2rstRgt,h2rstBas,h2rstTop,3,">") END SUB ! --- horizontal slider 1 --- 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 = white LET h3name$= "terms" LET h3form$= "--%.##" LET h3m = 1 LET h3a = 0 LET h3olda = -999 LET h3axis = h2axis + 45 LET h3wLft = h1wLft LET h3wRgt = h3wLft + 200 LET h3fLft = 0 LET h3fRgt = 20 LET h3mlt = 1 LET h3STik = 0 ! short tick marks LET h3LTik = 1 ! long tick marks LET h3Label= 5 ! labels LET h3First= 0 ! first tick mark CALL h3SliderVariables ! ---------- Text Output Rects ---------- LET tlnspc = 22 ! --- text rectangle 1 --- LET t1BasLn = w2Bas + 40 LET t1Lft = w2Lft - 2 LET t1Rgt = workRgt LET t1Bas = t1BasLn + 5 LET t1Top = t1BasLn - 15 LET t1Label$= "u(x,t) = " SUB t1Label CALL SetTextFont(1,12,"bold") CALL SuperSubScriptRJ(t1Lft,t1BasLn,t1Label$,white) END SUB SUB t1Set CALL SetTextFont(1,12,"bold") CALL t1Clear CALL SuperSubScriptLJ(w1Lft,w2Bas+40,f$,white) 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 + tlnspc LET t2Lft = t1Lft LET t2Rgt = t1Rgt LET t2Bas = t2BasLn + 5 LET t2Top = t2BasLn - 15 LET t2Label$= "phi = " SUB t2Label CALL SetTextFont(1,12,"bold") CALL SuperSubScriptRJ(t2Lft,t2BasLn,t2Label$,white) END SUB SUB t2Set CALL SetTextFont(1,12,"bold") CALL t2Clear LET t$ = "atn(w/k)" 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 ! --- text rectangle 3 --- LET t3BasLn = t2BasLn + tlnspc LET t3Lft = t1Lft LET t3Rgt = t1Rgt LET t3Bas = t3BasLn + 5 LET t3Top = t3BasLn - 15 LET t3Label$= "A = " SUB t3Label CALL SetTextFont(1,12,"bold") CALL SuperSubScriptRJ(t3Lft,t3BasLn,t3Label$,white) END SUB SUB t3Set CALL SetTextFont(1,12,"bold") CALL t3Clear LET t$= "cos(phi)" CALL PlotTextLJ(t3Lft,t3BasLn,t$,white) END SUB SUB t3Clear BOX CLEAR t3Lft-2,t3Rgt,t3Bas,t3Top END SUB SUB t3Init CALL t3Label CALL t3Set END SUB ! --- Check boxes --- DIM chkbox$(1:5),chkbox(1:5) DATA "term 1","term 2","term 3","term 4","term 5" MAT READ chkbox$ MAT chkbox= Zer LET chkcnt= 5 LET chksiz= 10 LET chkstp= 16 LET chklft= w1Rgt - 65 LET chkrgt= chklft + chksiz LET chktop= w1Bas+105 LET chkbas= chktop + chkcnt*chkstp ! --- Check boxes --- DIM modelbox$(1:3),modelbox(1:3) DATA "uniform","spike","ramp" MAT READ modelbox$ MAT modelbox= Zer LET modelcnt= 3 LET modelsiz= 10 LET modelstp= 16 LET modellft= w1Midx LET modelrgt= modellft + modelsiz LET modeltop= chktop LET modelbas= modeltop + modelcnt*modelstp ! --- Finite and Fourier --- ! ------ Fourier Series Functions ---- DEF fsqrterm(t,x,term) LET termval = term*2 - 1 ! odd series LET c = k*t LET fsqrterm= (4/pi) * e^( -termval*termval*c) * sin(termval*x) / termval END DEF DEF fsqr(t,x,terms) LET out= 0 LET c = k*t FOR term= 1 to terms LET termval= term*2 - 1 LET out = out + (1/termval) * e^(-termval*termval*c) * sin(termval*x) NEXT term LET fsqr= 4/pi * out END DEF DEF fdltterm(t,x,term) LET termval= term*2 - 1 LET c = k*t LET odd = mod(term,2) IF odd=1 then LET fdltterm= + (1/pi) * e^(-termval*termval*c) * sin(termval*x) ELSE LET fdltterm= - (1/pi) * e^(-termval*termval*c) * sin(termval*x) END IF END DEF DEF fdlt(t,x,terms) LET out= 0 LET c = k*t FOR term= 1 to terms LET termval= term*2 - 1 LET odd = mod(term,2) IF odd=1 then LET out= out + e^(-termval*termval*c) * sin(termval*x) ELSE LET out= out - e^(-termval*termval*c) * sin(termval*x) END IF NEXT term LET fdlt= 1/pi * out END DEF DEF frmpterm(t,x,term) LET c = k*t LET odd= mod(term,2) IF odd=1 then LET frmpterm= + (2/pi) * e^(-term*term*c) * sin(term*x) / term ELSE LET frmpterm= - (2/pi) * e^(-term*term*c) * sin(term*x) / term END IF END DEF DEF frmp(t,x,terms) LET out= 0 LET c = k*t FOR term= 1 to terms LET odd= mod(term,2) IF odd=1 then LET out= out + e^(-term*term*c) * sin(term*x) / term ELSE LET out= out - e^(-term*term*c) * sin(term*x) / term END IF NEXT term LET frmp= 2/pi * out END DEF ! ---- Finite Elements Arrays ---- DIM ulist(0:100),oldu(0:100) DIM sqrArray(0:200,0:100) DIM dltArray(0:200,0:100) DIM rmpArray(0:200,0:100) LET xsteps= 100 LET maxu = xsteps ! maximum array element LET k= 1 MAT redim ulist(0:maxu) MAT redim oldu(0:maxu) MAT redim sqrArray(0:200,0:maxu) MAT redim dltArray(0:200,0:maxu) MAT redim rmpArray(0:200,0:maxu) ! ---- Initialize Finite Elements Models ---- SUB LoadFiniteArrays FOR model= 0 to 2 SELECT CASE model CASE 0 ! square MAT ulist = 1 LET ulist(0) = 0 LET ulist(maxu)= 0 CASE 1 ! delta MAT ulist = ZER LET n = maxu/2 LET ulist(n)= n/pi CASE 2 ! ramp FOR i= 0 to maxu LET ulist(i)= i/maxu NEXT i LET ulist(maxu)= 0 END SELECT CALL CopyList(0,ulist) SET COLOR yellow CALL FiniteCurve(model,0) ! 50 averages here is a .02 step in time LET tsteps= 50 ! # of times array is processed LET s = 4/pi^2 FOR n= 1 to curves FOR i= 1 to tsteps ! update ulist array MAT oldu= ulist FOR xn= 1 to maxu-1 LET ulist(xn)= oldu(xn) + s*(oldu(xn-1) + oldu(xn+1) - 2*oldu(xn)) NEXT xn NEXT i CALL CopyList(n,ulist) CALL FiniteCurve(model,n) NEXT n CALL Plane1Refresh NEXT model END SUB SUB CopyList(t,list()) SELECT CASE model CASE 0 FOR i= 0 to maxu LET sqrArray(t,i)= list(i) NEXT i CASE 1 FOR i= 0 to maxu LET dltArray(t,i)= list(i) NEXT i CASE 2 FOR i= 0 to maxu LET rmpArray(t,i)= list(i) NEXT i END SELECT END SUB SUB StoreList(t,list(),array(,)) PRINT #1: "data " FOR i= 0 to maxu LET array(t,i)= list(i) PRINT #1: list(i); IF iw1Lft and mxw1Top and mymx then CALL RollOverSum(mx,sum) BOX SHOW temp$ at w1Lft,w1Bas SET COLOR 5 LET wy= max(w1Top,w1Wndy(sum)) PLOT mx,w1y0; mx,wy LET cleanup= 1 LET oldmx = mx END IF ELSE IF cleanup=1 then BOX SHOW temp$ at w1Lft,w1Bas LET cleanup= 0 END IF LOOP until ms=2 IF mx>=h1wLft-5 and mx<=h1wRgt+5 and my>h1wTop and myh1anmlft and mxh1anmtop and myh1lstlft and mxh1lsttop and myh1rstlft and mxh1rsttop and my=h2wLft-5 and mx<=h2wRgt+5 and my>h2wTop and myh2anmlft and mxh2anmtop and myh2lstlft and mxh2lsttop and myh2rstlft and mxh2rsttop and mychklft and mxchktop and mymodellft and mxmodeltop and mywirelft and mxwiretop and mysqLft and mxsqTop and myinfLft and mxinfTop and myqLft and mxqTop and my0 then LET sum= 0 SELECT CASE model CASE 0 LET x= w1Fncx(mx)*pi FOR term= 1 to 5 IF chkbox(term)=1 then LET y = fsqrterm(t,x,term) LET sum= sum+y END IF NEXT term CASE 1 LET x= w1Fncx(mx)*pi FOR term= 1 to 5 IF chkbox(term)=1 then LET y = fdltterm(t,x,term) LET sum= sum+y END IF NEXT term CASE 2 LET x= w1Fncx(mx)*pi FOR term= 1 to 5 IF chkbox(term)=1 then LET y = frmpterm(t,x,term) LET sum= sum+y END IF NEXT term END SELECT END IF END SUB SUB h1SliderAction IF t<>oldt then LET oldt= t LET t= min(max(t,0),1) CALL h1Mark(t) CALL Plane1Refresh LET n= t*50 CALL FiniteCurve(model,n) CALL GraphUofX(t) CALL GraphUofXWire(model,t) END IF END SUB SUB h2SliderAction IF k<>oldk then LET oldk= k CALL h2Mark(k) CALL GraphUofX(t) CALL GraphUofXWire(model,t) END IF END SUB ! --- Tool specific drawing and typesetting routines SUB GraphUofX(t) ! 1 term IF termflag>0 then SELECT CASE model CASE 0 FOR pix= 0 to w1wWid/2 LET wx1= w1Midx+pix LET wx2= w1Midx-pix !FOR wx= w1Lft to w1Rgt step 1 LET x = w1Fncx(wx1)*pi FOR term = 1 to 5 IF chkbox(term)=1 then SET COLOR red+(term-1) !LET y= fn(t,x,term) LET y= fsqrterm(t,x,term) PLOT wx1,w1Wndy(y) PLOT wx2,w1Wndy(y) END IF NEXT term !NEXT wx NEXT pix CASE 1 FOR pix= 0 to w1wWid/2 LET wx1= w1Midx+pix LET wx2= w1Midx-pix LET x = w1Fncx(wx1)*pi FOR term = 1 to 5 IF chkbox(term)=1 then SET COLOR red+(term-1) !LET y= fn(t,x,term) LET y= fdltterm(t,x,term) PLOT wx1,w1Wndy(y) PLOT wx2,w1Wndy(y) END IF NEXT term !NEXT wx NEXT pix CASE 2 FOR pix= 0 to w1wWid/2 LET wx1= w1Midx+pix LET wx2= w1Midx-pix LET x1 = w1Fncx(wx1)*pi LET x2 = w1Fncx(wx2)*pi FOR term = 1 to 5 IF chkbox(term)=1 then SET COLOR red+(term-1) !LET y= fn(t,x,term) LET y= frmpterm(t,x1,term) PLOT wx1,w1Wndy(y) LET y= frmpterm(t,x2,term) PLOT wx2,w1Wndy(y) END IF NEXT term ! NEXT wx NEXT pix END SELECT END IF END SUB ! SUB GraphUofXWire(t) ! 1 term ! IF wireFlag=1 then ! FOR wx= w1Lft to w1Rgt step ! LET x = w1Fncx(wx)*pi ! LET u = f(t,x,terms) ! LET clr= 20 + round((u - w1fBas)/w1fHgt * 20) ! SET COLOR clr ! PLOT wx,w2Top; wx,w2Bas ! NEXT wx ! END IF ! END SUB SUB GraphUofXWire(model,t) ! 1 term LET n= t*50 IF wireFlag=1 then LET c= pi/maxu ! now plot ulist array FOR xn= 0 to maxu LET x = c*xn SELECT CASE model CASE 0 LET u= sqrArray(n,xn) CASE 1 LET u= dltArray(n,xn) CASE 2 LET u= rmpArray(n,xn) END SELECT LET clr= 20 + round((u - w1fBas)/w1fHgt * 20) LET clr= min(clr,40) SET COLOR clr LET wx = w1Wndx(x/pi) LET lft= wx-3 LET rgt= wx+3 IF xn=0 then LET lft= w1Lft IF xn=maxu then LET rgt= w1Rgt BOX AREA lft,rgt,w2Bas,w2Top NEXT xn END IF END SUB ! LET c= pi/maxu ! now plot ulist array ! SELECT CASE model ! CASE 0 ! square ! FOR xn= 0 to maxu ! LET x= c*xn ! LET y= sqrArray(n,xn) ! PLOT w1Wndx(x/pi),w1Wndy(y); ! NEXT xn SUB Animate ! 10 terms FOR tAnim= t to 1 step .02 LET tAnim= round(tAnim,2) CALL h1Mark(tAnim) CALL Plane1Refresh CALL FiniteCurve(model,tAnim*50) CALL GraphUofX(tAnim) !CALL GraphUofXWire(model,tAnim) CALL Delay(1/64) GET MOUSE: mx,my,ms IF ms=3 then LET t= min(tAnim,1) LET t= round(t,2) EXIT SUB END IF NEXT tAnim END SUB ! SUB AnimateK ! 10 terms ! LET wait= 1/30 ! FOR k= 0 to 4 step h2STik ! LET k= round(k,2) ! CALL h2Mark(k) ! CALL GraphUofX(t) ! CALL GraphUofXWire(model,t) ! PAUSE wait ! GET MOUSE: mx,my,ms ! IF ms=3 then EXIT SUB ! NEXT k ! END SUB ! --- graphing and interface parts --- SUB Buttons ! --- Button parameters and methods LET bhgt = 17 LET wirelft= w2Lft-57 LET wirergt= wirelft+35 LET wiretop= w2Top-2 LET wirebas= w2Bas+3 CALL SetTextFont(1,9,"bold") CALL DrawButton(wirelft,wirergt,wirebas,wiretop,5,"wire") LET sqlft= w1Midx LET sqrgt= sqlft+50 LET sqtop= chktop LET sqbas= sqtop+bhgt-2 CALL DrawButton(sqlft,sqrgt,sqbas,sqtop,5,"uniform") LET dllft= w1Midx LET dlrgt= dllft+50 LET dltop= sqbas+2 LET dlbas= dltop+bhgt-2 CALL DrawButton(dllft,dlrgt,dlbas,dltop,5,"spike") LET rmlft= w1Midx LET rmrgt= rmlft+50 LET rmtop= dlbas+2 LET rmbas= rmtop+bhgt-2 CALL DrawButton(rmlft,rmrgt,rmbas,rmtop,5,"taper") END SUB ! --- SUB SetCheckBox(n,state) LET top= chktop + (n-1)*chkstp LET bas= top+chksiz IF state=1 then SELECT CASE n CASE 1 SET COLOR red CASE 2 SET COLOR yellow CASE 3 SET COLOR green CASE 4 SET COLOR cyan CASE 5 SET COLOR blue END SELECT ELSE SET COLOR 1 END IF BOX AREA chklft+2,chkrgt-2,bas-2,top+2 LET termflag= 0 FOR i= 1 to 5 IF chkbox(i)=1 then LET termflag=1 NEXT i END SUB SUB DrawCheckBoxes CALL SetTextFont(1,12,"bold") CALL PlotTextCJ(chklft+35,chktop-27,"Fourier",white) CALL PlotTextCJ(chklft+35,chktop-12,"Terms",white) FOR btn= 1 to chkcnt LET top = chktop + (btn-1)*chkstp LET bas = top + chksiz CALL SetTextFont(1,12,"bold") LET txt$= chkbox$(btn) CALL PlotTextLJ(chkrgt+5,bas-1,txt$,white) SET COLOR 4 BOX AREA chklft-1,chklft+chksiz+1,bas+1,top-1 SET COLOR 1 BOX AREA chklft,chklft+chksiz,bas,top NEXT btn IF checktitle$<>"" then CALL SetTextFont(1,12,"bold") CALL PlotTextLJ(chklft+17,chktop-13,checktitle$,black) SET COLOR black PLOT chklft,chktop-8; chklft+txtwid+18,chktop-8 END IF END SUB ! --- SUB SetModelBox(n,state) FOR i= 1 to 3 LET top= modeltop + (i-1)*modelstp LET bas= top+modelsiz SET COLOR black BOX AREA modellft+2,modelrgt-2,bas-2,top+2 NEXT i LET top= modeltop + (n-1)*modelstp LET bas= top+modelsiz !IF state=1 then SET COLOR white ! ELSE ! SET COLOR black ! END IF BOX AREA modellft+2,modelrgt-2,bas-2,top+2 ! LET termflag= 0 ! FOR i= 1 to 3 ! IF modelbox(i)=1 then LET termflag=1 ! NEXT i END SUB SUB DrawModelBoxes CALL SetTextFont(1,12,"bold") CALL PlotTextCJ(modellft+35,modeltop-27,"Heat",white) CALL PlotTextCJ(modellft+35,modeltop-12,"Distribution",white) FOR btn= 1 to modelcnt LET top = modeltop + (btn-1)*modelstp LET bas = top + modelsiz LET txt$= modelbox$(btn) CALL PlotTextLJ(modelrgt+5,bas-1,txt$,white) SET COLOR 4 BOX AREA modellft-1,modellft+modelsiz+1,bas+1,top-1 SET COLOR 1 BOX AREA modellft,modellft+modelsiz,bas,top NEXT btn IF checktitle$<>"" then CALL SetTextFont(1,12,"bold") CALL PlotTextLJ(modellft+17,modeltop-13,checktitle$,black) SET COLOR black PLOT modellft,modeltop-8; modellft+txtwid+18,modeltop-8 END IF END SUB END SUB ! --- end of tool code -------------------