summaryrefslogtreecommitdiff
path: root/Master/texmf/asymptote/ode.asy
blob: 5135dd9443fa0fbbea36d2923166237e186fc9a3 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
real euler(real y, real f(real x, real y), real a, real b=a, int n=0,
           real h=0, bool dynamic=false, real tolmin=0, real tolmax=0)
{
  if(h == 0) {
    if(b == a) return y;
    if(n == 0) abort("Either n or h must be specified");
    else h=(b-a)/n;
  }
  real x=a;
  for(int i=0; i < n; ++i) {
    y += h*f(x,y);
    x += h;
  }
  return y;
}