%this function represents a possible differential equation system. You will
%have to modify this to suit your needs for this project.

function xprime=diff_eq(t,x,m,k,a)
xprime(1,1)=m*x(2);
xprime(2,1)=-k*x(1)+a*x(2);