-
Notifications
You must be signed in to change notification settings - Fork 3
/
Theta.m
9 lines (7 loc) · 986 Bytes
/
Theta.m
1
2
3
4
5
6
7
8
9
function Theta=Theta(Mat,X,Y,Mesh_X,Mesh_Y)
if ((1.0/12.0)*(Mat(X+Mesh_X,Y+Mesh_Y)-Mat(X-Mesh_X,Y+Mesh_Y)+4*Mat(X+Mesh_X,Y)-Mat(X-Mesh_X,Y)+Mat(X+Mesh_X,Y-Mesh_Y)-Mat(X-Mesh_X,Y-Mesh_Y))~=0 || (1.0/12.0)*(Mat(X+Mesh_X,Y+Mesh_Y)-Mat(X+Mesh_X,Y-Mesh_Y)+4*Mat(X,Y+Mesh_Y)-Mat(X,Y-Mesh_Y)+Mat(X-Mesh_X,Y+Mesh_Y)-Mat(X-Mesh_X,Y-Mesh_Y))~=0)
Theta=atan(((1.0/12.0)*(Mat(X+Mesh_X,Y+Mesh_Y)-Mat(X+Mesh_X,Y-Mesh_Y)+4*Mat(X,Y+Mesh_Y)-Mat(X,Y-Mesh_Y)+Mat(X-Mesh_X,Y+Mesh_Y)-Mat(X-Mesh_X,Y-Mesh_Y)))/((1.0/12.0)*(Mat(X+Mesh_X,Y+Mesh_Y)-Mat(X-Mesh_X,Y+Mesh_Y)+4*Mat(X+Mesh_X,Y)-Mat(X-Mesh_X,Y)+Mat(X+Mesh_X,Y-Mesh_Y)-Mat(X-Mesh_X,Y-Mesh_Y))));
elseif ((1.0/12.0)*(Mat(X+Mesh_X,Y+Mesh_Y)-Mat(X-Mesh_X,Y+Mesh_Y)+4*Mat(X+Mesh_X,Y)-Mat(X-Mesh_X,Y)+Mat(X+Mesh_X,Y-Mesh_Y)-Mat(X-Mesh_X,Y-Mesh_Y))==0 && (1.0/12.0)*(Mat(X+Mesh_X,Y+Mesh_Y)-Mat(X+Mesh_X,Y-Mesh_Y)+4*Mat(X,Y+Mesh_Y)-Mat(X,Y-Mesh_Y)+Mat(X-Mesh_X,Y+Mesh_Y)-Mat(X-Mesh_X,Y-Mesh_Y))==0)
Theta=0;
end
end