function [xdot] = xdotfunc(t,H,deltaT)

xdot=H/2*(2*pi/deltaT)*sin(2*pi*t/deltaT);