%the shooting method on undamped h.o.
clear all
close all
clc

A=3;
b=10;
y0=[1 A];
tspan=[0 b];
yb=2;
tol=1e-7;
[T Y]=ode45('ho',tspan,y0);

high=0;
low=3;

for i=1:10000
    mid=(high+low)/2;
    y0=[1 mid];
    [T Y]=ode45('ho',tspan,y0);
    if abs(yb-Y(end,1))<tol
        break;
    else
        if Y(end,1)<yb
            low=mid;
        else
            high=mid;
        end
    end
    Y(end,1)
end