Hey everyone I have a paper that has a method for computing rise and set times of a satellite given a closed form solution. It is a complicated sinusoidal function and the paper has a method to calculate a single rise set over a single period using newton rhapson method. I want to find a way to calculate all zeros over a given interval of time.
Here is the equation
Here is a picture of one of the initial intervals of the plot
I now newton approach will work given a good inital guess. My first thought was to create a loop and when it finds a zero to move forward to another set, but I find it hard to believe there is not a method already developed.
I have come across ode event detection in MATLAB, or a lot of other solutions that seem to be for polynomials but wanted to see if there were better approaches a
Here is a link to the paper https://arc.aiaa.org/doi/abs/10.2514/3.2057