genD function

Generate Bates and Watts D Matrix