function no_result = driverII()

% Simple driver for BC_caseII.m

% First, plot the function
omega = (0:100);
g = zeros(101,1);
for i = 1:101
	g(i) = BC_caseII(omega(i));
end

figure(3) 
hold on
plot(omega,g)
title('y(L) as a function of omega, Case II')
plot([0,100],[0,0],'r--')


% Then, find a root of BC_caseII
omega1 = bisect('BC_caseII', 30, 35)
omega2 = bisect('BC_caseII', 60, 65)
omega3 = bisect('BC_caseII', 90, 100)
w1 = rk2( [0,1,omega1], 0, 0.001, 1000, 'caseII' );
w2 = rk2( [0,1,omega2], 0, 0.001, 1000, 'caseII' );
w3 = rk2( [0,1,omega3], 0, 0.001, 1000, 'caseII' );
x = 0:0.001:1;
figure(4), title('First 3 eigenmodes y(x), Case II')
hold on
plot(x, w1(:,1),'r')
plot(x, w2(:,1),'g')
plot(x, w3(:,1),'b')

