Processing math: 100%

Mathematica에서 shooting method를 이용해서 비선형 미분방정식의 해를 구하자. 

f=1x(1a)fa=x2f2(f21)

경계조건은 f(0)=0, f()0, a(0)=0을 만족시켜야 한다. x가 증가하면 f(x)는 빠르게 1로 수렴하므로 오른쪽 경계조건은 f(15)1로 잡아도 충분하다. x=0에서의 apparent singularity를 Mathematica가 처리할 수 있도록 방정식을 변형시켜주어야 한다.

 

 

 

728x90
,