Unpacking a Mathematica expression returned by DSolve

Question

I was trying to solve a system of differential equations in Mathematica and had troubles understanding what the solution looked like. So I wanted help to unpack it.

I had a system of two coupled differential equations y1[x] , y2[x] and the solution Mathematica spits out is the following

{{y2 -> Function[{x}, 
    C[1] InverseFunction[
       Inactive[Integrate][
         1/((Conjugate[K[1]] - C[1] Conjugate[C[1] K[1]]) K[1]^2), {K[
           1], 1, #1}] &][-((I x α)/ℏ) + C[2]]], 
  y1 -> Function[{x}, 
    InverseFunction[
      Inactive[Integrate][
        1/((Conjugate[K[1]] - C[1] Conjugate[C[1] K[1]]) K[1]^2), {K[
          1], 1, #1}] &][-((I x α)/ℏ) + C[2]]]}}

I do not understand function, inverse function, and inactive integrate methods. I was wondering how can I understand the solutions qualitatively and quantitatively.

EDIT : The original code about the two differential equations :

eq1 = I*ℏ*
   D[y1[x], x] - α*(Conjugate[y1[x]]*y1[x] - Conjugate[y2[x]]*y2[x])*
   y1[x]
eq2 = I*ℏ*
   D[y2[x], x] - α*(Conjugate[y1[x]]*y1[x] - Conjugate[y2[x]]*y2[x])*
   y2[x] 
eqns = {eq1 == 0, eq2 == 0};

Answer 1

To answer the question about the form of the solution:

In solving DEs, DSolve often has to integrate (obviously) an equation. Sometimes the integral cannot be solved. So it ends up with an equation $F ( x , y ) = 0$ in terms of an integral that is inactivated because it couldn’t be solved. Next DSolve tries to solve for $y$. If the equation is of the form $x = f ( y )$, we get $y = f^{ − 1} ( x )$. If $f^{ − 1}$ cannot be computed (usually can’t when it involves an inactive integral), we get solution for $y$ in terms of an inverse function.

Here's a way to unwind the steps to get the solution in implicit form, still in terms of an unevaluated integral. We'll transform the expressions with three or four replacement rules. The last rule is optional.

(* Convert y -> Function[{x}, a] to y[x] == a *)
{HoldPattern[dv_ -> Verbatim[Function][{iv__}, body_]] :> dv[iv] == body} 

(* Solve each equation for the InverseFunction[..][..] in it *)
{eq_Equal :> First@ Solve[eq, Cases[eq, _InverseFunction[_], Infinity]]}

(* Rewrite  $f^{-1}(a) \rightarrow z$  as  $f(z) = a$  *)
{{Verbatim[InverseFunction][f_][a_] -> z_} :> f[z] == a}

(* Change dummy variable K[1] to something more traditional *)
{K[1] -> ξ}

Apply the rules in succession with ReplaceAll (/.):

sol /.
    {HoldPattern[dv_ -> Verbatim[Function][{iv__}, body_]] :> 
      dv[iv] == body} /.
   {eq_Equal :> 
     First@ Solve[eq, Cases[eq, _InverseFunction[_], Infinity]]} /.
  {{Verbatim[InverseFunction][f_][a_] -> z_} :>
     f[z] == a} /.
 {K[1] -> ξ}

Answer 2

Trying to solve a Cauchy problem by

ClearAll["Global`*"];
eq1 = I*\[HBar]*D[y1[x],x] - \[Alpha]*(Conjugate[y1[x]]*y1[x] - Conjugate[y2[x]]*y2[x])*y1[x];
eq2 = I*\[HBar]* D[y2[x],x] - \[Alpha]*(Conjugate[y1[x]]*y1[x] - Conjugate[y2[x]]*y2[x])* y2[x];
eqns = {eq1 == 0, eq2 == 0, y1[1] == I, y2[1] == 2};
DSolve[eqns, {y1, y2}, x]

, I obtain {} and a warning "DSolve::bvnul: "For some branches of the general solution, the given boundary conditions lead to an empty solution."" The same with other initial values.

The ParametricNDSolve command succeeds:

\[HBar] = 1; \[Alpha] = 2;
eqns = {eq1 == 0, eq2 == 0, y1[1] == a, y2[1] == b};
sol = ParametricNDSolve[eqns, {y1, y2}, {x, 1, 2}, {a, b}];
y1 = y1[I, 2] /. sol;
y1[1.5] /. sol

-0.14112 - 0.989992 I