Added previous lnp code.  TODO- make stimulus an argument, make nonlinearity a parameterized sigmoid, and decide how to fix firing rate