function [t, y] = euler_solver(f, dt, tf, y0) % find number of steps, round up using ceil() % discretize time, t, using number of steps % define y vector using initial condition % loop through time and perform Euler's method end