I need a solution handwritten, in matlab, or mathematica, thank you!

1 (a). Use the Runge-Kutta method for systems to approximate the solutions of the following systems of first-order differential equations and compare the results to the actual solutions. (a) u1=3u1+2u2(2t2+1)e2t, for 0t1, with u1(0)=1 u2=4u1+u2+(t2+2t4)e2t, for 0t1 with u2(0)=1;h=0.2 actual solutions u1(t)=31e5t31et+e2t and u2(t)=31e5t+32et+t2e2t