E1 = 3300;
E2 = 29700;
v1 = 0.41;
v2 = 0.28;
E=((1-v1^2)/E1)+((1-v2^2)/E2);
H = 7.2087e-07;
beta = 6.5;
sigma = 0.00331;
psi=((H*beta)/(E*sigma));
row = 0.1;
z1 = 0.00168;
z2 = 0.00387;
z3 = 0.00265;
c=2*z2-z1-z3;
C=((-beta*log(row))^2)/(c*sigma);
deltac=C/2;
gamma1 = 22.14;
gamma2 = 22.14;
gamma12 = 0;
gamma=gamma1+gamma2-gamma12;
phi=((gamma*beta)/(E*(sigma^2)));
deltac1 = 12.25207;
h = 0.00289;
syms Z
P=0.532*(1/psi)*beta*((1/((C^1/2)*(-log(row))))*(int(exp(-Z^2/2)*(Z - h)^(3/2), Z, h, (h - deltac1)))-(3.76*phi^(1/2)*(1/((C^1/2)*(-log(row))^(1/2))))*(int(exp(-Z^2/2)*(Z - h)^(3/4), Z, h, (h - deltac1))))+(2.5)*((1/C)*((int(exp(-Z^2/2)*(Z - h), Z, (h+deltac), inf)))-(phi/psi)*(1/C)*((int(exp(-Z^2/2), Z, (h+deltac), inf))));
Pval = vpa(P);
disp(Pval)
Best Answer