% Guo and Julien's method, in the closure
% Example in the Closure

t0 = cputime;

h = 2.38; % m/s

M = 0.2171; % m/s
N = 1.8321; % m/s
K = 0.1775; % kg/m^3
z = 0.3469;

dm = (0.105 + 0.074)./2.*1e-3; % 0.0895e-3 m
E = 2.*dm./h; % 0.07521

J = einstein(z,E); 

qT = K.*M.*h.*(4.61.*(1-E).^z./E.^(z-1) + N./M.*J(1) + J(2)) % kg/s.m

t = cputime - t0