function xdot=insulin(t,x)                     
xdot=zeros(2,1)
xdot(1)=x(2)^2-x(1)*x(2);
xdot(2)=-x(1)*x(2)+1;
endfunction

x0=[1;10];      // Anfangswert
t0=0;           // Anfangszeitpunkt    
t=[0:0.01:20];  // Simulationsintervall, das ausgegeben werden soll (muss hinter t0 liegen)

x=ode(x0,t0, t, insulin); // Simulationsrechnung erzeugt Matrix x, die die gesamte Trajektorie enthaelt.

plot2d(t,x')