Numerical solution in Mathematica

  • Mathematica
  • Thread starter FanOfGR
  • Start date
  • #1
2
0
Hello,

I have a problem with numerical solution of the system of ODE's (geodesic eqs.) in Mathematica. The only relevant command I have found is NDSolve. It works perfectly, unless my equations contain functions which are again computed using numerical methods.

Simple example illustrating the problem.
Let's define function

F[x_] := Module[{}, FindRoot[ x y == 1, {y, 1/x} ][[1,2]] ];

I don't want to write F[x_] = 1/x, because I want Mathematica to compute it numerically. This function gives, of course, correct values, e.g. F[2]=0.5.

Now I want to solve equation y' = F with initial condition y(1)=0 using

sol = NDSolve[ {y'[x]== F[x], y[1]==0}, y, {x, 1, 3} ];

If I define F[x_] = 1/x, it works. If I define it by FindRoot, it doesn't. It seems to me that NDSolve tries to work with F[x] in symbolic form, but then FindRoot says, that 1/x is not a number.

I hope that I described the problem clearly. Please, if anybody knows the solution, tell me, I am getting crazy :)

Thank you
 

Answers and Replies

  • #2
64
0
The main problem is that, Mathematica tries to solve your problem analytically first. So it plugs in symbolic x, and your function can't handle it.

Please refer to this question in http://www.voofie.com/concept/Mathematica/" [Broken]:

http://www.voofie.com/content/163/ndsolve-in-mathematica-with-numerical-function/" [Broken]

For ways to solve your problem.
 
Last edited by a moderator:

Related Threads on Numerical solution in Mathematica

Replies
3
Views
2K
Replies
5
Views
2K
Replies
8
Views
2K
Replies
6
Views
10K
Replies
4
Views
2K
  • Last Post
Replies
3
Views
594
Replies
1
Views
2K
Top