@Preben Alsholm 

You are right (i'm capitan obvious)! I use 

dsolve({ic_Im, ic_Re, `Eqn_ψIm`, `Eqn_ψRe`}, numeric, method = bvp[middefer], abserr = .1, maxmesh = 100000)

with condition in 1000

and i have result. I don't know right or wrong, but i have it!

Thank you!

@Preben Alsholm 

I tried to use psi in 10, for example. But nothing changes.