function B=LocalB(x,y,z,x0,y0,z0,h)
B=zeros(6,24);
B(1,1)=-(1.0*(h - 1.0*y + y0)*(h - 1.0*z + z0))/h^3;
B(1,2)=0.0;
B(1,3)=0.0;
B(1,4)=((h - 1.0*y + y0)*(h - 1.0*z + z0))/h^3;
B(1,5)=0.0;
B(1,6)=0.0;
B(1,7)=-(1.0*(y - 1.0*y0)*(h - 1.0*z + z0))/h^3;
B(1,8)=0.0;
B(1,9)=0.0;
B(1,10)=((y - 1.0*y0)*(h - 1.0*z + z0))/h^3;
B(1,11)=0.0;
B(1,12)=0.0;
B(1,13)=-(1.0*(z - 1.0*z0)*(h - 1.0*y + y0))/h^3;
B(1,14)=0.0;
B(1,15)=0.0;
B(1,16)=((z - 1.0*z0)*(h - 1.0*y + y0))/h^3;
B(1,17)=0.0;
B(1,18)=0.0;
B(1,19)=-(1.0*(y - 1.0*y0)*(z - 1.0*z0))/h^3;
B(1,20)=0.0;
B(1,21)=0.0;
B(1,22)=((y - 1.0*y0)*(z - 1.0*z0))/h^3;
B(1,23)=0.0;
B(1,24)=0.0;
B(2,1)=0.0;
B(2,2)=-(1.0*(h - 1.0*x + x0)*(h - 1.0*z + z0))/h^3;
B(2,3)=0.0;
B(2,4)=0.0;
B(2,5)=-(1.0*(x - 1.0*x0)*(h - 1.0*z + z0))/h^3;
B(2,6)=0.0;
B(2,7)=0.0;
B(2,8)=((h - 1.0*x + x0)*(h - 1.0*z + z0))/h^3;
B(2,9)=0.0;
B(2,10)=0.0;
B(2,11)=((x - 1.0*x0)*(h - 1.0*z + z0))/h^3;
B(2,12)=0.0;
B(2,13)=0.0;
B(2,14)=-(1.0*(z - 1.0*z0)*(h - 1.0*x + x0))/h^3;
B(2,15)=0.0;
B(2,16)=0.0;
B(2,17)=-(1.0*(x - 1.0*x0)*(z - 1.0*z0))/h^3;
B(2,18)=0.0;
B(2,19)=0.0;
B(2,20)=((z - 1.0*z0)*(h - 1.0*x + x0))/h^3;
B(2,21)=0.0;
B(2,22)=0.0;
B(2,23)=((x - 1.0*x0)*(z - 1.0*z0))/h^3;
B(2,24)=0.0;
B(3,1)=0.0;
B(3,2)=0.0;
B(3,3)=-(1.0*(h - 1.0*x + x0)*(h - 1.0*y + y0))/h^3;
B(3,4)=0.0;
B(3,5)=0.0;
B(3,6)=-(1.0*(x - 1.0*x0)*(h - 1.0*y + y0))/h^3;
B(3,7)=0.0;
B(3,8)=0.0;
B(3,9)=-(1.0*(y - 1.0*y0)*(h - 1.0*x + x0))/h^3;
B(3,10)=0.0;
B(3,11)=0.0;
B(3,12)=-(1.0*(x - 1.0*x0)*(y - 1.0*y0))/h^3;
B(3,13)=0.0;
B(3,14)=0.0;
B(3,15)=((h - 1.0*x + x0)*(h - 1.0*y + y0))/h^3;
B(3,16)=0.0;
B(3,17)=0.0;
B(3,18)=((x - 1.0*x0)*(h - 1.0*y + y0))/h^3;
B(3,19)=0.0;
B(3,20)=0.0;
B(3,21)=((y - 1.0*y0)*(h - 1.0*x + x0))/h^3;
B(3,22)=0.0;
B(3,23)=0.0;
B(3,24)=((x - 1.0*x0)*(y - 1.0*y0))/h^3;
B(4,1)=-(1.0*(h - 1.0*x + x0)*(h - 1.0*z + z0))/h^3;
B(4,2)=-(1.0*(h - 1.0*y + y0)*(h - 1.0*z + z0))/h^3;
B(4,3)=0.0;
B(4,4)=-(1.0*(x - 1.0*x0)*(h - 1.0*z + z0))/h^3;
B(4,5)=((h - 1.0*y + y0)*(h - 1.0*z + z0))/h^3;
B(4,6)=0.0;
B(4,7)=((h - 1.0*x + x0)*(h - 1.0*z + z0))/h^3;
B(4,8)=-(1.0*(y - 1.0*y0)*(h - 1.0*z + z0))/h^3;
B(4,9)=0.0;
B(4,10)=((x - 1.0*x0)*(h - 1.0*z + z0))/h^3;
B(4,11)=((y - 1.0*y0)*(h - 1.0*z + z0))/h^3;
B(4,12)=0.0;
B(4,13)=-(1.0*(z - 1.0*z0)*(h - 1.0*x + x0))/h^3;
B(4,14)=-(1.0*(z - 1.0*z0)*(h - 1.0*y + y0))/h^3;
B(4,15)=0.0;
B(4,16)=-(1.0*(x - 1.0*x0)*(z - 1.0*z0))/h^3;
B(4,17)=((z - 1.0*z0)*(h - 1.0*y + y0))/h^3;
B(4,18)=0.0;
B(4,19)=((z - 1.0*z0)*(h - 1.0*x + x0))/h^3;
B(4,20)=-(1.0*(y - 1.0*y0)*(z - 1.0*z0))/h^3;
B(4,21)=0.0;
B(4,22)=((x - 1.0*x0)*(z - 1.0*z0))/h^3;
B(4,23)=((y - 1.0*y0)*(z - 1.0*z0))/h^3;
B(4,24)=0.0;
B(5,1)=0.0;
B(5,2)=-(1.0*(h - 1.0*x + x0)*(h - 1.0*y + y0))/h^3;
B(5,3)=-(1.0*(h - 1.0*x + x0)*(h - 1.0*z + z0))/h^3;
B(5,4)=0.0;
B(5,5)=-(1.0*(x - 1.0*x0)*(h - 1.0*y + y0))/h^3;
B(5,6)=-(1.0*(x - 1.0*x0)*(h - 1.0*z + z0))/h^3;
B(5,7)=0.0;
B(5,8)=-(1.0*(y - 1.0*y0)*(h - 1.0*x + x0))/h^3;
B(5,9)=((h - 1.0*x + x0)*(h - 1.0*z + z0))/h^3;
B(5,10)=0.0;
B(5,11)=-(1.0*(x - 1.0*x0)*(y - 1.0*y0))/h^3;
B(5,12)=((x - 1.0*x0)*(h - 1.0*z + z0))/h^3;
B(5,13)=0.0;
B(5,14)=((h - 1.0*x + x0)*(h - 1.0*y + y0))/h^3;
B(5,15)=-(1.0*(z - 1.0*z0)*(h - 1.0*x + x0))/h^3;
B(5,16)=0.0;
B(5,17)=((x - 1.0*x0)*(h - 1.0*y + y0))/h^3;
B(5,18)=-(1.0*(x - 1.0*x0)*(z - 1.0*z0))/h^3;
B(5,19)=0.0;
B(5,20)=((y - 1.0*y0)*(h - 1.0*x + x0))/h^3;
B(5,21)=((z - 1.0*z0)*(h - 1.0*x + x0))/h^3;
B(5,22)=0.0;
B(5,23)=((x - 1.0*x0)*(y - 1.0*y0))/h^3;
B(5,24)=((x - 1.0*x0)*(z - 1.0*z0))/h^3;
B(6,1)=-(1.0*(h - 1.0*x + x0)*(h - 1.0*y + y0))/h^3;
B(6,2)=0.0;
B(6,3)=-(1.0*(h - 1.0*y + y0)*(h - 1.0*z + z0))/h^3;
B(6,4)=-(1.0*(x - 1.0*x0)*(h - 1.0*y + y0))/h^3;
B(6,5)=0.0;
B(6,6)=((h - 1.0*y + y0)*(h - 1.0*z + z0))/h^3;
B(6,7)=-(1.0*(y - 1.0*y0)*(h - 1.0*x + x0))/h^3;
B(6,8)=0.0;
B(6,9)=-(1.0*(y - 1.0*y0)*(h - 1.0*z + z0))/h^3;
B(6,10)=-(1.0*(x - 1.0*x0)*(y - 1.0*y0))/h^3;
B(6,11)=0.0;
B(6,12)=((y - 1.0*y0)*(h - 1.0*z + z0))/h^3;
B(6,13)=((h - 1.0*x + x0)*(h - 1.0*y + y0))/h^3;
B(6,14)=0.0;
B(6,15)=-(1.0*(z - 1.0*z0)*(h - 1.0*y + y0))/h^3;
B(6,16)=((x - 1.0*x0)*(h - 1.0*y + y0))/h^3;
B(6,17)=0.0;
B(6,18)=((z - 1.0*z0)*(h - 1.0*y + y0))/h^3;
B(6,19)=((y - 1.0*y0)*(h - 1.0*x + x0))/h^3;
B(6,20)=0.0;
B(6,21)=-(1.0*(y - 1.0*y0)*(z - 1.0*z0))/h^3;
B(6,22)=((x - 1.0*x0)*(y - 1.0*y0))/h^3;
B(6,23)=0.0;
B(6,24)=((y - 1.0*y0)*(z - 1.0*z0))/h^3;
end