function phi2 = flux2(h,z)
% Problem 2.a second function
A1=100; Bz1=0.0159; Bz2=0.0321; H=380;
A2 = A1 * sin(Bz1*h) / sinh(Bz2*(H - h));
phi2 = A2 * sinh(Bz2*(H-z));