Pa = 1000;
Pb = 500;
g = 9.8;
H = 1.0;

h = 0.1;
N = 1000;

xini = [0.01, 0.01];

P = [1, 0 , Pa*g/(Pb*H), 0];
t = h:h:N*h;
x = Ecuacion_Segundo_Orden(xini, h, N, P);
plot(t,x)