%   
% Written by:
% -- 
% John L. Weatherwax                2006-08-28
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----

close all; clc; clear; 

allAs = [ -1,1 ];
allBs = [ -2,2 ];
allDeltaTs = [ 0.01, 0.1, 1. ];

qC = 1; 
L_Qc_LT = [ 0, 0; 0, qC ];
Nsteps = 100; 

fprintf('%10s %10s %10s: %10s %10s %10s\n','a','b','dt','q11','q12','q22');
for a = allAs
  for b = allBs
    for dt = allDeltaTs
      v = example_4_2_2_compute_qkm1( a, b, L_Qc_LT, dt, Nsteps );
      fprintf('%10.4f %10.4f %10.4f: %10.4f %10.4f %10.4f\n',a,b,dt, v(1,1),v(1,2),v(2,2));
    end
  end
end