  We investigate in detail the $\gamma\gamma$ fusion production mechanisms of
two neutral Higgs bosons ($h^0A^0$, $ H^0A^0$, $h^0H^0$ and $H^0H^0$) within
the framework of the mSUGRA-inspired minimal supersymmetric standard
model(MSSM) at an $e^+e^-$ linear colliders, which provide a probe of the
trilinear Higgs self-couplings. We calculate the dependence of the production
rates on Higgs boson masses, the ratio of the vacuum expectation values $\tan
\beta$ and the CMS energy $\sqrt{s}$. We find that the cross section for the
$H^0H^0$ production at LC can reach $0.2 fb$, while the cross section of
$A^0H^0$ production is only $10^{-4}\sim 10^{-3} fb$ under our parameters.
