Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- For[k = 5, k <= 7, k = k + 0.1,
- sol =
- ParametricNDSolve[
- {Y''[x] + k^2 (1 - 0.1 Conjugate[Y[x]] Y[x]) Y[x] == 0,
- Y[1] == N E^(I k), Y'[1] == I k N E^(I k)},
- Y[x], {x, 0, 1}, {N},
- AccuracyGoal -> Infinity, MaxSteps -> Infinity];
- Ns = FindRoot[(Y'[N][0] + I k Y[N][0] == 2 I k) /. sol, {N, 0}];
- Print[Ns]]
- sol = ParametricNDSolveValue[
- {
- Y''[x]+k^2 (1-0.1 Conjugate[Y[x]] Y[x]) Y[x]==0,
- Y[1] == n E^(I k),
- Y'[1] == I k n E^(I k)
- },
- {Y[0], Y'[0]},
- {x, 0, 1},
- {k, n}
- ];
- ns[k_] := n /. FindRoot[sol[k, n] . {k I, 1} == 2 k I, {n, .5+.5I}]
- Table[ns[k], {k, 5, 7, .1}]
Add Comment
Please, Sign In to add comment