Coupled non-linear PDEs using DifferentialEquations.jl and the MOL

Thx for the reply.

I see what you mean, with an issue similar to this topic: Best practice for avoiding line break bugs in long expressions

So I’ve put everything on 1 line again (which is not great for readibility but anyway) and… It doesn’t change anything :sweat_smile:

I don’t know why, but julia seems to understand my notation if it is in a loop, like this work:

test = zeros(2)
test2 = zeros(2)

for i in 1:2
    test[i] = 
        1 +
        2
    test2[i] = 1 + 2
    println("test[i]", test[i])
    println("test2[i]", test2[i])
end

shows:

test[i]3.0
test2[i]3.0
test[i]3.0
test2[i]3.0

So anyway, I am back to my problem.