function [z]=jetvari(k,om,U,h)

z=((om-0*k).*(1+tanh(k*h))-(k*U)).^2+(k*U).^2.*(tanh(k*h));