function [t,l,E]=runchua(a)
% by Dr. Ute Feldmann (26.10.12)
% runcchua(9);% 9-FP,11-P1,13.5-P2,14-Scroll,15-Doublescroll
%a=9;b=14.286;
b=28;
m0=-1/7;m1=2/7;
%y0=[rand*30+5,rand*35-30,rand*40-5];
y0=[.1,0,0];
options = odeset('OutputFcn',@odephas3);%3D-Ausgabe
%options = odeset('OutputFcn',@odeplot);%Zeitreihen-Ausgabe
[t,l]=ode23(@chuaD,[0 100],y0,options,a,b,m0,m1);%mit graphischer Ausgabe
%[t,l]=ode23(@chuaD,[0 100],y0,[],a,b,m0,m1);%keine graphische Ausgabe
pause(1)
plot(t,l(:,[1,2]));
legend('Uc1','Uc2')
pause(1)
ls=l(end-900:end,:); % nur das Ende der Zeitreihe
plot3(ls(:,1),ls(:,2),ls(:,3));
xlabel('Uc1');ylabel('Uc2');zlabel('Il');

% --------------Eigenwerte in jeder linearen Region angucken
A0=[-a*m0 a 0;1 -1 -1;0 b 0]
Eigenwerte0=eig(A0)
A1=[-a*m1 a 0;1 -1 -1;0 b 0]
Eigenwerte1=eig(A1)