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/60e6169f-ccc0-45b7-81c3-e09ab02b8a24%40googlegroups.com.

Reply via email to