I'm trying to solve a set of differential equations that all depend on a parameter, $\kappa$. I can use the system of ODEs to reduce the four equations into one second order differential equation, for $y[x;\kappa]$. I've seen certain tricks to solving equations such as $d^2y/dx^2 = \kappa y[x]$.But I can not put my equation in that form. So the real trick is,how do I solve for $\kappa$ and $y[x]$ with only one equation.
My idea is as follows, I will look at the limit of the differential equation at $x =0$ (ODE1), and the limit of it at the other boundary $x = a$, (ODE2).Then I will technically have two differential equations. Then I pick a value for $\kappa$, from some physics that I know about the problem, i.e. set $\kappa = constant$. After this I solve ODE1 in the domain $ x= 0$ to $x0$ and I solve ODE2 in the domain $x0$ to $x=a$. Finally I check whether or not the function $y[x]$ matches at the point $x0$. If it doesn't match I try this again until I find the value of $\kappa$ that makes the two sides match up.
Its the last part that I don't know how to implement? How do I use the mis-match to get me a new value of $\kappa$? I've read several sources about this idea, under "shooting method", but I really don't get how I should update the values of $\kappa$.
Thanks for reading, and any help is much appreciated.
1st Edit Hey Marty Cohen,
Thanks for the response. I was thinking about using the secant method, but my confusion lies in this: The nth root for the secant method is given by
$x_n = x_{n-1} - f(x_{n-1})* \frac{x_{n-2} - x_{n-1}}{f(x_{n-2}) - f(x_{n-1})}$,
so for me this boils down to
$\kappa_3 = \kappa_{2} - f(\kappa_{2})* \frac{\kappa_{2} - \kappa_{1}}{f(\kappa_{2}) - f(\kappa_{1})}$,
and I'm ok with this so far. I just need to pick two values of $\kappa$ to get the third. But what I don't get is this function $f(\kappa$). Even if I integrate both ODE1 and ODE2 and check the mismatch at a fixed point ($x0$), so something like
$\kappa_3 = \kappa_{2} - f(x0,\kappa_{2})* \frac{\kappa_{2} - \kappa_{1}}{f(x0,\kappa_{2}) - f(x0,\kappa_{1})}$,
the question is : Which f do I use, is it $f1$ which is the solution of ODE1 or $f2$ which is the solution of ODE2. It was because of the mistmatch in functions that I thought I would need to use something like Newton-Raphson where I would have to build derivatives.
Sorry it took me so long to get to the question. Also this problem is related to bit of research I'm doing, so no worries about what my prof wants. I'm just doing this for me, and I really need to understand it, in order to move the research forward.
Thanks again.