Hi Greg,
Thanks for your reply.
I think my issues start when I try to split the big vector v, into 4 small vectors (one for each equation:v1,v2,v3,v4).
After that, the boundary conditions(vv1,vv2,vv3,vv4) is not correctly implemented (it's need to be created dynamically).
In conclusion, there...