function f = crps(mu,sig,x);
% calculate the continuous ranked probability srores for normal predictive
% distribution
xscaled = (x-mu)/sig;
fv = sig*(1/sqrt(pi)-2*normpdf(xscaled)-xscaled*(2*normcdf(xscaled)-1));
f = -fv;
