To improve the current status of home multilayer optical coating design with low speed and poor efficiency when a large layer number occurs, the accurate calculation and fast realization of merit function's gradient and Hesse matrix is pointed out. Based on the matrix method to calculate the spectral properties of multilayer optical coating, an analytic model is established theoretically. And the corresponding accurate and fast computation is successfully achieved by programming with Matlab. Theoretical and simulated results indicate that this model is mathematically strict and accurate, and its maximal precision can reach floating-point operations in the computer, with short time and fast speed. Thus it is very suitable to improve the optimal search speed and efficiency of local optimization methods based on the derivatives of merit function.It has outstanding performance in multilayer optical coating design with a large layer number.The optimality and speed of multilayer optical coating design has always been one of the focuses in the field of thin film optics. Dozens of optical coating optimization methods have been put forward with proved practical significance. And all these methods can be mainly distinguished into two types. One is called refinement or local method which needs a starting design and aims at speeding as fast as possible [1] . The other can be called synthesis or non-local method which can create a starting design by itself and mainly aims at enhancing optimal search capability [2][3][4][5][6][7] . In most of local methods based on merit function optimization, the accurate calculation and fast realization of merit function's gradient and Hesse matrix plays a key role in determining the actual speed of film design [8] . However, for the complexity of the calculation of merit function's Hesse matrix of multilayer coatings, many kinds of film design software only take advantage of gradient information to realize first order optimization. Researches about second order optimization method are seldom reported. And there is obvious speed difference between second and first order methods especially when a large layer number occurs. And then low speed and poor efficiency of home film design software is caused.In order to improve the computation efficiency, based on the matrix method to calculate the spectral properties of multilayer optical coating, an analytic model is established theoretically for absolutely accurate calculation of the merit function's gradient and Hesse matrix with respect to the design parameters. And the corresponding fast computation is successfully achieved by programming with Matlab, with the maximal precision of floating-point operations in the computer. The model is very suitable to improve the optimal search speed and efficiency of local optimization methods. And it has outstanding performance in multilayer optical coating design with a large layer number. The concept of modern film design is generally featured with constructing a merit function, with the structure and param...