Extrap1

clear all;
k = 2; F = [1,2,3];
a = 0; b = pi; nSteps = 100; tau = (b-a)/nSteps; x0 = 2;
t_cur = a; x_cur = x0;
x_soln = zeros(nSteps+1,1);
x_soln(1) = x_cur;

for step = 1:nSteps

    x_cur = psi_Extrap(t_cur, x_cur, tau, k, F);
    x_soln(step+1) = x_cur;
    t_cur = t_cur + tau;

end

figure(1)
plot((a:tau:b),x_soln); axis equal;