%openf91, system av diff.ekv, ode23
t0=0;tslut=10;
y0=[1 0]';
[T,Y]=ode23(@fosc,[t0,tslut],y0);
plot(T,Y)

