help-octave
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

ode23s, requested step-size too small


From: Torstein Fjermestad
Subject: ode23s, requested step-size too small
Date: Fri, 05 Apr 2013 18:54:29 +0200
User-agent: Roundcube Webmail/hz/2m6f.frs/uio-1

Dear all,

I am trying to solve a set of ordinary differential equations with ode23s.

In the output I get the following error message:



ode23s: requested step-size too small at t = 0, h = 4.65661e-10, err = -3.87203e-06 ode23s: requested step-size too small at t = 0.000770534, h = -0.000766662, err = 4.49164e-23 ode23s: requested step-size too small at t = 2.37561e-37, h = -4.2693e-37, err = -4.62411e-41 ode23s: requested step-size too small at t = 6.73057e-48, h = 1.57155e-50, err = 3.98936e-47
ode23s: requested step-size too small at t = 3.98936e-47, h =
error: called from 'ode23s' in file /home/torstein/octave/odepkg-0.8.4/ode23s.m near line 193, column 10
error: called from:
error: /usr/local/share/octave/3.6.4/m/miscellaneous/run.m at line 55, column 7


I am calling the ode23s function in the following way:


opt=odeset("InitialStep", 1e1,"RelTol", 1e-9, "AbsTol", 1e-12,"NormControl", "on","MaxStep",1e1);
[vt,y]= ode23s(f,[0 endTime],Z0,opt);

Is there a way to get rid of the error message?

The problem I am trying to solve is a microkinetic analysis on a network of elementary reactions. The rate constants are in the range from 10^-7 to 10^12 so the system is quite stiff.

Thanks in advance for your help.

Yours sincerely,

Torstein Fjermestad


reply via email to

[Prev in Thread] Current Thread [Next in Thread]