varD function

Dominance variance