Thank you. It's an easy problem to reproduce. You can cut-paste that code into the official online Sage Cell Server
https://sagecell.sagemath.org/ and see for yourself what I meant. Best regards Toan On Thursday, April 16, 2020 at 10:17:47 PM UTC+7, vdelecroix wrote: > > Hi, > > What do you mean by crash? Could you copy paste the *full* traceback > that you obtain? > > Le 16/04/2020 à 17:04, Toan T Nguyen a écrit : > > Hi all, > > > > I'm trying to solve a system of 1st order ODE using desolve_tides_mpfr() > > arbitrary precision module. My script crashes at the external gcc call. > > So I try this very simple example I found from a tutorial in a website. > It > > also crashes at the same gcc call. I test it in the SageMath online cell > > server, it crashes as well, so it's not just my code or my setup. > > Could somebody help with this example? Is this a bug, or the script is > > wrong? > > > > Thank you > > > > Toan > > > > var('t,x,y,z') > > s = 10 > > r = 28 > > b = 8/3 > > f(t,x,y,z)= [s*(y-x),x*(r-z)-y,x*y-b*z] > > x0 = > > > -13.7636106821342005250144010543616538641008648540923684535378642921202827747268115852940239346395038284 > > > > y0 = > > > -19.5787519424517955388380414460095588661142400534276438649791334295426354746147526415973165506704676171 > > > > z0 = 27 > > T = > > > 15.586522107161747275678702092126960705284805489972439358895215783190198756258880854355851082660142374 > > > > sol = desolve_tides_mpfr(f, [x0, y0, z0],0 , T, T, 1e-100, 1e-100, 100) > # > > optional - tides > > print(sol) > > > -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-support+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/sage-support/ba0eed69-0a08-4616-a27f-0c4a085f91ca%40googlegroups.com.
