In recent years, the concept of fractal analysis is the best nonlinear tool towards understanding the complexities in nature. Especially, fractal interpolation has flexibility for approximation of nonlinear data obtained from the engineering and scientific experiments. Random fractals and attractors of some iterated function systems are more appropriate examples of the continuous everywhere and nowhere differentiable (highly irregular) functions, hence fractional calculus is a mathematical operator which best suits for analyzing such a function. The present study deals the existence of fractal interpolation function (FIF) for a sequence of data { (xn, yn) : n≥ 2 } with countable iterated function system, where xn is a monotone and bounded sequence, yn is a bounded sequence. The integer order integral of FIF for sequence of data is revealed if the value of the integral is known at the initial endpoint or final endpoint. Besides, Riemann–Liouville fractional calculus of fractal interpolation function had been investigated with numerical examples for analyzing the results. © 2016, Springer International Publishing.