MATLAB program for quadrature in 2D

被引:70
作者
Shampine, L. F. [1 ]
机构
[1] So Methodist Univ, Dept Math, Dallas, TX 75275 USA
关键词
integral over plane region; MATLAB; vectorization; singularities; adaptive quadrature;
D O I
10.1016/j.amc.2008.02.012
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
We discuss here the algorithms of TwoD, a MATLAB program for approximating integrals over generalized rectangles and sectors. Capabilities of the language are exploited to make TwoD very easy to use. Vectorization is exploited by a novel algorithm to make the program efficient and reliable. TwoD provides for moderate singularities on boundaries of the region. (C) 2008 Elsevier Inc. All rights reserved.
引用
收藏
页码:266 / 274
页数:9
相关论文
共 13 条
[1]   AN ADAPTIVE ALGORITHM FOR THE APPROXIMATE CALCULATION OF MULTIPLE INTEGRALS [J].
BERNTSEN, J ;
ESPELID, TO ;
GENZ, A .
ACM TRANSACTIONS ON MATHEMATICAL SOFTWARE, 1991, 17 (04) :437-451
[2]   Algorithm 764: Cubpack++: A C++ package for automatic two-dimensional cubature [J].
Cools, R ;
Laurie, D ;
Pluym, L .
ACM TRANSACTIONS ON MATHEMATICAL SOFTWARE, 1997, 23 (01) :1-15
[3]  
DAVIS PJ, 1984, METHODS NUMERICAL IN
[4]  
KAHANER D, 1984, SOURCES DEV MATH SOF
[5]   TWODQD AN ADAPTIVE ROUTINE FOR TWO-DIMENSIONAL INTEGRATION [J].
KAHANER, DK ;
RECHARD, OW .
JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS, 1987, 17 (1-2) :215-234
[6]  
Kronrod A.S., 1965, NODES WEIGHTS QUADRA
[7]   WHEN NOT TO USE AN AUTOMATIC QUADRATURE ROUTINE [J].
LYNESS, JN .
SIAM REVIEW, 1983, 25 (01) :63-87
[8]  
*NUM ALG GROUP INC, NAG FORTR 90 LIB REL
[9]  
Piessens R., 1983, QUADPACK SUBROUTINE
[10]   AUTOMATIC COMPUTATION OF IMPROPER INTEGRALS OVER A BOUNDED OR UNBOUNDED PLANAR REGION [J].
ROBINSON, I ;
DEDONCKER, E .
COMPUTING, 1981, 27 (03) :253-284