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

close all; clc; clear all; 

dOverEta = linspace( 0, 5, 100 );

qs = [ 1.1, 1.2, 1.5, 1.8, 2, 3, 5 ]; 

figure(); hold on;
for ii=1:length(qs)
  q   = qs(ii);
  uij = 1 ./ ( 1 + dOverEta.^(1/(q-1)) );

  plot( dOverEta, uij, '-' );
end

% Do the exponential decay in uij:
uij = exp( - dOverEta );

plot( dOverEta, uij, '-r', 'linewidth', 3 );

xlabel('d(x,theta)/eta'); ylabel('uij');

%saveas( gcf, '../../WriteUp/Graphics/Chapter14/chap_14_prob_17_plot.eps', 'epsc' );