Hi, I am having trouble bypassing or including a singularity for a generalized one dimensional flow model. Here is the ODE I am solving for, but as expected there is a singularity where the Mach number=1. I am relatively new to Mathematica, so I am struggling with the syntax a bit. Anyway, here is the ODE I am solving for [M]'[x] == ( 70 M[x] (1 + 0.125 M[x]^2) (1 + 1.25 M[x]^2))/((200 + 140 x) (1 - M[x]^2)) for x=0 to 20 Is there a way to exclude the solution for M[x]=1 or there is another form of the equation that can handle the singularity, but it is only valid for M=1, so is there some way to use that equation for M=1 and combine it with the original ODE?