Function uses the third-order Runge–Kutta scheme

function y = odesRK3(f, t, y0)
n=length(t);
y=nan(length(y0),n);
y(:,1)=y0(:);
for k=1:n-1
    h=t(k+1)-t(k);
    F1=h*f(t(k),y(:,k));
    F2=h*f(t(k)+h/2,y(:,k)+(F1/2));
    F3=h*f(t(k)+0.75*h,y(:,k)+(0.75*F1));
    y(:,k+1)=y(:,k)+(2*F1+3*F2+4*F3)/9;
end

[Torsten]

And do you also have a question ?

[John D'Errico]

This must be like a quiz show. We are given the answer, and now we need to guess the question.

I will guess, "What color is the bluebird?"