The Stack Overflow podcast is back! Listen to an interview with our new CEO.
2 of 6 added 1 character in body

Using shooting method to solve coupled ODE for domain walls

How can one solve a coupled equation for f[x] and a[x] with known boundary conditions for f[x] and a'[x] at some -xmax and xmax. The following approach doesn't work.

\[Kappa]=0.2;xmax=10;
sol = NDSolve[{1/\[Kappa]^2 f''[x] == (a[x]^2 - 1) f[x] + f[x]^3, a''[x] == 
a[x] f[x]^2, f[-xmax] == 0, f[xmax] == 1, a[-xmax] == 0, a[xmax] == 1}, {f,  
a}, {x, -xmax, xmax}, Method -> {"Shooting",  "StartingInitialConditions" ->  
{a'[-xmax] == 1/Sqrt[2], f'[-xmax] == 0}}]