function dy  = F(t,y)
w  = 2.*pi;
dy = (0 1 ; -w 0);