function [hl,hv]=enth(t,x,y)
global nc c1 c2 c3
hl=0;
hv=0;
for j=1:nc
    hl=hl+x(j)*c3(j)*t;
    hv=hv+y(j)*(c1(j)+c2(j)*t);
end
   