function E = chainEnergy(p)

E = 

end
