Operational matrix approach for solution of integro-differential equations arising in theory of anomalous relaxation processes in vicinity of singular point