Function to create damping matrix - Jeremy Kim

Needs to be tested.