Abstract. In this paper, we implement numerical solution of differential equations of frac- tional order based on hybrid functions consisting of block-pulse function and rationalized Haar functions. For this purpose, the properties of hybrid of rationalized Haar functions are presented. In addition, the operational matrix of the fractional integration is obtained and is utilized to convert compu...