% 
% Written by:
% -- 
% John L. Weatherwax, Ph.D.         2003-01-28
% 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.
% 
%-----

x0 = 3.9;
u0 = 0.01;

x = linspace( eps, 1, 100 );
u = u0*( (x0./x).^12 - 2*(x0./x).^6 );

figure;
plot( x, u, '-o' ); title( 'Lennard-Jones potential' );