Aproximacion de la reentrada en la atmosfera

BREVE MEMORIA

a) LA OBTENCIÓN DE LAS ECUACIONES ADIMENSIONALES
          b) EL PROCESO DE DISCRETIZACIÓN Y LA OBTENCIÓN DE LAS ECUACIONES DEL SISTEMA F(X)=0
          c) EL CÁLCULO DELA MATRIZ JACOBIANA Y LA ESTRATEGIA DE ASIGNACIÓN DE LOS ELEMENTOS DE LA MISMA

a) Para la obtención de las ecuaciones adimensionales se han tomado valores caracteristicos de longitud ,velocidad y tiempo.
lc = R
Vc = (g0R)1/2
tc = (R/g0)1/2

Utilizando este cambio de variable las ecuaciones se simplifican notablemente ayudando obeter su solución.b) En el proceso de discretizacion y obtención del sistema de ecuaciones se parte de el sitema dado en el enunciado del problema. A partir de ese punto se integran las ecuacions diferencialesutilizando el método del trapecio simple de forma que se obtienen cuatro ecuaciones, una por cada variable, siendo las tres restantes análogas a la siguiente:
F(1) = x(1) - v(1) -h*(f1(x,B,β,λ)+f1(v,B, β, λ))/2

c) El cálculo de la matriz jacobiano se ha realizado a mano la ecuaciones del sistema F(X)=0. Para introducir los elementos en la matriz se han asignado uno a uno los valores acada termino de la matriz
J(i,J)=…

Programa principal

Pseudocódigo

Datos: kmax, n, B, Beta, Landa, Grados, At, E
Resultados: U(n)
Real: U(n), B, beta, landa, G, F(n),grados, At, X(n), V(n), E
Enteros: n, i, kmax
Subprogramas: FUNC(X,V,At,B,Landa,Beta,F), JACOB(X,At,B,Landa,Beta,J), SUB_SIST(X,V,At,B,Landa,Beta,kmax,E),

kmax=10000
n=4
B =1.518d4
Beta = 9.52d2
Landa = 3d-1
Grados = -7d0 %Cambiamos el valor inicial de Gamma aquí
G = Grados*arccos(-1d0)/180d0
At = 1d-3
E = 1d-5

i = 0

U(1) = 1.32743278d0U(2) = G
U(3) = 0.0188146d0
U(4) = 0d0

X = U
V = U

Escribir 't', 'U(1)', 'U(2)', 'U(3)', 'U(4)'

i = i + 1

hacer
llamar_a SUB_SIST(X,V,At,B,Landa,Beta,kmax,E)...