We determine the stability of an axisymmetric geometrically thin and optically thickrnaccretion disc around a magnetized neutron star in the region where gas pressure isrndominant. Analysis is made with the inclusion of pressure from coulombic sources. Thernopacity in the middle region is mainly due to electron scattering whereas that in thernouter region is mainly due to free-free emission. Starting from the vertically integratedrnnon-relativistic hydrodynamics equations we set up the basic equations which govern thernstructure of the disc and for the stability analysis of the disc model we have kept the timerndependencies in the equations.rnAlthough we include the e ect of coulombic pressure it is thermally stable as if there wererngas pressure only and our graphical solution given at the end of this manuscript showsrnthat the disc is viscously unstable