The code is clearer now and we have removed all data structures that were not used. Also the bending energy is still not right -> Eb(0) is not 0 as it should be by construction.