Hello, I have numerically solved the system below: dalpha/dx = -beta/(c_1) dbeta/dx = -(c_2)*beta/(c_1)+(c_3)*gamma*dgamma/dx dgamma/dx = (c_4)*[alpha^(c_5)]/c_2 where c_1, c_2, etc. are specified constants. If one plots beta as a function of x, there's a peak. I would like to use some approximations (perhaps linearization?) to get the distance between a region where dbeta/dx is approximately zero and where there is a peak in beta, after given an initial value of alpha, beta, and gamma where dbeta/dx ~ 0. Any suggestions?