Re: Using quad() for multidimensional integration

[Vinayak Dutt]

Ted Harding writes:
#( Re Message From: Vinayak Dutt )
#> 
#Since posting my previous comment on the above, I ran my little test program
#
#    ( for  int_0^1  int_0^x  y dy dx )
#
#on matlab (an old version).
#
#The result, after a few seconds, was a stream of messages all like
#
#     Recursion level limit reached in quad. Singularity likely.
#
#and, when interrupted, it backed up through several layers of calls to
#"quadstp".
#
#This adds support to my suspicion that this problem may be due to running
#out of stack space.

 I think you might be right. I ran a simpler double integral which has a faster
convergence, and, it did work. So this might be a stack space problem.


--vinayak-
/*
 * vinayak dutt
 * ex-graduate student, ultrasound research
 * mayo graduate school, rochester mn
 *
 * e-mail: address@hidden
 *         address@hidden
 *         address@hidden
 *
 */
#include "disclaimer.h"