%  This function solves the fifth part of homework problem [H6.9]
%  --------------------------------------------------------------
%
echo on
diary nsds_hw6_9e.dry
%
%  We want to determine the distribution of eigenvalues
%  of the Jacobian at time 0
%
[tvec,x0,options] = poiseuille4(0,0,'init');
J0 = poiseuille4(0,x0,'jacobian');
l0 = eig(J0);
%
%  Plot the results
%
plot(real(l0),imag(l0),'x')
grid on
title('ttt')
xlabel('rellamb')
ylabel('imllamb')
print -deps nsds_hw6_9e.eps
%
diary off
echo off
%
return