Results 1 to 2 of 2

Thread: stiff ode solving

  1. #1
    Sep 2009

    Unhappy stiff ode solving

    Dear friends,

    My problem is a kind of stiff ordinary differential equation of 11 first order equations (state space equations). I have tried to solve it by ode15s and ode23s but there are some situations that the integrator is not successful and fails to integrate, the message which is displayed is shown below please let me know how to solve the problem.

    Warning: Failure at t=2.694000e+002.
    Unable to meet integration tolerances
    without reducing the step size below the
    smallest value allowed (9.094947e-013) at
    time t.
    > In ode15s at 753
    Follow Math Help Forum on Facebook and Google+

  2. #2
    Super Member
    Aug 2008
    Quote Originally Posted by hadidehghanzade View Post
    please let me know how to solve the problem.
    Take for example, the equation:

    $\displaystyle y'=\frac{3x^2}{3y^2-4},\quad y(1)=0$

    If I try to numerically integrate that I get:

    Clear[x, y]
    NDSolve[{Derivative[1][y][x] == 
        (3*x^2)/(3*y[x]^2 - 4), y[1] == 0}, 
      y, {x, 1, 2}]
    During evaluation of In[105]:=
     NDSolve::ndsz:At x == 1.597809453897525`, 
    step size is effectively zero;
     singularity or stiff system suspected. >>
    {{y -> InterpolatingFunction[]}}
    The (full) solution is below. Why is the integrator giving me that message?
    Attached Thumbnails Attached Thumbnails stiff ode solving-bending-reed.jpg  
    Last edited by shawsend; Sep 25th 2009 at 06:50 AM.
    Follow Math Help Forum on Facebook and Google+

Similar Math Help Forum Discussions

  1. Are these ODE's stiff?
    Posted in the Differential Equations Forum
    Replies: 1
    Last Post: Nov 11th 2009, 02:36 AM
  2. Solving x'Px = v
    Posted in the Advanced Algebra Forum
    Replies: 0
    Last Post: Dec 11th 2008, 03:21 PM
  3. stiff differential equations [prove]
    Posted in the Calculus Forum
    Replies: 3
    Last Post: Mar 31st 2008, 08:58 AM
  4. Solving
    Posted in the Algebra Forum
    Replies: 1
    Last Post: Oct 3rd 2007, 08:28 AM
  5. Replies: 3
    Last Post: Oct 11th 2006, 09:15 PM

Search Tags

/mathhelpforum @mathhelpforum