function no_result = driverI()

% Simple driver for BC_caseI.m

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

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


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

