function g = BC_caseI( omega ) % % return the value of w_1(x) at x = L = 1 meter. % L = 1; % lenght of domain n_steps = 1000; % number of steps to be used dx = L/n_steps; % set initial condtions: w_1(0) = phi(0) = 0 % w_2(0) = phi'(0) = 1 % w_3(0) = omega w0 = [ 0.0, 1.0, omega ]; w = rk2( w0, 0.0, dx, n_steps, 'caseI' ); % solve ODE g = w(n_steps+1,1); % w_1(L) = g