function f = schrof(t,y)
f = [-2*y(2)*y(1);...
     y(1)^2+y(3)^2-y(2)^2-1;...
    -2*(y(2)+y(1))*y(3)];
