function bessel_backward

h=0.01;
n=0;

t=0:h:10;
[l,m]=size(t);

x(:,1)=[0;1];

for k=2:m
    A=[1, -h*(n/t(k)-t(k));-h/t(k),1];
    x(:,k)=inv(A)*x(:,k-1);
end

figure(1)
plot(t,x(2,:),'b',t,bessel(n,t),'r')