Numerical Integration

[Ivan Sutoris]

On Sat, Dec 6, 2008 at 6:26 PM, Marcin Ślęczka <marcin.sleczka at wp.pl> wrote:
>
> I'm a beginner on using the Octave. I've a problem with numerical integration
> of such a function:
> \int_0^1 dt \frac{1}{t^{3/2}}\exp[-iR(t)]\frac{1}{(1+R^2(t))^2} where
> R(t)=\int_0^t dx e^{-x^2}\sin(at+b) (tex code)
> I've no idea how to integrate this. Every help will be usefull.
> Thanks for help

Hi

You can integrate univariate function with 'quad' command. You need to
define function under integral (e.g. by programming it in separate
m-file), see 'help quad' for details. Looking at your function, R(t)
seems to be a multiply of error function which is available directly
in Octave (see 'help erf'), so you don't have to compute the second
integral.

Regards
Ivan Sutoris