function f=gravityspring(t,a)  %a is a 2x1 column vector

%{
Function for computation system
version of y"+9y=-1
%}

f(1,1)=a(2);
f(2,1)=-9*a(1)-1;