  We obtain a stringent upper limit on axion-photon coupling from the 1038 days
gamma-ray observations of the TeV blazar Markarian 421. The long-term VHE
gamma-ray spectra are measured by the collaborations Fermi-LAT and HAWC from
2015 June to 2018 July. We show the best-fit SEDs of Markarian 421 under the
null and axion hypotheses. Then we set the axion-photon limit in the $\{m_a, \,
g_{a\gamma}\}$ plane. The 99% $\rm C.L.$ upper limit set by Markarian 421 is
$g_{a\gamma} \lesssim 4.0\times 10^{-12} \rm \, GeV^{-1}$ for the axion mass
$[1.0\times10^{-9} \, {\rm eV} \lesssim m_a \lesssim 1.0\times10^{-8} \, {\rm
eV}]$. It is the most stringent upper limit in this axion mass region.
