# que pasa si Delta t no tienda a cero sino a 1
# reescribir esto para otra deltat


clear all
clc
close all

N = 500;
x1 = zeros(N, 1);

x1(1)     = 1;
x1(1 + 1) = 2 + x1(1);

n = 1;


while n <= N - 2

  x1(n + 2) = 2*x1(n+1) - 2*x1(n)

  n += 1;
end

plot(x1)