1000 'DIANEQNS.TXT=TEST EQUATIONS FOR MERGING WITH DIANALOG.RENUMBER THE LINES TO SUIT BY REMOVING LEFT DIGIT/S OF EACH LINE NUMBER. 1005 =========================================================================== 1010 X(1,2)=2*X(2,0)-5*X(1,0) 1011 X(2,2)=2*X(1,0)-2*X(2,0) 1012 X(3,5)=X(3,1)+1 1900 '2 2nd ORDER SPRINGS+ A 5th order.SEE MAL BOOK TO RUN versus TIME.DELETE EQN ON LINE 1012 IN ORDER TO GRAPH THE 2 SPRINGS EQNS REMAINING.ANSWER QUESTIONS WITH:2,2,0,.1,2,5,0,2,0,0,0,CONT,40,4,0,4,20,1 2000 =========================================================================== 2010 V1=SIN(100*3.14159*T):IF V1<0 THEN V1=0 2011 X(1,2)=500000!*V1-102*X(1,1)-2200*X(1,0) 2170 FOR J=1 TO Q+1:PRINT"X";J;" SCALE MULTIPLIER M(";J;")=";:INPUT M(J) 2590 X(2,0)=V1:FOR J=0 TO Q+1:V(J)=INT(V+M(J)*X(J,0)+Y(J)+.5):NEXT J:'CAN PLOT OTHER VARIABLES 2640 FOR J=0 TO Q+1:IF J=0 AND N0=10 THEN N0=0:GOTO 690 2745 'CAN PUT TIME T OR DESIRED DIFFERENTIALS,VARIABLES,FUNCTIONS OF THEM OR AXES TRANSFORMS,INSTEAD OF X(2,0).CAN CHANGE SCALE FACTORS DURING RUN WITH BREAK,CONT. 2750 VY=362-K6*16-(YO-YM*T)*W1*.875/W0:'THIS LINE SET UP FOR A PLOT versus TIME 2900 'TO GRAPH ELECTRICAL FILTER DIFF.EQUATION IN DIANALOG,ANSWER QUESTIONS WITH FOLLOWING NUMBERS:1,2,0,.0005,0,0,0,2,CONT,0,.5,0,400,50,4:ALSO SEE MAL BOOK TO RUN TIME PLOT. 3000 =========================================================================== 3010 X(1,2)=-X(1,0) 3011 X(2,2)=-X(2,0) 3900 'TO GRAPH CIRCLE DIFF.EQUATION IN DIANALOG,ANSWER QUESTIONS WITH FOLLOWING NUMBERS:2,2,0,.1,2,0,1,2,1,0,0,2,CONT,40,25,0,25,25,1 4000 =========================================================================== 4010 X(1,2)=1+X(2,1) 4011 X(2,2)=-X(1,1) 4900 'TO GRAPH ELECTRON DIFF.EQUATION IN DIANALOG,ANSWER QUESTIONS WITH FOLLOWING NUMBERS:2,2,0,.1,2,0,4,2,0,0,0,2,CONT,40,4,0,1,35,1 5000 =========================================================================== 5010 X(1,2)=X(1,0)*X(2,1)^2-2/X(1,0)^2 5011 X(2,2)=-2*X(1,1)*X(2,1)/X(1,0) 5740 XA=X(1,0)*COS(X(2,0)):VX=(V+.5+XO+XM*XA)*W1/W0:'ORBIT AXES TRANSFORMS 5750 YA=X(1,0)*SIN(X(2,0)):VY=362-(YO+YM*YA)*W1*.875/W0 5900 'TO GRAPH ORBIT DIFF.EQUATION IN DIANALOG,ANSWER QUESTIONS WITH FOLLOWING NUMBERS:2,2,0,.2,2,3(4.5 FOR ESCAPE VELOCITY),0,2,0,.333333333(.2 FOR ESCAPE),0,2,CONT,10,4(2 FOR ESCAPE),40,2,20,3:FORMULA IS X(1,0)=9/(2+COS X(2,0)) 6000 =========================================================================== 6008 T0=4:T1=6: 6009 INPUT "T=";T 6010 IF TT1 THEN Y=O:GOTO 12 6011 Y=1 6012 F=Y*(T+1-T0):PRINT"F=";F 6014 'PUT D.Es HERE 6019 GOTO 9:'DELETE THIS LINE WHEN SOLVING EQUATIONS 6900 'SINGLE BURST DRIVING FUNCTIONS FOR DIANALOG etc.Change Variables (other than T) not to clash with DIANALOG. 7000 =========================================================================== 7008 T0=3:T1=6:T2=1:N=1 7009 INPUT "T=";T 7010 IF T>=N*T1-T2 THEN N=N+1:Y=1 7011 IF T>T0-T2+(N-1)*T1 THEN Y=0:GOTO 13 7012 Y=1 7013 F=Y*(T+1+T2-(N-1)*T1):PRINT"F=";F 7014 'PUT D.Es HERE 7019 GOTO 9:'DELETE THIS LINE WHEN SOLVING EQUATIONS 7900 'REPETITIVE BURST DRIVING FUNCTIONS FOR DIANALOG etc.PRESET Y=0,1 IF INITIAL FUNCTION to be ON,OFF.Change Variables (other than T) not to clash with DIANALOG. 8000 ========================================================================== 8010 X(1,2)=K(0)*(1-X(1,0)*X(1,0))*X(1,1)-X(1,0) 8020 'For PHASE PLOT,replace DIANALOG LINE 750 with following line by deleting 8 from 8750. 8750 VY=317-K6*8-(YO+YM*X(1,1))*6.005 8900 'TO GRAPH VAN DER POL's 2nd order DIFF.EQUATION IN DIANALOG,ANSWER QUESTIONS WITH FOLLOWING NUMBERS:1,2,0,.01,1,1,2,2,30,5,0,5,25,2 9000 ========================================================================== 9010 K(1)=26:X(1,2)=-K(0)*X(1,1)-13*X(1,0)+K(1) 9020 'For TIME PLOT,replace DIANALOG LINE 750 with following line by deleting 9 from 9750. 9750 VY=317-K6*8-(YO-YM*T)*6.005 9900 'TO GRAPH DAMPED VIBRATING SPRING MASS 2nd order DIFF.EQUATION IN DIANALOG,ANSWER QUESTIONS WITH FOLLOWING NUMBERS:1,2,0,.025,0,0,1,2,0,20,0,8,47,1 10000 ========================================================================== 10010 ' IBM COMPATIBLE LISTING