function f = KdVInit(x,c)
    global m;
%f = 0.5*c*(cosh(0.5*sqrt(c).*x)).^(-2);
    f = c./(2*cosh(0.5*sqrt(c).*x).^2);
end