function phi1 = flux1(z)
%  Problem 2.a first function
A1=100; Bz1=0.0159;
phi1 = A1 * sin(Bz1*z);