rk works, new is to test different structure with stepper
fixed theta<0 at n=3/2.