A numerical method for solution of integro-differential equations of fractional order