моделирование движения сверхзвукового самолета

[Азаматkai]

дана блок схема для фортрана, надо написать программу в паскале


F(V,θ)=(g(T-Wsinθ-D))/WVsinθ; (2.4)
D=1/2 C_D ρSV^2; (2.5)
G(V,θ)=(g(L-Wcosθ))/(WV^2 sinθ); (2.6)
L=1/2 C_L ρSV^2.


├ █(█(█(g=32.174 фут/〖сек〗^2 -ускорение свободного падения@W=20000 фунтов -вес самолета)@C_D=0.02-коэффициент сопротивления@C_L=0.1-коэффициент подъемной силы)@S=160 〖фут〗^2-площадь крыльев)} (2.8)
(T=(1+1/2 e^(-(V•〖10〗^(-3)-1.5)^2 ) )(〖10〗^4-0.27Z) фунтов-тяга двигателя (2.9)
├ █(ρ=ρ_0 e^(-Z/(2.36•〖10〗^4 )) фунт•с^2/〖фут〗^4@ρ_0=2,3769•〖10〗^(-3 ) фунт•с^2/〖фут〗^4 )} плот



Исходя из приведенных соотношений, процесс решения сводится к следующему:
Используется исправленный метод Эйлера для вычисления значений V_1, θ_1.
V_1=V_0+h/2 [F(V_0,θ_0 )+F(V_0+hF(〖V_0,θ〗_0 ),θ_0+hG(V_0,θ_0))]; (2.11)
θ_1=θ_0+h/2 [G(V_0,θ_0 )+G(V_0+hF(〖V_0,θ〗_0 ),θ_0+hG(V_0,θ_0 ) ) ]. (2.12)
Для m=1,2,… значения V_(m+1) и θ_(m+1) предсказываются с помощью формул
V_(m+1)^((0))=V_(m-1)+2hF(V_m,θ_m ), (2.13)
θ_(m+1)^((0))=θ_(m-1)+2hF(V_m,θ_m ). (2.14)

Коррекция значений V_(m+1) и θ_(m+1):
V_(m+1)^((i))=V_m+h/2 [F(V_m,θ_m )+F(V_(m+1)^((i-1) ),θ_(m+1)^((i-1)))], (2.15)
θ_(m+1)^((i))=θ_m+h/2 [G(V_m,θ_m )+G(V_(m+1)^((i-1) ),θ_(m+1)^((i-1)))]. (2.16)

Оценка ошибки ограничения
E_T (V_(m+1) )=1/5 (V_m^((0) )-V_m^((i) ) ), (2.17)
E_T (θ_(m+1) )=1/5 (θ_m^((0) )-θ_m^((i) ) ). (2.18)

[solivan]

Блок-схема - это хорошо, а что является исходными данными, что за исправленный метод Эйлера, корреляции, оценки и т.д. и т.п.
Как-то у меня много вопросов...

[Азаматkai]

вот все формулы