%% lambert_plot plots the lambert function.
%
a = -3.0;
b = 3.0;
x = linspace ( a, b, 101 );
y = lambert ( x );

clf ( );
hold ( 'on' );
plot ( x, y, 'linewidth', 3 );
plot ( [a,b], [0,0], 'k-', 'linewidth', 3 );
hold ( 'off' );

grid ( 'on' );
xlabel ( '<-- x -->' );
ylabel ( '<-- lambert(x) -->' );
title ( 'lambert function' );

print ( '-dpng', 'lambert.png' );
