Results 1 to 2 of 2

Math Help - Numerical solution of a first order ODE

  1. #1
    Newbie
    Joined
    Apr 2010
    Posts
    1

    Numerical solution of a first order ODE

    Hello,

    My question seems to be simple. I would like to numerically solve the following first order ODE to obtain v(x):

    v'(x) = b.[v(x) - f(x)] , given boundary condition v(+infinity) = 0 [b is a known constant]. OR equivalently realize operator b/(b-D) , where D is d/dx

    These are the problems:
    1)) f(x) is not known explicitly (f(x) is sampled at a dense x grid).
    2)) f(x) has a discontinuity, i.e. there is a large jump in df(x)/dx
    3)) This is a [final value problem], so I need an algorithm that works backwards instead of the usual forward Runge-Kutta algorithms

    I know that using the integrating factor, we obtain the following solution
    v(x) = b.exp(b.x) . integral(from -infinity to x) [exp(-b.y).f(y) dy ].
    However, given we have a dense x grid, this means I need to evaluate the above integral many times, which is computationally very slow.


    Can you please help me? any hints for solving the ODE or realize the (b/(b-D)) operator will be highly appreciated...
    Follow Math Help Forum on Facebook and Google+

  2. #2
    Grand Panjandrum
    Joined
    Nov 2005
    From
    someplace
    Posts
    14,972
    Thanks
    4
    Quote Originally Posted by tamzam600 View Post
    Hello,

    My question seems to be simple. I would like to numerically solve the following first order ODE to obtain v(x):

    v'(x) = b.[v(x) - f(x)] , given boundary condition v(+infinity) = 0 [b is a known constant]. OR equivalently realize operator b/(b-D) , where D is d/dx

    These are the problems:
    1)) f(x) is not known explicitly (f(x) is sampled at a dense x grid).
    2)) f(x) has a discontinuity, i.e. there is a large jump in df(x)/dx
    3)) This is a [final value problem], so I need an algorithm that works backwards instead of the usual forward Runge-Kutta algorithms

    I know that using the integrating factor, we obtain the following solution
    v(x) = b.exp(b.x) . integral(from -infinity to x) [exp(-b.y).f(y) dy ].
    However, given we have a dense x grid, this means I need to evaluate the above integral many times, which is computationally very slow.


    Can you please help me? any hints for solving the ODE or realize the (b/(b-D)) operator will be highly appreciated...
    I would consider changing the variable to z=1/x then you will have an initial value problem.

    CB
    Follow Math Help Forum on Facebook and Google+

Similar Math Help Forum Discussions

  1. Numerical solution of a (system of?) DE
    Posted in the Advanced Math Topics Forum
    Replies: 3
    Last Post: May 27th 2010, 02:15 PM
  2. Numerical solution in mathematica
    Posted in the Math Software Forum
    Replies: 0
    Last Post: May 12th 2010, 10:22 AM
  3. Help with numerical solution of PDE
    Posted in the Differential Equations Forum
    Replies: 3
    Last Post: March 20th 2010, 01:37 PM
  4. numerical solution
    Posted in the Calculus Forum
    Replies: 3
    Last Post: May 7th 2008, 04:30 AM
  5. Numerical solution w/o derivative
    Posted in the Calculus Forum
    Replies: 3
    Last Post: October 31st 2007, 04:47 PM

Search Tags


/mathhelpforum @mathhelpforum