We present the results of analysis and modelling of the eclipsing binary system, KIC 10661783. The Fourier analysis of the Kepler light curve, corrected for the binary effects, reveals 750 frequency peaks in both p and g-mode regions. Those with the highest amplitudes concentrate in the range of 20-30 d$^{-1}$. To reproduce observed spectrum of frequencies we construct seismic models accounting for the mode instability. In order to obtain instabilities in the g-mode regime we modify the opacity tables data near the Z-bump. In order to reproduce system parameters, we construct evolutionary models including binary evolution.