Hello
I demand your help .
I'm going to study the accretion-ejection disks using the DVODE-f90 code.
when I compile my example. the calcul stops at T=3,21.
the above message, R1 = 0.3215402479903D+01
At t = 0.3215D+01 y = 0.321568D+01 0.158443D+14 0.193170D+39
Error halt: ISTATE = -1
Can you help me how to integrate my equations system for large T (INFINITY)
you find in the attached files my source file and the dvode_f90 subroutine.
I demand your help .
I'm going to study the accretion-ejection disks using the DVODE-f90 code.
when I compile my example. the calcul stops at T=3,21.
the above message, R1 = 0.3215402479903D+01
At t = 0.3215D+01 y = 0.321568D+01 0.158443D+14 0.193170D+39
Error halt: ISTATE = -1
Can you help me how to integrate my equations system for large T (INFINITY)
you find in the attached files my source file and the dvode_f90 subroutine.