simulation seems to work, found bug in ellipitical ints on Bz and Br