octave:1> looptest2 Options for LSODE include: keyword value ------- ----- absolute tolerance 1.49012e-08 relative tolerance 1.49012e-08 integration method stiff initial step size -1 maximum order -1 maximum step size -1 minimum step size 0 step limit 100000 P0 = 5.0000e+16 T0 = 1000000 ne0 = 5.0000e+10 EH = 0.0020000 Npoints = 1000 s1 = 10000 sN = 5.0000e+09 Tchromosphere = 3.1623e+04 gradTchromosphere = 0 FCchromosphere = 0