| # | 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 |