边值问题的解
psol = ParametricNDSolveValue[{y''[t] + Sin[y[t]] == 0, y[0] == 0, y'[0] == s}, y, {t, 0, 10}, {s}];Plot[Evaluate@Table[psol[s][t], {s, 0, 2, 0.05}], {t, 0, 10}, ImageSize -> Medium]pfun = ParametricNDSolveValue[{y''[t] + Sin[y[t]] == 0, y[0] == 0, y'[0] == s}, y[10], {t, 0, 10}, {s}];
val = Table[FindRoot[pfun[s], {s, s0}], {s0, {0.9, 1.9, 2}}]Show[Plot[Evaluate@Table[psol[a][t], {a, 0, 2, .05}], {t, 0, 10}, PlotRange -> All, PlotPoints -> 500, PlotStyle -> Lighter[Gray, 0.5]], Plot[Evaluate[psol[s][t] /. val], {t, 0, 10}, PlotStyle -> Thick], ImageSize -> Medium]