harj21.m

function pal = harj21( x, y )
    global z
    pal = exp(-y/z)*sin(x);
end