subplot(1,1,1);
mesh(tau,fcvec,Ltaufc)

% Uncomment the line below to plot with weighting by p(tau, fc)
%mesh(tau,fcvec,Ltaufc.*ptau(tau,fcvec))

view(-20,60)
