% Define right hand side in ODE y' = f(t,y)
f = @(t,y) sin(t^2 + 1/sqrt(abs(y)));

% Define initial and final times
tmin=0;
tmax=20;

% Define initial value
yzero=1;

% Optional: Choose number of time steps N
% N = 100000;
% h = (tmax-tmin)/N;

% Solve ODE
[t y] = ode45(f,[tmin,tmax],yzero);

% Optional: Specify number of timesteps
% [t y] = ode45(f,[tmin:h:tmax],yzero);

% Plot results
plot(t,y), xlabel('t'), ylabel('y')