demo/romberg.dem

   1 define_variable(pi,apply('ev,[%pi,numer:true]),float)$
   2
   3 define_variable(freq,3.3,float,"the frequency in radians/second")$
   4
   5 v(x):=(modedeclare(x,float),x*abs(sin(freq*x)))$
   6
   7 h1(freq):=romberg('v,0,pi)$
   8
   9
  10 /* an alternative syntax is */
  11
  12 h2(freq):=romberg(t*abs(sin(freq*t)),t,0,pi)$
  13
  14
  15 h1(5.5);
  16
  17 h2(5.5);
  18
  19 translate(h1,v,h2)$
  20
  21 h1(5.5);
  22
  23 h2(5.5);
  24
  25 compile(functions);
  26
  27 h1(5.5);
  28
  29 h2(5.5);
  30