mirror of
https://github.com/mbirth/gwbasic.git
synced 2024-09-19 16:53:26 +01:00
49 lines
1.2 KiB
QBasic
49 lines
1.2 KiB
QBasic
10 'Programm: SONNE.BAS
|
||
20 'Funktion: Simulation der Planeten-
|
||
30 ' gravitation
|
||
40 'Sprache: GW-Basic
|
||
50 'Autor: Axel Werner
|
||
60 '(C)1991 DMV Widuch GmbH & Co. KG
|
||
70 DEFINT I,S,A,F
|
||
80 INPUT "A: ",A
|
||
90 IF A>5 THEN A=5
|
||
100 IF A<1 THEN A=1
|
||
110 FOR I=1 TO A
|
||
120 PRINT"XP("I"),YP("I"),ZP("I"): ";
|
||
130 INPUT"",XP(I),YP(I),ZP(I)
|
||
140 PRINT"DXP("I"),DYP("I"),DZP("I"): ";
|
||
150 INPUT"",DXP(I),DYP(I),DZP(I)
|
||
160 PRINT"M("I"): ";:INPUT"",M(I)
|
||
170 NEXT
|
||
180 INPUT"MS: ",M(0)
|
||
190 XP(0)=0:YP(0)=0:ZP(0)=0
|
||
200 SCREEN 2:CLS:KEY OFF
|
||
210 WINDOW SCREEN(-320,-250)-(320,250)
|
||
220 G=.006672
|
||
230 FOR I=0 TO A
|
||
240 NXP=XP(I)+DXP(I):NYP=YP(I)+DYP(I)
|
||
250 NZP=ZP(I)+DZP(I)
|
||
260 FOR S=0 TO A
|
||
270 IF I=S THEN 350
|
||
280 X2=XP(S)-XP(I):Y2=YP(S)-YP(I)
|
||
290 Z2=ZP(S)-ZP(I)
|
||
300 D=SQR(X2*X2+Y2*Y2+Z2*Z2)
|
||
310 GK=G*M(S)*M(I)/(D*D)
|
||
320 B=GK/(2*M(I))
|
||
330 X1=X2*B/D:Y1=Y2*B/D:Z1=Z2*B/D
|
||
340 NXP=NXP+X1:NYP=NYP+Y1:NZP=NZP+Z1
|
||
350 NEXT
|
||
360 DXP(I)=NXP-XP(I):DYP(I)=NYP-YP(I)
|
||
370 DZP(I)=NZP-ZP(I):GOTO 400
|
||
380 X%=XP(I):Y%=YP(I):Z%=ZP(I)
|
||
390 F=0:GOSUB 450
|
||
400 X%=NXP:Y%=NYP:Z%=NZP:F=1:GOSUB 450
|
||
410 XP(I)=NXP:YP(I)=NYP:ZP(I)=NZP
|
||
420 NEXT
|
||
430 IF INKEY$="" THEN 230
|
||
440 END
|
||
450 X%=X%+Z%/3
|
||
460 Y%=Y%-Z%/3
|
||
470 PSET(X%,Y%),F
|
||
480 RETURN
|
||
|