Stop at singularity with lsode

[jchtt]

Hi,

I am trying to solve a system of two coupled differential equations with 
Octave, namely those for Hathaway's circular pursuit curve and I am running 
into the problem that there of course is a singularity in case of a 
catch-event. However, I don't know when this is going to happen, so I can't 
just restrict the time frame before issuing the lsode command. Is there any way 
to tell lsode to stop at a singularity and/or at least produce output up to the 
point where the step limit is exceeded?

Thanks in advance,

jchtt