0
|
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 |