Descifrado del código fuente de Matlab
Función [] = SOLVE(x_0, y_0, A, k, N, step_t)
Función [ u2 ] = D( u1, A1, A2, h)
Función [ f ] = F( v, B1, B2)
v 1 = v(1);
v2 = v(2);
f 1 = b 1 B2 * pecado(v2-v 1);
F2 = b 1 B2 * pecado(v 1-v2);
f =[f 1;F2]; p> p>
Fin
K1 = F( u1, A1, A2
K2 = F( u1 h/2 * K1, A1, A2); /p >
K3 = F( u1 h/2 * K2, A1, A2
K4 = F( u1 h * K3, A1, A2); U2 = u 1 h/6 *(k 1 2 * K2 2 * K3 K4);
Fin
p _ 0 =[x _ 0]; /p >
PLOT = P_0;
Porque i = 1: 1: N
P_1 = D(P_0,A,k,step_t);
PLOT = [PLOT,P_1];
P_0 = P_1;
Fin
PLOT_X = PLOT( 1, );
PLOT_Y = PLOT(2,);
Figura (1);
plot(PLOT_X,PLOT_Y);
Fin