- square of a negative was reason for NANs

- plot PopSpring(t+1)~popFall