clear all;
clc;

mat = [1000, 0, 0; %Ox, Txy, Txz
       0, -1000, 500; %Tyx, Oy, Tyz
       0, 500, 0]; %Tzx, Tzy, Oz    kpsi
E = 30000; %kpsi OR 206.8 GPa
G = 11700; %kpsi
pois = .28;
%1 kpsi = 6894.76 kPa

PrinStress = eig(mat);
normPrinStress = sort(PrinStress, 'descend')


maxShear = (max(normPrinStress)-min(normPrinStress))/2

appNormStrain = [mat(1,1)/E - (pois*mat(2,2))/E - (pois*mat(3,3))/E; %x
               -(pois*mat(1,1))/E + mat(2,2)/E - (pois*mat(3,3))/E;   %y
               -(pois*mat(1,1))/E - (pois*mat(2,2))/E + mat(3,3)/E]   %z

appShearStrain = [mat(1,2)/G; %xy
              mat(3,2)/G;     %zy
              mat(1,3)/G]     %xz

prinStrain = [normPrinStress(1,1)/E - (pois*normPrinStress(2,1))/E - (pois*normPrinStress(3,1))/E; %x
               -(pois*normPrinStress(1,1))/E + normPrinStress(2,1)/E - (pois*normPrinStress(3,1))/E; %y
               -(pois*normPrinStress(1,1))/E - (pois*normPrinStress(2,1))/E + normPrinStress(3,1)/E] %z