function xDot = stiffFcn2(t,Xvec)
xDot = [ Xvec(2); -1001*Xvec(2) - 1000*Xvec(1) ];
