So I have a DE which has no analytical solution and thus solved using NDSolve, involving some parameters like this:
I was hoping to use RegionPlot to illustrate the region of a,b which allows the solution to have always positive values over the region:
However, with some a,b, a singularity would be encountered and that also should be counted as not the correct region. But catching this is tricky. Also, the code above won't work anyway because somehow NDSolve would be evaluated before the values of a and b are inserted, effectively breaking everything. How would I get the correct result?