# | Line 29 | Line 29 | double HarmonicBond::bond_force( double r_ab ){ | |
---|---|---|
29 | dr = r_ab - d0; | |
30 | dr2 = dr * dr; | |
31 | ||
32 | < | c_potential_E = k0 * dr2; |
33 | < | force = -2.0 * k0 * dr; |
32 | > | c_potential_E = 0.5* k0 * dr2; |
33 | > | force = - k0 * dr; |
34 | return force; | |
35 | ||
36 | } |
– | Removed lines |
+ | Added lines |
< | Changed lines |
> | Changed lines |