function w = h(t,y)
w = 2*y -1;