Added odesolver?
Uses odesolver, currently wrong initial eta(1)!!! Must rewrite three functions