On the box dimension of recurrent fractal interpolation functions defined with Matkowski contractions