function y=phi(x)
    y=exp(-x);
end