%eck.m
function V=eck(x)
a=1;b=8*a;
V=0.5*(a*exp(x)./(1+exp(x)) + b*exp(x)./(1+exp(x)).^2);
