# | 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 = 0.5* k0 * dr2; |
32 | > | c_potential_E = 0.5 * k0 * dr2; |
33 | force = - k0 * dr; | |
34 | return force; | |
35 |
– | Removed lines |
+ | Added lines |
< | Changed lines |
> | Changed lines |