We have investigated the numerical methods to calculate model atmosphere for the analysis of spectral lines emitted from the sun and stars. Basic equations used in our calculations are radiative transfer, statistical equilibrium and charge-particle conservations. Transfer equation has been solved to get emitting spectral line profile as an initial value problem using Adams-Bashforth-Moulton method with accuracy as high as 12th order. And we have calculated above non linear differential equations simultaneously as a boundary value problem by finite difference method of 3 points approximation through Feautrier elimination scheme. It is found that all computing programs coded by above numerical methods work successfully for our model atmosphere.