restart; # Sub in the first two terms of the w1 solution near the middle
# singularity wold := c0 + c0*(1+6*R^2)/10*eta^2 + w(eta)*eta^3; NiM+SSV3b2xkRzYiLChJI2MwR0YlIiIiKihGJ0YoLCZGKEYoKiRJIlJHRiUiIiMiIidGKEkkZXRhR0YlRi0jRigiIzUqJi1JIndHRiU2I0YvRihGLyIiJEYo -eta*(1-R^2*eta^2)*diff(wold,eta,eta) + 2*diff(wold,eta) + eta*wold; NiMsLCooSSRldGFHNiIiIiIsJkYnRicqJkkiUkdGJiIiI0YlRishIiJGJywqKiZJI2MwR0YmRicsJkYnRicqJEYqRisiIidGJyNGJyIiJiomLUklZGlmZkdJKnByb3RlY3RlZEdGODYkLUkid0dGJjYjRiUtSSIkR0Y4NiRGJUYrRidGJSIiJEYnKiYtRjc2JEY6RiVGJ0YlRitGMiomRjpGJ0YlRidGMkYnRiwqKEYvRidGMEYnRiVGJyNGK0Y0KiZGQkYnRiVGQEYrKiZGOkYnRiVGK0YyKiZGJUYnLChGL0YnKihGL0YnRjBGJ0YlRisjRiciIzUqJkY6RidGJUZARidGJ0Yn collect(%, diff(w(eta),eta)); NiMsLComLCYqJEkkZXRhRzYiIiIkIiIjKiZGJ0YpLCYiIiJGLSomSSJSR0YoRipGJ0YqISIiRi0hIidGLS1JJWRpZmZHSSpwcm90ZWN0ZWRHRjQ2JC1JIndHRig2I0YnRidGLUYtKihGJ0YtRixGLSwoKiZJI2MwR0YoRi0sJkYtRi0qJEYvRioiIidGLSNGLSIiJiomLUYzNiRGNi1JIiRHRjQ2JEYnRipGLUYnRilGLSomRjZGLUYnRi1GP0YtRjAqJkYnRi0sKEY8Ri0qKEY8Ri1GPUYtRidGKiNGLSIjNSomRjZGLUYnRilGLUYtRi0qKEY8Ri1GPUYtRidGLSNGKkZBKiZGNkYtRidGKkY/