Write a MATLAB implementation that derives the numerical solution of IVP $($P$)$ by a third order RK method. The implementation should take the necessary values as input, assuming the exact solution is $y(t)=e2t,$ find the exact and approximate values at the mesh points and find the error by using the formula

$E=ma{x}_{0}kN|y_k-y(t_k)|$

Give the numerical analysis for the different values of step sizes and make the following table:

$$table$[[,$Error$],[,N=20,N=40,N=80,N=160,N=320$