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

close all; drawnow; clc; clear; clear functions; 

alpha = 0.1; w = 1 - alpha; 

for l = 1:4, 
  cl2 = 1 + alpha/( ( 1 + w )^3 ) * ( ( 1 + 4 * w + 5 * w^2 ) + 2 * l * ( 1 - w ) * ( 1 + 3 * w ) + 2 * l^2 * ( 1 - w^2 ) ); 
  fprintf( 'l = %10d; cl^2 = %10.6f; cl = %10.6f\n', l, cl2, sqrt(cl2) ); 
end