10 CLS 1
20 SCREEN 3
30 X0=100:Y01=100:Y02=270:Z=70
40 INPUT "αの値";A
50 LINE (X0,Y01-85)-(X0,Y01+85),7
60 LINE (X0,Y02-80)-(X0,Y02+120),7
70 LINE (X0,Y01)-(X0+400,Y01),7
80 LINE (X0,Y02)-(X0+400,Y02),7
90 LINE (X0+180,Y01-5)-(X0+180,Y01+5)
100 LINE (X0+360,Y01-5)-(X0+360,Y01+5)
110 LINE (X0+180,Y02-5)-(X0+180,Y02+5)
120 LINE (X0+360,Y02-5)-(X0+360,Y02+5)
130 FOR I=1 TO 360 
140 THETA=I*3.141592#/180
150 '速度計算と描画
160 Y1=(-1)*SIN(THETA)-((A/2)*SIN(2*THETA))/SQR(1-(A*SIN(THETA))^2)
170 PSET (X0+I,Y01-Y1*Z),7
180 'LINE -(X0+I,Y01-Y1*Z),7
190 '加速度計算と描画
200 Y2=(-1)*COS(THETA)-((A*COS(2*THETA)+A^3*SIN(THETA)^4)/SQR(1-(A*SIN(THETA))^2)^3)
210 PSET (X0+I,Y02-Y2*Z),7
220 'LINE -(X0+I,Y02-Y2*Z),7
230 NEXT I
240 WIDTH 80,25
250 '速度の文字表示
260 LOCATE 4,1:PRINT"Xp/(rθ)"
270 LOCATE 4,0:PRINT".     ."
280 LOCATE 11,7:PRINT"0"
290 LOCATE 34,7:PRINT"π"
300 LOCATE 56,7:PRINT"2π"
310 LOCATE 62,7:PRINT"θ"
320 '加速度の文字表示
330 LOCATE 2,12:PRINT"Xp/(rθ^2)"
340 LOCATE 2,11:PRINT"..    ."
350 LOCATE 11,17:PRINT"0"
360 LOCATE 34,17:PRINT"π"
370 LOCATE 56,17:PRINT"2π"
380 LOCATE 62,17:PRINT"θ"