% 
% Written by:
% -- 
% John L. Weatherwax, Ph.D.         2003-01-15
% MIT Lincoln Laboratory
% 244 Wood Street
% Lexington, MA 02420-9176 USA
% 
% email: weatherwax@ll.mit.edu
% Voice: (781) 981-5370       Fax: (781) 981-0721
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----

b = linspace( eps, 10, 100 );
jArray = [ 1/2, 3/2, 5/2 ];

close all;

for j=1:length( jArray ),
  jPlusOneHalf = jArray(j)+1/2;
  mOverNDelta = jPlusOneHalf*coth( b*jPlusOneHalf ) - (1/2)*coth( b/2 );
  figure;
  plot( b, mOverNDelta, '-o' ); title( '\frac{M}{N \delta_{\mu}' );
end