added damping coeff: arguments are from,delta,plot = 0, damp = 0, which method = 0. plot=1 for plotting,whichmethod = 1 for legacy, damp = x for damping != 0