Solve using python code. Sympy, Numpy.

Define a function, forward_euler, that uses forward Euler to solve ODEs of the form: y' + py = q and accepts the following arguments in order specified: p is the coefficient of y(x), q is right-hand side of differential equation, y is the initial value of y at t = 0 (default should be 0.0), N is the number of equally-spaced times after t = 0 to be used for the approximation (default should be 25), and t_max is the end point for approximation (default should be 10.0). The function should return y, an np. array of the approximate y(t) values, and t, an np. array of the correspond times at which y(t) is evaluated