Matematické Fórum

Skrytý text:
Sub aaa()
Rem vystartujeme z R0=6378 000 m, rychlost v0=11200/2 m/s smerom hore
Rem kapa= 6.67 × 10-11 N m2 kg-2 ,,Hmotnosť Zeme je 5974. 1024 kg,
Rem casovy interval vypoctu dt=1 sek, zastavit pri tstop ..pozor vyrobi sheets s tymto poctom riadkov (tstop)
Rem s=6378000+v0*t-(1/2)*a*t*t... a= kapa*Mz/(r*r)
Rem prvy interval, na zaciat v0=v1, a1=kapa*Mz/(r0*r0),s=v1*t-(1/2)*a1*t*t, r1=r0+s, v2=v1-a1t
Rem druhy interval, v2, a2=kapa*Mz/(r1*r1),s=v2*t-(1/2)*a2*t*t,, r2=r1+s, v3=v2-a2t
tstop = 10000: Dim v(100000): Dim a(100000): Dim r(100000): v(0) = 11200 / 2: r(0) = 6378000: kapa = 6.67 * 10 ^ -11: mz = 5.974 * 10 ^ 24
dt = 1: Sheets.Add
For i = 0 To tstop
a(i + 1) = kapa * mz / (r(i) ^ 2): s = v(i) * dt - (1 / 2) * a(i + 1) * dt ^ 2: r(i + 1) = r(i) + s: v(i + 1) = v(i) - a(i + 1) * dt
Cells(i + 1, 1) = i: Cells(i + 1, 2) = a(i): Cells(i + 1, 3) = v(i): Cells(i + 1, 4) = r(i)
Next i

End Sub