We present a non-perturbative calculation for the pion decay constant with quenched Kogut-Susskind quarks. Numerical simulations are carried out at $beta = 6.0$ and 6.2 with various operators extending over all flavors. The renormalization correction is applied for each flavor by computing non-perturbative renormalization constants, and it is compared with a perturbative calculation. We also study the behavior of $f_pi$ in the continuum limits for both non-perturbative and perturbative calculations. The results in the continuum limit is also discussed.