dvald = 2.7*10^-3;

invpsi = (dvald/(2*g2))-(pi/z2)+((pi+4*x2*tan(alphat))/(2*z2)) + invat;

invAvec4 = [0.01356 0.01378];
Y4 = [19.4 19.6];
alphaSearch4 = 0.013719847817667;

alpha4 = interp1(invAvec4,Y4,alphaSearch4);
angle4 = alpha4*pi/180;

L2 = (dvald + 2*(g2/cos(angle4)))*10^3;
L3 = (dvald + 2*(g3/cos(psimin3)))*10^3;