Código fuente de Vmatlab

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];

Fin

K1 = F( u1, A1, A2

K2 = F( u1 + h/2 * K1, A1, A2);

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); ;y_0];

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