M = 10;
C = 1;
K = 5;
F = 0;
xini = [2, 2];
h = 0.1;
N = 1000;

P = [(M/h^2 + C/h + K), (-2*M/h^2 - C/h) , M/h^2, F];
t = h:h:N*h;
x = Ecuacion_Segundo_Orden(xini, h, N, P);
plot(t,x);
