Crank-Nicholson method for rate equations in powder random lasers. (21st January 2015)