function f = peakfkt(x)

    f = 0.51+25*exp(-81*(11*x+45).^2);
    


