Mathematica에서 shooting method를 이용해서 비선형 미분방정식의 해를 구하자.
f′=1x(1−a)fa′=−x2f2(f2−1)
경계조건은 f(0)=0, f(∞)→0, a(0)=0을 만족시켜야 한다. x가 증가하면 f(x)는 빠르게 1로 수렴하므로 오른쪽 경계조건은 f(15)→1로 잡아도 충분하다. x=0에서의 apparent singularity를 Mathematica가 처리할 수 있도록 방정식을 변형시켜주어야 한다.

728x90
'Physics > Research' 카테고리의 다른 글
Ricci tensor 계산 (0) | 2012.12.21 |
---|---|
Shooting method for the boundary value problems: Maple (0) | 2012.08.01 |
Relaxation Methods for Differential Equations (0) | 2012.07.26 |
Wald’s Entropy, Area & Entanglement (0) | 2012.02.25 |
Topics in f(R) THEORIES OF GRAVITY (0) | 2012.02.24 |