1000 'DIFSIMEQ.TXT:TEST EQUATIONS FOR MERGING WITH DIFFSIMP.RENUMBER THE LINES TO SUIT BY REMOVING LEFT DIGITS OF EACH LINE NUMBER. 1010 PI=3.141592:F=(PI*SIN(5*X(0,0))/8-EXP(PI*SQR(X(0,0)+1)/8)+SQR(X(0,0)+PI*5/8))*SQR(1+PI*(X(0,0)-1)^2/8)/((EXP(-1*X(0,0)*X(0,0))-PI/4)*(1+PI*(X(0,0)+1)^2/8))+PI/8:RETURN:'1 ALG TRANSCENDENTAL:C=-1 2000 ------------------------------------------------------------------------- 2010 F=2*X(0,0)-3*X(1,0)-17:RETURN:'2 ALG LINS:C=0,0:X=4,-3 2011 F=X(0,0)+X(1,0)-1:RETURN 3000 ------------------------------------------------------------------------- 3010 F=X(0,0)^2+X(0,0)*X(1,0)+X(1,0)^2-3:RETURN:'2 ALG TRANS:C=0,0:Xo=1.999,1.0636;X1=-.9537,.93457: 01--,-+ 10-- 3011 F=SIN(X(0,0))-X(1,0)^2:RETURN 4000 ------------------------------------------------------------------------- 4010 F=X(0,1)-X(0,0)-T:RETURN:'1st ORDER:C(J)=0:X=EXP T-T-1=.718281828 @T=1 4011 F=X(1,5)-X(1,1)-1:RETURN:'5th ORDER:C(J)=0,0,1,1,1,:X=EXP T-T-1=.718281828 @T=1 5000 -------------------------------------------------------------------------- 5010 F=X(0,2)^2+X(1,0)^2+X(0,1)-EXP(2*T)+2*T-3:RETURN:'2nd+Oth ORDER:C=1,0: Xo=EXP T-t^2-=2.255252 ;X1=EXP(T/2)-EXP(-T/2)=1.042191 @T=1 5011 F=X(1,0)-SQR(X(0,1)+2*T+1):RETURN 6000 -------------------------------------------------------------------------- 6010 F=X(0,2)-SQR((2+X(1,0)*X(0,1))/8):RETURN:'CATENARY:2nd+1st ORDER:C=2,0,0: Xo=EXP(T/2)+EXP(-T/2)=2.255252 ;X1=EXP(T/2)-EXP(-T/2)=1.042191 @T=1 6011 F=X(1,1)-SQR(1+X(0,1)^2):RETURN 7000 --------------------------------------------------------------------------- 7010 F=X(0,2)-2*X(1,0)+5*X(0,0):RETURN:'2 SPRINGS:2 2nd ORDER LINS:C=5,0,0,0:Xo=4COS[T SQR 6]+COS T=-2.539321 @T=1 7011 F=X(1,2)-2*X(0,0)+2*X(1,0):RETURN:'X1=-2COS[T SQR 6]+2COS T=2.620416 @T=1 8000 --------------------------------------------------------------------------- 8010 F=X(0,2)-X(0,0)*X(1,1)^2+2/X(0,0)^2:RETURN:'ORBITS:2 2nd ORDERS:C=3,0,0,.33333333:Xo=9/[2+COS X1] 8011 F=X(1,2)+2*X(0,1)*X(1,1)/X(0,0):RETURN: 9000 --------------------------------------------------------------------------- 9010 GOSUB 740:F=X(0,2)-SIN((X(0,2)-(1+T)*SIN(T))/X(0,0)):RETURN:'IMPLICIT SIN 2nd ORDER:C=0,-1:X=-SIN T=-.8414709848 @T=1:Use VARIABLEs LIMITERS below on lines 740-741. 9740 IF X(0,0)=0 THEN X(0,0)=.1 9741 RETURN 10000 -------------------------------------------------------------------------- 10010 GOSUB 740:F=LOG(X(0,2))-(X(0,2)^2*LOG(X(0,0))+T*EXP(2*T))/(2*X(0,1)^2):RETURN:'IMPLICIT EXP 2nd ORDER:C=1,1:X=EXP T=2.718281828 @T=1:Use VARIABLEs LIMITERS below on lines 740-743. 10740 IF X(0,0)<.1 THEN X(0,0)=.1 10741 IF X(0,1)=0 THEN X(0,1)=.1 10742 IF X(0,2)<.1 THEN X(0,2)=.1 10743 RETURN 12000 ' IBM COMPATIBLE LISTING