Calculate the Flux of a PDE solution

Hi there,

In previous versions it was possible to calculate the flux of a PDE solution c*grad(u) with the command:

[CGXU,CGYU]=pdecgrad(P,T,C,U)

In version 2016, the help leads one to:

[gradx,grady] = evaluateGradient(results,xq,yq);

but there isn't an example to estimate the flux. Any ideas on how that can be done?

Thanks!