[differentialequations.jl] stuck at initial value when using function

You need du .= A*u so that du is modified instead of replaced.