
%Tälläisiä siis pitäisi tulla
X=linspace(-1.5,10,100);
iA=airy(-X);
iB=airy(2,-X);
plot(X,real(iA),'k--',X,real(iB),'k.-.')

hold 

%kokeillaan, tuleeko
opt0=odeset;
ode45(@(t,x,a) te42d(t,x,2),[-1.5,50],[airy(0);airy(2,0)],opt0);


%Event
figure
opt5=odeset('event',@tapaus,'OutputFcn',@odephas2);
ode45(@(t,x,a) te42d(t,x,2),[-1.5,50],[airy(0);airy(2,0)],opt5);



