  We consider a series of massive scaling limits m_1 -> infty, q -> 0, lim m_1
q = Lambda_{3} followed by m_4 -> infty, Lambda_{3} -> 0, lim m_4 Lambda_{3} =
(Lambda_2)^2 of the beta-deformed matrix model of Selberg type (N_c=2, N_f=4)
which reduce the number of flavours to N_f=3 and subsequently to N_f=2. This
keeps the other parameters of the model finite, which include n=N_L and
N=n+N_R, namely, the size of the matrix and the "filling fraction". Exploiting
the method developed before, we generate instanton expansion with finite g_s,
epsilon_{1,2} to check the Nekrasov coefficients (N_f =3,2 cases) to the lowest
order. The limiting expressions provide integral representation of irregular
conformal blocks which contains a 2d operator lim frac{1}{C(q)} : e^{(1/2)
\alpha_1 \phi(0)}: (int_0^q dz : e^{b_E phi(z)}:)^n : e^{(1/2) alpha_2 phi(q)}:
and is subsequently analytically continued.
