I have asked about precision and accuracy settings in the past and my most recent post got no replies. So, I am going to ask the question differently this time.
I have a function that I need to numerically integrate. It is a function with two parametric parameters. In performing the integration varying the parameters two things happen.
At some threshold of one parameter, the result seems to be similar to a bifurcation. The results abruptly take a step to a new value and as I narrow the domain to "zero in" on the jump point, it does not go away even as I increase the precision and accuracy requirements to extreme levels until it simply never runs to completion. I have never worked with functions that exhibited a bifurcation. How does Mathmatica deal with these, which leads to the next question?
As I increase the accuracy and precision varying the second parameter, the Plot goes from smooth to broken up. As the accuracy and precision is increased further and using the Table and ListPlot functions, the calculated points begin to look like data with a noise component. In other words, the plotted points look as if there is some statistical component and the occasional outlier seems to be consistent with the RMS limits given the vertical scale used in the plot. Also, with enough plotted points, the "smooth" curve from the less precise calculations seems to be a regression fit to the higher precision data points.
I am puzzled by what I may be looking at. It seems as though my function is "funky" to begin with but the transitions in the way the integral is numerically calculated may be introducing unwanted noise into the results. If these descriptions make sense to anyone, I would appreciate some discussion.