function y=lin1p(t,x)

y= -50*(x + cos(t));