10 '-----ピストン加速度の近似------
20 CLS 3:SCREEN 3
30 X=100:Y=200 : NN=6 : TIME=1000
40 PAI=3.141592653#
50 RAD=PAI/180
60 ZOOM=80
70 A=.7#
80 LINE(X-50,Y)-(X+400,Y)
90 LINE(X,Y-150)-(X,Y+150)
100 '-----厳密解-----
110 FOR T=0 TO 360 STEP 1
120 THETA=T*RAD
130 V1=-(COS(THETA)+(((A*COS(2*THETA))+(A^3*SIN(THETA)^4))/(SQR(1-(A*SIN(THETA))^2)^3)))
140 IF T=0 THEN 150 ELSE 160
150 PSET(X+T,Y-V1*ZOOM):GOTO 170
160 LINE -(X+T,Y-V1*ZOOM)
165 FOR W=0 TO TIME STEP 1 : NEXT W
170 NEXT T
180 LOCATE 15,3 : PRINT "Enter キーを押してください"
190 IF INKEY$ = "" THEN
200 GOTO 190
210 ELSE
220 GOTO 250 							
230 END IF 
240 '------近似解------
250 FOR T=0 TO 360 STEP 1
260 THETA=T*RAD
270 B(1)=(4*(((1/4)*A)+((1/16)*A^3)+((15/512)*A^5)+((35/2048)*A^7)+((735/65536)*A^9)+((2079/262144)*A^11))*COS(2*THETA))
280 B(2)=(16*(((-1/64)*A^3)-((3/256)*A^5)-((35/4096)*A^7)-((105/16384)*A^9)-((10395/2097152)*A^11))*COS(4*THETA))
290 B(3)=(36*(((1/512)*A^5)+((5/2048)*A^7)-((315/131072)*A^9)-((1155/524288)*A^11))*COS(6*THETA))
300 B(4)=(64*(((-5/16384)*A^7)-((35/65536)*A^9)-((693/1048576)*A^11))*COS(8*THETA))
310 B(5)=(100*(((7/131072)*A^9)-((63/524288)*A^11))*COS(10*THETA))
320 B(6)=(144*(((-21/2097152)*A^11))*COS(12*THETA))
330 C(1)=-COS(THETA)-B(1)
340 C(2)=-COS(THETA)-(B(1)+B(2))
350 C(3)=-COS(THETA)-(B(1)+B(2)+B(3))
360 C(4)=-COS(THETA)-(B(1)+B(2)+B(3)+B(4))
370 C(5)=-COS(THETA)-(B(1)+B(2)+B(3)+B(4)+B(5))
380 C(6)=-COS(THETA)-(B(1)+B(2)+B(3)+B(4)+B(5)+B(6))
385 IF I=0 THEN 390 ELSE 411
390 IF T=0 THEN 400 ELSE 410
400 PSET(X+T,Y-C(NN)*ZOOM),4:GOTO 420
410 LINE-(X+T,Y-C(NN)*ZOOM),4
415 'FOR W=0 TO TIME STEP 1 ; NEXT W
420 NEXT T
430 END