Coulson (ZINDO), Mulliken $(MP2/6-31G^*)$ and Natural $(MP2/6-31G^*)$ population analyses of several large molecules were performed by the Fragment Reassociation (FR) method. The agreement between the conventional ZINDO (or conventional MP2) and FR-ZINDO (or FR-MP2) charges of these molecules was excellent. The standard deviations of the FR-ZINDO net atomic charges from the conventional ZINDO net atomic charges were 0.0008 for $C_{10}H_{22}$ (32 atoms), 0.0012 for $NH_2-C_{16}O_2H_{28}-COOH$ (53 atoms), 0.0014 for $NH_3^+-C_{16}O_2H_{28}-COOH$ (54 atoms), 0.0017 for $NH_2-C_{16}O_2H_{28}-COO^-$ (52 atoms), 0.0019 for $NH_3^+-C_{16}O_2H_{28}-COO^-$ (53 atoms), 0.0024 for a conjugated model $(O=CH-(CH=CH)_{15}-C=O-(CH=CH)_{12}-CH=CH_2)$, 118 atoms), 0.0038 for aglycoristocetin $(C_{60}N_7O_{19}H_{52}^+$, 138 atoms), 0.0023 for a polypropylene model complexed with a zirconocene catalyst $(C_{68}H-{121}Zr^+$, 190 atoms) and 0.0013 for magainin $(C_{112}N_{29}O_{28}SH_{177}$, 347 atoms), respectively. The standard deviations of the FR-MP2 Mulliken (or Natural) partial atomic charges from the conventional ones were 0.0016 (or 0.0016) for $C_{10}H_{22}$, 0.0019 (or 0.0018) for $NH_2-C_{16}O_2H_{28}-COOH$ and 0.0033 (or 0.0023) for $NH_3^+-C_{16}O_2H_{28}-COO^-$, respectively. These errors were attributed to the shape of molecules, the choice of fragments and the degree of ionic characters of molecules as well as the choice of methods. The CPU time of aglycoristocetin, conjugated model, polypropylene model complexed with zirconocene and magainin computed by the FR-ZINDO method was respectively 2, 4, 6 and 21 times faster than that by the normal ZINDO method. The CPU time of $NH_2-C_{16}O_2H_{28}-COOH\;and\;NH_3^+-C_{16}O_2H_{28}-COO^-$ computed by the FR-MP2 method was, respectively, 6 and 20 times faster than that by the normal MP2 method. The largest molecule calculated by the FR-ZINDO method was B-DNA (766 atoms). These results will enable us to compute atomic charges of huge molecules near future.