diff --git a/doc/compute_heat_flux.html b/doc/compute_heat_flux.html
index 7a1ba30d71..2d11ac21f4 100644
--- a/doc/compute_heat_flux.html
+++ b/doc/compute_heat_flux.html
@@ -97,7 +97,7 @@ be multiplied by V/(kB T^2) times the sample interval and the
 appropriate unit conversion factors.  For real <A HREF = "units.html">units</A> in
 LAMMPS,
 </P>
-<PRE>lamda (KCal/(mol fmsec Ang K)) = V/(k_B*(T^2)) x "thermo" output frequency x timestep 
+<PRE>kappa (KCal/(mol fmsec Ang K)) = V/(k_B*(T^2)) x "thermo" output frequency x timestep x Integral value 
 </PRE>
 <P>where
 </P>
diff --git a/doc/compute_heat_flux.txt b/doc/compute_heat_flux.txt
index b3a4df786d..edacb42cef 100644
--- a/doc/compute_heat_flux.txt
+++ b/doc/compute_heat_flux.txt
@@ -94,7 +94,7 @@ be multiplied by V/(kB T^2) times the sample interval and the
 appropriate unit conversion factors.  For real "units"_units.html in
 LAMMPS,
 
-lamda (KCal/(mol fmsec Ang K)) = V/(k_B*(T^2)) x "thermo" output frequency x timestep :pre
+kappa (KCal/(mol fmsec Ang K)) = V/(k_B*(T^2)) x "thermo" output frequency x timestep x Integral value :pre
 
 where