format long
um=linspace(-0.03,0.03,1001);
[m,num]=size(um);
T=-100;
a1=(T-115)*4.124*10^5;
a11=-2.097*10^8;
a12=7.974*10^8;
a111=1.294*10^9;
a112=-1.950*10^9;
a123=-2.5*10^9;
a1111=3.863*10^10;
a1112=2.5293*10^10;
a1122=1.637*10^10;
a1123=1.367*10^10;
Q11=0.11;
Q12=-0.043;
Q44=0.059;
s11=8.3*10^-12;
s12=-2.7*10^-12;
s44=9.24*10^-12;
ar11=1/4*(2*a11+a12)+1/24*(2*(Q11-Q12+Q44)^2/(s11-s12+s44)+(4*Q11+8*Q12+Q44)^2/(4*s11+8*s12+s44));
ar33=1/3*(a11+a12)+(2*Q11+4*Q12-Q44)^2/6*(4*s11+8*s12+s44);
ar13=2*a11+1/6*((2*Q11+4*Q12-Q44)*(4*Q11+8*Q12+Q44)/(4*s11+8*s12+s44)+(2*Q11-2*Q12-Q44)^2/(s11-s12-s4