Modified-euler method and midpoint method

Question

dy/dt = te^3t - 2y, 0 ≤ t ≤ 1, y(0) = 0

Approximate the solution to the above initial value problem using

(a) Modified-Euler Method

(b) Midpoint Method

(c) Heun's Method

(d) 4-stage Runge-Kutta Method

with a time step of h = 0.1. For each method, tabulate the approximate solution and the error at each time step using the exact solution to this problem

y(t) = 1/5 te^3t -1/25 e^3t +1/25e^-2t

Plot the approximate solutions and the exact function for comparison. In another graph, plot the error distributions. Note that you may need to give a separate plot for the error of 4-stage Runge-Kutta scheme to see the trend in a larger scale.