

fi = field_integral_raw(0,0,4000/0.511,Bmap,0.003);

figure;set(gcf,'color','w');
plot(fi(:,1),fi(:,2)); hold on;
plot(fi(:,1),fi(:,3),'color','red'); hold off;
xlabel('z position (m)'); ylabel('second integral (T m^2)');
legend('x','y','location','best');
enhance_plot('arial',20,2,8)

figure;set(gcf,'color','w');
plot(fi(:,1),fi(:,4)); hold on;
plot(fi(:,1),fi(:,5),'color','red'); hold off;
xlabel('z position (m)'); ylabel('first integral (T m)');
legend('x','y','location','best');
enhance_plot('arial',20,2,8)
