[ alpha and R ] - [ double vs quad ] - [ Umax(delta) ] - [ Umax(cte) ] |

[top]

Here we plot alpha(t) and R(t) (for Tc=0.5):

[top]

ubar (approx inside the core) in red

psi (approx outside of the core) in blue

The largest curves (resp. double and quad) corresponds to highest possible u0 obtained numerically.

[top]

M=1+|u|^{0.5+delta}

The numerics stopped when:

oscillations: when umax(t) reached a maximum and then started oscillating

crash: when the numerics crashed (stiffness/min step reached)

[top]

M=cte+|u|^{0.5}

The numerics stopped when:

oscillations: when umax(t) reached a maximum and then started oscillating

crash: when the numerics crashed (stiffness/min step reached)