In this work, we proposed an effective method based on cubic and pantic B-spline scaling functions to solve partial differential equations of fractional order. Our method is based on dual functions of B-spline scaling functions. We derived the operational matrix of fractional integration of cubic and pantic B-spline scaling functions and used them to transform the mentioned equations to a syste...