h=0.1;
t = -1:h:0;
m = length(t);
x(m) = 0;
x(m-1) = -(-1)*h + x(m);

for n=1:m-2
	x(m-n-1) = 2*x(m-n) - x(m-n+1) - (h^2)*x(m-n)*(1+t(m-n)^2);
end
plot(t,x,'g');
axis ([-2, 2, -2, 2])

hold on;

h=0.1;
t = 0:h:1;
x(1) = 0;
x(2) = (-1)*h + x(1);

for n=2:length(t)-1
	x(n+1) = 2*x(n) - x(n-1) - (h^2)*x(n)*(1+t(n)^2);
end
plot(t,x,'g');
axis ([-2, 2, -2, 2])

h=0.02;
t1 = -1:h:0;
m = length(t1);
x1(m) = 0;
x1(m-1) = -(-1)*h + x1(m);
for n=1:m-2
	x1(m-n-1) = 2*x1(m-n) - x1(m-n+1) - (h^2)*x1(m-n)*(1+t1(m-n)^2);
end
plot(t1,x1,'r');
axis ([-2, 2, -2, 2])


h=0.02;
t1 = 0:h:1;
x1(1) = 0;
x1(2) = (-1)*h + x1(1);

for n=2:length(t1)-1
	x1(n+1) = 2*x1(n) - x1(n-1) - (h^2)*x1(n)*(1+t1(n)^2);
end
plot(t1,x1,'r');
axis ([-2, 2, -2, 2])
