RESIDUALVARIANCE scan with a single dipole and computes the RV
 at each grid location.

 Use as
   [dipout] = residualvariance(dip, grad, headmodel, dat, ...)