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