Med energikilde langs z-aksen 

 

 

Plot 


Samlet energioptag for hele bevægelsen :
 

> Int(Int(Int(add(Jp[a], a = 1 .. 5), u = 0 .. 1), v = 0 .. 1), t = 0 .. Pi); -1