%
%  We want to determine the distribution of eigenvalues
%  of the Jacobian at time 0
%
[tvec,x0,options] = poiseuille4(0,0,'init');

out3 = 

              AbsTol: []
                 BDF: []
              Events: []
         InitialStep: []
            Jacobian: []
           JConstant: []
            JPattern: []
                Mass: []
        MassConstant: []
        MassSingular: []
            MaxOrder: 2
             MaxStep: []
         NonNegative: []
         NormControl: 'on'
           OutputFcn: []
           OutputSel: []
              Refine: []
              RelTol: 1.0000e-005
               Stats: []
          Vectorized: []
    MStateDependence: []
           MvPattern: []
        InitialSlope: []

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