First we initialize our initial value problem we want to solve.
odefun = @(t,y) cos(y).^2*sin(t);
tspan = [0 16*pi];
y0=1;
We then use the ode45 function without any specified options to solve this problem. To compare it later we plot the trajectory.
[t,y] = ode45(odefun, tspan, y0);
plot(t,...
