function dydx = f(y,t)
 dydx = (sin(t))^2-cos(t);
endfunction