[X,Y,Z] = meshgrid(-3:.02:3,-3:.02:3,2);
[ue,un,uv] = TDdispFS(X,Y,Z,[-1 0 0],[1 -1 -1],[0 1.5 .5],-1,2,3,.25);

[X,Y,Z] = meshgrid(-3:.02:3,-3:.02:3,-5);
[ue,un,uv] = TDdispHS(X,Y,Z,[-1 0 0],[1 -1 -1],[0 1.5 -2],-1,2,3,.25);

[X,Y,Z] = meshgrid(-3:.02:3,-3:.02:3,2);
[Stress,Strain] = TDstressFS(X,Y,Z,[-1 0 0],[1 -1 -1],[0 1.5 .5],...
-1,2,3,.33e11,.33e11);

[X,Y,Z] = meshgrid(-3:.02:3,-3:.02:3,-5);
[Stress,Strain] = TDstressHS(X,Y,Z,[-1 0 0],[1 -1 -1],[0 1.5 -2],...
-1,2,3,.33e11,.33e11);