Write an algorithm for RK$4$ method to find the approximate solution of the IVP corresponding to a system of $1$st order linear ODE $(\}.$Write a MATLAB code to solve the following IVP by taking the step parameter

$y_1^{\text{'}}(t)=y_2(t)$

$y_2^{\text{'}}(t)=-0\mathrm{.}75y_1(t)-2y_2(t)$

with initial value $y_1(0)=3$ and $y_2(0)=-2\mathrm{.}5$ in $0,1.$ Demonstrate MATLAB code by plotting the approximate and exact solutions at the respective node points.