Runge kuta
Enviado por Ninoka • 27 de Noviembre de 2018 • 6.791 Palabras (28 Páginas) • 309 Visitas
...
> k3:=f(X[1]+h,Y[1]+0.67*k2*h);
> k4:=f(X[1]+h,Y[1]+k3*h);
> Y[2]:=Y[1]+0.125*(k1+3*k2+3*k3+k4)*h;
k1 := 0.8285531011
k2 := 1.047494400
k3 := 1.047494400
k4 := 1.047494400
Y[2] := 0.01827353671
> X[1];
> Y[1];
> k1:=f(X[1],Y[1]);
0.1600000000
0.008072269328
k1 := 0.8285531011
> X[1]+h;
> Y[1]+0.5*k1*h;
> k2:=f(X[1]+h,Y[1]+0.5*k1*h);
0.1700000000
0.01221503483
k2 := 1.047494400
> X[1]+h;
> Y[1]+0.67*k2*h;
> k3:=f(X[1]+h,Y[1]+0.67*k2*h);
0.1700000000
0.01509048181
k3 := 1.047494400
> X[1]+h;
> Y[1]+k3*h;
> k4:=f(X[1]+h,Y[1]+k3*h);
0.1700000000
0.01854721333
k4 := 1.047494400
> 0.125*(k1+3*k2+3*k3+k4);
> Y[2]:=Y[1]+0.125*(k1+3*k2+3*k3+k4)*h;# Altura del espesador
# en la segunda iteración.
1.020126738
Y[2] := 0.01827353671
> XY:=[[X[0],Y[0]],[X[1],Y[1]],[X[2],Y[2]]]:
>
> X[3]:=X[2]+h; # Calculo de la concentración correspondiente al
# tercer intervalo o tercera iteración.
X[3] := 0.1800000000
> k1:=f(X[2],Y[2]);
> k2:=f(X[2]+h,Y[2]+0.5*k1*h);
> k3:=f(X[2]+h,Y[2]+0.67*k2*h);
> k4:=f(X[2]+h,Y[2]+k3*h);
> Y[3]:=Y[2]+0.125*(k1+3*k2+3*k3+k4)*h;
k1 := 1.047494400
k2 := 1.328888147
k3 := 1.328888147
k4 := 1.328888147
Y[3] := 0.03121067600
> X[2];
> Y[2];
> k1:=f(X[2],Y[2]);
0.1700000000
0.01827353671
k1 := 1.047494400
> X[2]+h;
> Y[2]+0.5*k1*h;
> k2:=f(X[2]+h,Y[2]+0.5*k1*h);
0.1800000000
0.02351100871
k2 := 1.328888147
> X[2]+h;
> Y[2]+0.67*k2*h;
> k3:=f(X[2]+h,Y[2]+0.67*k2*h);
0.1800000000
0.02717708730
k3 := 1.328888147
> X[2]+h;
> Y[2]+k3*h;
> k4:=f(X[2]+h,Y[2]+k3*h);
0.1800000000
0.03156241818
k4 := 1.328888147
> 0.125*(k1+3*k2+3*k3+k4);
> Y[3]:=Y[2]+0.125*(k1+3*k2+3*k3+k4)*h;# Altura del espesador en la
# tercera iteración.
1.293713928
Y[3] := 0.03121067600
> XY:=[[X[0],Y[0]],[X[1],Y[1]],[X[2],Y[2]],[X[3],Y[3]]]:
>
> X[4]:=X[3]+h; # Calculo de la concentración correspondiente al
# cuarto intervalo o cuarta iteración.
X[4] := 0.1900000000
>
...