You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
19 lines
400 B
19 lines
400 B
11 months ago
|
function intB = GaussIntB(h)
|
||
|
intB = zeros(6,24);
|
||
|
x0 = 0; y0 = 0; z0=0;
|
||
|
numgp=2;
|
||
|
gpx = getgp(x0,x0+h,numgp);
|
||
|
gpy = getgp(y0,y0+h,numgp);
|
||
|
gpz = getgp(z0,z0+h,numgp);
|
||
|
for ix = 1:numgp*2
|
||
|
x=gpx(ix);
|
||
|
for iy = 1:numgp*2
|
||
|
y=gpy(iy);
|
||
|
for iz = 1:numgp*2
|
||
|
z=gpz(iz);
|
||
|
intB = intB + LocalB(x,y,z,x0,y0,z0,h);
|
||
|
end
|
||
|
end
|
||
|
end
|
||
|
intB = intB.*h^3/(numgp*2)^3;
|
||
|
end
|