2Â¥: Originally posted by xzczd at 2015-10-14 09:50:32
NDSolve»áû·´Ó¦°Ë³ÉÊÇÄãÓï·¨´íÁË¡£ÊÔÊÔÕâ¸ö£º

NDSolve
Plot[Abs /. nsol // Evaluate, {t, 0, 10}]



ÖÁÓÚDSolve£¬½ØÖ¹10.2°æ£¬MathematicaµÄDSolveÊÇÓÐÐí¶à²»¾¡ÈçÈËÒâµÄµØ·½µÄ¡££¨²»¹ý¼´½«ÍƳöµÄ10.3°æËÆ ...

NDSolve[{I x'[t] == 3/10 x[t] - 4/10 Exp[I t/2] y[t],
  I y'[t] == -4/10 Exp[-I t/2] x[t] - 3/10 y[t], x@0 == 0, y@0 == 1}, {x, y}, {t, 0, 10}]
Plot[Abs[{x[t], y[t]}] /. nsol // Evaluate, {t, 0, 10}]

eqn = {I x'[t] == 3/10 x[t] - 4/10 Exp[I t/2] y[t],
   I y'[t] == -4/10 Exp[-I t/2] x[t] - 3/10 y[t], x@0 == 0, y@0 == 1};
representWithy = Solve[eqn[[2]], x[t]][[1, 1, -1]];
eqn[[{1, 3, 4}]] /. x -> (Function[t, #] &@representWithy)
ysol = DSolve[{%, eqn[[-1]]}, y, t];
(*y[t]=*)y[t] /. ysol
(*x[t]=*)representWithy /. ysol

2Â¥: Originally posted by xzczd at 2015-10-14 09:50:32
NDSolve»áû·´Ó¦°Ë³ÉÊÇÄãÓï·¨´íÁË¡£ÊÔÊÔÕâ¸ö£º

NDSolve
Plot[Abs /. nsol // Evaluate, {t, 0, 10}]



ÖÁÓÚDSolve£¬½ØÖ¹10.2°æ£¬MathematicaµÄDSolveÊÇÓÐÐí¶à²»¾¡ÈçÈËÒâµÄµØ·½µÄ¡££¨²»¹ý¼´½«ÍƳöµÄ10.3°æËÆ ...

3Â¥: Originally posted by xzczd at 2015-10-14 09:52:54
°¡£¬ysol = DSolve[{%,eqnÕâÐеÄeqn[]ÊǶàÓàµÄ£¬²»¹ýÕâ¸ö²»Ó°Ïì½á¹û¡£...