This paper applies an ordinary differential equation (ODE) solver of interpolating matrix method to calculate large deflection of thin circular plate with varible thickness.