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>

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