I am trying to solve a set of coupled ODEs where the solution variable is a complex array. When I run the following code, I get an error which says:

**LoadError: InexactError: Int64(Int64, 1731//50000)**

The Stacktrace sarts from

**[1] convert(::Type{Int64}, ::Rational{Int64}) at ./rational.jl:85
[2] OrdinaryDiffEq.Vern9ConstantCache(::Type, ::Type) at /Users/home/.julia/packages/OrdinaryDiffEq/YudGf/src/tableaus/verner_tableaus.jl:2539**

I have not used Int64 or Rational anywhere, and this error does not occur when I try to solve, e.g., the Lorenz equations. So I think this might be a problem with how OrdinaryDiffEq handles complex numbers.

All my real numbers are Float64, and the real and imaginary parts of my Complex numbers are also Float64.

I have also confirmed that it doesn’t matter what calculations I actually place inside the function

**VortexStreet**which is being solved by

**solve**; even if I do no calculations inside this function, the same error occurs.

What could be the problem?

```
using DifferentialEquations
function VortexStreet(dzdt,z,p,t)
for i = 1:length(z)
for j = 1:length(z)
if j != i
dzdt[i] += cot(z[i]-z[j]);
end
end
end
end
a = 1.0; b = 0.56; L = 2a;
n = 3;
z0 = Array{Complex}(undef,2n);
dzdt = Array{Complex}(undef,2n);
for i in 1:n
z0[2(i-1)+1] =complex(0.5*a+(i-1)*L,b*0.5);
z0[2(i-1)+2] =complex(1.5*a+(i-1)*L,b*0.5);
end
for i in 1:2n
dzdt[i] = 0.0;
end
tspan = (0.0,0.1);
problem = ODEProblem(VortexStreet,z0,tspan);
solution = solve(problem);
```