I have system of first order ordinary diff equations, indipendent variable is x cordinate. I know asymptotic solution in left and right region (i.e. when x->-infinity or x->infinity, e.g. when abs(x)>1000), it's const plus exponentially falling function. I want to find numerical solution in middle region, witch will have given asymptotic in left and right region. If I give initial value at left (i.e. at x=-1000) numerical solution blow up at right (I also have exponentially growing functions on right). How to do this?