clear all; close all;

tol=10^(-5);
xspan=[-1 1];
n0=100;
beta0=99;
A=1;
x0=[0 A];

for modes=1:7
  beta=beta0;
  dbeta=1;
    
  for j=1:1000
     [t,y]=ode45('shoot2',xspan,x0,[],beta,n0);
 
     if abs(y(end,1))<tol
         break    
     end
     
     if (-1)^(modes+1)*y(end,1) > 0
         beta=beta-dbeta;
     else 
         beta=beta+dbeta/2;
         dbeta=dbeta/2;
     end
    
  end

  plot(t,y(:,1)), hold on
  beta
  beta0=beta-1;
  
     
end