From owner-help-octave at bevo dot che dot wisc dot edu  Tue Nov 28 11:40:00 1995
Subject: Re: How to quad("sin(a*sin(x))*sin(x)",0,pi/2)?
From: Ian Searle <ian at flash dot ds dot boeing dot com>
To: niles at axp745 dot gsfc dot nasa dot gov
Cc: help-octave at bevo dot che dot wisc dot edu
Date: Tue, 28 Nov 1995 09:36:01 -0800


> for (i in 1:4)
  {
    quadr ("sin(p1*sin(x)).*sin(x)", 0, pi/2,1.e-6,0,i).Q
  }
    0.691  
    0.906  
    0.533  
   -0.104  

The above is a piece of Rlab code that does what you want (I
think). The point is: I translated quadr from Matlab source. It is
written by Kirill K. Pankratov (kirill at plume dot mit dot edu), and is a very
nice function. It ought to run in Octave with little or no mods ? I
don't have the original Matlab source handy, otherwise I would have
appended it.

-Ian Searle