function dy = dervect(t,y)
dy = zeros(2,1);    % a column vector
dy(1) = y(2);
dy(2) = sin(t)-10*y(1) -2*y(2);


