Mercurial > forge
comparison extra/integration/test/run2dtests.m @ 0:6b33357c7561 octave-forge
Initial revision
| author | pkienzle |
|---|---|
| date | Wed, 10 Oct 2001 19:54:49 +0000 |
| parents | |
| children |
comparison
equal
deleted
inserted
replaced
| -1:000000000000 | 0:6b33357c7561 |
|---|---|
| 1 format long | |
| 2 % | |
| 3 % x^2+y^2 integrated from -1 to 1 in x and -1 to 1 in y = 8/3 | |
| 4 % | |
| 5 [bx,by,w]=grule2d(2,2); | |
| 6 v=gquad2d('gxy',-1,1,-1,1,bx,by,w) | |
| 7 % or | |
| 8 v=quad2dg('gxy',-1,1,-1,1) | |
| 9 % or | |
| 10 v=gquad2d('gxy',-1,1,-1,1,2,2) | |
| 11 % or | |
| 12 n= gquadnd('hx',[-1;-1],[1;1],[2;2]) | |
| 13 % or | |
| 14 n=quadndg('hx',[-1;-1],[1;1]) | |
| 15 correct_ans=8/3 | |
| 16 % | |
| 17 % x^2+y^2 integrated from 0 to 2 in x and 0 to 2 in y = 32/3 | |
| 18 % | |
| 19 [bx,by,w]=grule2d(2,2); | |
| 20 vol3=gquad2d('gxy',0,2,0,2,bx,by,w) | |
| 21 % or | |
| 22 v=quad2dg('gxy',0,2,0,2) | |
| 23 % or | |
| 24 v=gquad2d('gxy',0,2,0,2,2,2) | |
| 25 % or | |
| 26 n = gquadnd('hx',[0;0],[2;2],[2;2]) | |
| 27 % or | |
| 28 n = quadndg('hx',[0;0],[2;2]) | |
| 29 correct_ans=32/3 | |
| 30 % | |
| 31 % x^2+y^2 integrated from 0 to 3 in x and 0 to 3 in y = 54 | |
| 32 % | |
| 33 [bx,by,w]=grule2d(2,2); | |
| 34 v=gquad2d('gxy',0,3,0,3,bx,by,w) | |
| 35 % or | |
| 36 v=gquad2d('gxy',0,3,0,3,2,2) | |
| 37 % or | |
| 38 n = gquadnd('hx',[0;0],[3;3],[2;2]) | |
| 39 correct_ans=54 | |
| 40 % | |
| 41 % x^2+y^2 integrated from 0 to 3 in x and 0 to 3 in y = 54 | |
| 42 % with the general area of intergration (functional limits in x) | |
| 43 % | |
| 44 [bx2,wx,by2,wy]=grule2dgen(2,2); | |
| 45 v=gquad2dgen('gxy','gliml','glimh',0,3,bx2,wx,by2,wy) | |
| 46 % or | |
| 47 v=gquad2dgen('gxy','gliml','glimh',0,3,2,2) | |
| 48 % or | |
| 49 v=quad2dggen('gxy','gliml','glimh',0,3) | |
| 50 correct_ans=54 | |
| 51 % | |
| 52 % x^2+y^2 integrated from 0 to y in x and 0 to 2 in y = 16/3 | |
| 53 % | |
| 54 v=gquad2dgen('gxy','gliml','glimh2',0,2,bx2,wx,by2,wy) | |
| 55 % or | |
| 56 v=gquad2dgen('gxy','gliml','glimh2',0,2,2,2) | |
| 57 correct_ans=16/3 | |
| 58 % | |
| 59 % 1 integrated from -sqrt(4-y^2) to sqrt(4-y^2) in x | |
| 60 % and -2 to 2 in y = 4*pi -- area of circle with radius 2 or (pi r^2) | |
| 61 % | |
| 62 v=gquad2dgen('gxy1','lcrcl','lcrcu',-2,2,bx2,wx,by2,wy) | |
| 63 % or | |
| 64 v=gquad2dgen('gxy1','lcrcl','lcrcu',-2,2,2,2) | |
| 65 correct_ans=4*pi | |
| 66 % --- same problem better quadratue (more points) | |
| 67 % 1 integrated from -sqrt(4-y^2) to sqrt(4-y^2) in x | |
| 68 % and -2 to 2 in y = 4*pi -- area of circle with radius 2 or (pi r^2) | |
| 69 % | |
| 70 [bx3,wx3,by3,wy3]=grule2dgen(5,5); | |
| 71 v=gquad2dgen('gxy1','lcrcl','lcrcu',-2,2,bx2,wx,by2,wy) | |
| 72 % or | |
| 73 v=gquad2dgen('gxy1','lcrcl','lcrcu',-2,2,5,5) | |
| 74 correct_ans=4*pi | |
| 75 % | |
| 76 % sqrt(x^2+y^2) integrated from -sqrt(4-y^2) to sqrt(4-y^2) in x | |
| 77 % and -2 to 2 in y = 16*pi/3 | |
| 78 % | |
| 79 % Need higher order quadrature | |
| 80 [bx3,wx3,by3,wy3]=grule2dgen(10,10); | |
| 81 v=gquad2dgen('gxy2','lcrcl','lcrcu',-2,2,bx3,wx3,by3,wy3) | |
| 82 % or | |
| 83 v=gquad2dgen('gxy2','lcrcl','lcrcu',-2,2,10,10) | |
| 84 correct_ans=16*pi/3 | |
| 85 % | |
| 86 % 1/sqrt(1-x^2) integrated from -1 to 1 in x = pi | |
| 87 % | |
| 88 % Use Gauss-Chebyshev quadrature | |
| 89 [bpc,wfc]=crule(2); | |
| 90 a=gquad('gxy1',-1,1,1,bpc,wfc) | |
| 91 %or | |
| 92 a=quadc('gxy1',-1,1) | |
| 93 correct_ans=pi | |
| 94 % | |
| 95 % x^2/sqrt(1-x^2) integrated from -1 to 1 in x = pi/2 | |
| 96 % | |
| 97 % Use Gauss-Chebyshev quadrature | |
| 98 a=gquad('xsquar',-1,1,1,bpc,wfc) | |
| 99 a=quadc('xsquar',-1,1) | |
| 100 correct_ans=pi/2 |