The mechanism of successive phase transitions occurring at antiferroelectric smectics is well described by ANNNI model. Here the transfer matrix method is generalised so as to be applied to this model with long range interaction. The correlation function is obtained, which shows oscillating behaviour with varying wave number due to the complex eigenvalue of transfer matrix. In the framework of this method the critical temperature of the model is estimated by replacing an intralayer interaction by the molecular field.