function [x]=xfunct(t,H,deltaT) x=H/2*(1-cos(2*pi*t/deltaT));