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')