In models.R, moved (1-h) from

A%*%N1*(1-h)
to
build_Leslie(s = s*(1-h), f = f)

Exact same results, but easier to understand