Fixed error in original numerical solver function where I had incorrectly calculated the competition effects on species j in site q. It is slower but now is correct.