function aifaz=aifaz_fct(r,uz,t,E,A,l,w,s0)
c=abs(uz);
aifaz=(1/8*(-1/(w^2*r^2+4*c^2*t^2-4*w*r*c*t+4*l^2*t^2)^(1/2)*c*t*(2*w^2*r-4*w*c*t)+1/2/(w^2*r^2+4*c^2*t^2-4*w*r*c*t+4*l^2*t^2)^(1/2)*w*r*(2*w^2*r-4*w*c*t)+(w^2*r^2+4*c^2*t^2-4*w*r*c*t+4*l^2*t^2)^(1/2)*w+(4*(-2*w*c*t+2*w^2*r+1/2/((w^2*r^2+4*c^2*t^2-4*w*r*c*t+4*l^2*t^2)/l^2/t^2)^(1/2)*(w^2*r^2/t^2/l^4)^(1/2)*l*(2*w^2*r-4*w*c*t)+((w^2*r^2+4*c^2*t^2-4*w*r*c*t+4*l^2*t^2)/l^2/t^2)^(1/2)/(w^2*r^2/t^2/l^4)^(1/2)/l*w^2*r)/l^3/t^2/(w^2*r^2/t^2/l^4)^(1/2)-4*(-2*w*r*c*t+w^2*r^2+((w^2*r^2+4*c^2*t^2-4*w*r*c*t+4*l^2*t^2)/l^2/t^2)^(1/2)*(w^2*r^2/t^2/l^4)^(1/2)*l^3*t^2)/l^7/t^4/(w^2*r^2/t^2/l^4)^(3/2)*w^2*r)/(-2*w*r*c*t+w^2*r^2+((w^2*r^2+4*c^2*t^2-4*w*r*c*t+4*l^2*t^2)/l^2/t^2)^(1/2)*(w^2*r^2/t^2/l^4)^(1/2)*l^3*t^2)*l^5*t^4*(w^2*r^2/t^2/l^4)^(1/2)+1/(w^2*r^2+4*c^2*t^2+4*w*r*c*t+4*l^2*t^2)^(1/2)*c*t*(2*w^2*r+4*w*c*t)+1/2/(w^2*r^2+4*c^2*t^2+4*w*r*c*t+4*l^2*t^2)^(1/2)*w*r*(2*w^2*r+4*w*c*t)+(w^2*r^2+4*c^2*t^2+4*w*r*c*t+4*l^2*t^2)^(1/2)*w-(4*(-2*w*c*t-2*w^2*r+1/2/((w^2*r^2+4*c^2*t^2+4*w*r*c*t+4*l^2*t^2)/l^2/t^2)^(1/2)*(w^2*r^2/t^2/l^4)^(1/2)*l*(2*w^2*r+4*w*c*t)+((w^2*r^2+4*c^2*t^2+4*w*r*c*t+4*l^2*t^2)/l^2/t^2)^(1/2)/(w^2*r^2/t^2/l^4)^(1/2)/l*w^2*r)/l^3/t^2/(w^2*r^2/t^2/l^4)^(1/2)-4*(-2*w*r*c*t-w^2*r^2+((w^2*r^2+4*c^2*t^2+4*w*r*c*t+4*l^2*t^2)/l^2/t^2)^(1/2)*(w^2*r^2/t^2/l^4)^(1/2)*l^3*t^2)/l^7/t^4/(w^2*r^2/t^2/l^4)^(3/2)*w^2*r)/(-2*w*r*c*t-w^2*r^2+((w^2*r^2+4*c^2*t^2+4*w*r*c*t+4*l^2*t^2)/l^2/t^2)^(1/2)*(w^2*r^2/t^2/l^4)^(1/2)*l^3*t^2)*l^5*t^4*(w^2*r^2/t^2/l^4)^(1/2))/t/r/w-1/8*(-2*(w^2*r^2+4*c^2*t^2-4*w*r*c*t+4*l^2*t^2)^(1/2)*c*t+(w^2*r^2+4*c^2*t^2-4*w*r*c*t+4*l^2*t^2)^(1/2)*w*r+4*log(4*(-2*w*r*c*t+w^2*r^2+((w^2*r^2+4*c^2*t^2-4*w*r*c*t+4*l^2*t^2)/l^2/t^2)^(1/2)*(w^2*r^2/t^2/l^4)^(1/2)*l^3*t^2)/l^3/t^2/(w^2*r^2/t^2/l^4)^(1/2))*l^2*t^2+2*(w^2*r^2+4*c^2*t^2+4*w*r*c*t+4*l^2*t^2)^(1/2)*c*t+(w^2*r^2+4*c^2*t^2+4*w*r*c*t+4*l^2*t^2)^(1/2)*w*r-4*log(4*(-2*w*r*c*t-w^2*r^2+((w^2*r^2+4*c^2*t^2+4*w*r*c*t+4*l^2*t^2)/l^2/t^2)^(1/2)*(w^2*r^2/t^2/l^4)^(1/2)*l^3*t^2)/l^3/t^2/(w^2*r^2/t^2/l^4)^(1/2))*l^2*t^2)/t/r^2/w-t/E/A*(-l^2/r^2-c^2/r^2+1/12*w^2/t^2))*uz/r+(1/8*(-8/(w^2*r^2+4*c^2*t^2-4*w*r*c*t+4*l^2*t^2)^(1/2)*c*t^3*l+4/(w^2*r^2+4*c^2*t^2-4*w*r*c*t+4*l^2*t^2)^(1/2)*w*r*l*t^2+(4*(1/2/((w^2*r^2+4*c^2*t^2-4*w*r*c*t+4*l^2*t^2)/l^2/t^2)^(1/2)*(w^2*r^2/t^2/l^4)^(1/2)*l^3*t^2*(8/l-2*(w^2*r^2+4*c^2*t^2-4*w*r*c*t+4*l^2*t^2)/l^3/t^2)-2*((w^2*r^2+4*c^2*t^2-4*w*r*c*t+4*l^2*t^2)/l^2/t^2)^(1/2)/(w^2*r^2/t^2/l^4)^(1/2)/l^2*w^2*r^2+3*((w^2*r^2+4*c^2*t^2-4*w*r*c*t+4*l^2*t^2)/l^2/t^2)^(1/2)*(w^2*r^2/t^2/l^4)^(1/2)*l^2*t^2)/l^3/t^2/(w^2*r^2/t^2/l^4)^(1/2)-12*(-2*w*r*c*t+w^2*r^2+((w^2*r^2+4*c^2*t^2-4*w*r*c*t+4*l^2*t^2)/l^2/t^2)^(1/2)*(w^2*r^2/t^2/l^4)^(1/2)*l^3*t^2)/l^4/t^2/(w^2*r^2/t^2/l^4)^(1/2)+8*(-2*w*r*c*t+w^2*r^2+((w^2*r^2+4*c^2*t^2-4*w*r*c*t+4*l^2*t^2)/l^2/t^2)^(1/2)*(w^2*r^2/t^2/l^4)^(1/2)*l^3*t^2)/l^8/t^4/(w^2*r^2/t^2/l^4)^(3/2)*w^2*r^2)/(-2*w*r*c*t+w^2*r^2+((w^2*r^2+4*c^2*t^2-4*w*r*c*t+4*l^2*t^2)/l^2/t^2)^(1/2)*(w^2*r^2/t^2/l^4)^(1/2)*l^3*t^2)*l^5*t^4*(w^2*r^2/t^2/l^4)^(1/2)+8*log(4*(-2*w*r*c*t+w^2*r^2+((w^2*r^2+4*c^2*t^2-4*w*r*c*t+4*l^2*t^2)/l^2/t^2)^(1/2)*(w^2*r^2/t^2/l^4)^(1/2)*l^3*t^2)/l^3/t^2/(w^2*r^2/t^2/l^4)^(1/2))*l*t^2+8/(w^2*r^2+4*c^2*t^2+4*w*r*c*t+4*l^2*t^2)^(1/2)*c*t^3*l+4/(w^2*r^2+4*c^2*t^2+4*w*r*c*t+4*l^2*t^2)^(1/2)*w*r*l*t^2-(4*(1/2/((w^2*r^2+4*c^2*t^2+4*w*r*c*t+4*l^2*t^2)/l^2/t^2)^(1/2)*(w^2*r^2/t^2/l^4)^(1/2)*l^3*t^2*(8/l-2*(w^2*r^2+4*c^2*t^2+4*w*r*c*t+4*l^2*t^2)/l^3/t^2)-2*((w^2*r^2+4*c^2*t^2+4*w*r*c*t+4*l^2*t^2)/l^2/t^2)^(1/2)/(w^2*r^2/t^2/l^4)^(1/2)/l^2*w^2*r^2+3*((w^2*r^2+4*c^2*t^2+4*w*r*c*t+4*l^2*t^2)/l^2/t^2)^(1/2)*(w^2*r^2/t^2/l^4)^(1/2)*l^2*t^2)/l^3/t^2/(w^2*r^2/t^2/l^4)^(1/2)-12*(-2*w*r*c*t-w^2*r^2+((w^2*r^2+4*c^2*t^2+4*w*r*c*t+4*l^2*t^2)/l^2/t^2)^(1/2)*(w^2*r^2/t^2/l^4)^(1/2)*l^3*t^2)/l^4/t^2/(w^2*r^2/t^2/l^4)^(1/2)+8*(-2*w*r*c*t-w^2*r^2+((w^2*r^2+4*c^2*t^2+4*w*r*c*t+4*l^2*t^2)/l^2/t^2)^(1/2)*(w^2*r^2/t^2/l^4)^(1/2)*l^3*t^2)/l^8/t^4/(w^2*r^2/t^2/l^4)^(3/2)*w^2*r^2)/(-2*w*r*c*t-w^2*r^2+((w^2*r^2+4*c^2*t^2+4*w*r*c*t+4*l^2*t^2)/l^2/t^2)^(1/2)*(w^2*r^2/t^2/l^4)^(1/2)*l^3*t^2)*l^5*t^4*(w^2*r^2/t^2/l^4)^(1/2)-8*log(4*(-2*w*r*c*t-w^2*r^2+((w^2*r^2+4*c^2*t^2+4*w*r*c*t+4*l^2*t^2)/l^2/t^2)^(1/2)*(w^2*r^2/t^2/l^4)^(1/2)*l^3*t^2)/l^3/t^2/(w^2*r^2/t^2/l^4)^(1/2))*l*t^2)/t/r/w-2*t/E/A*l/r)*uz/l+(1/8*(-1/(w^2*r^2+4*c^2*t^2-4*w*r*c*t+4*l^2*t^2)^(1/2)*c*t*(8*c*t^2-4*w*r*t)-2*(w^2*r^2+4*c^2*t^2-4*w*r*c*t+4*l^2*t^2)^(1/2)*t+1/2/(w^2*r^2+4*c^2*t^2-4*w*r*c*t+4*l^2*t^2)^(1/2)*w*r*(8*c*t^2-4*w*r*t)+4*(-2*w*r*t+1/2/((w^2*r^2+4*c^2*t^2-4*w*r*c*t+4*l^2*t^2)/l^2/t^2)^(1/2)*(w^2*r^2/t^2/l^4)^(1/2)*l*(8*c*t^2-4*w*r*t))/(-2*w*r*c*t+w^2*r^2+((w^2*r^2+4*c^2*t^2-4*w*r*c*t+4*l^2*t^2)/l^2/t^2)^(1/2)*(w^2*r^2/t^2/l^4)^(1/2)*l^3*t^2)*l^2*t^2+1/(w^2*r^2+4*c^2*t^2+4*w*r*c*t+4*l^2*t^2)^(1/2)*c*t*(8*c*t^2+4*w*r*t)+2*(w^2*r^2+4*c^2*t^2+4*w*r*c*t+4*l^2*t^2)^(1/2)*t+1/2/(w^2*r^2+4*c^2*t^2+4*w*r*c*t+4*l^2*t^2)^(1/2)*w*r*(8*c*t^2+4*w*r*t)-4*(-2*w*r*t+1/2/((w^2*r^2+4*c^2*t^2+4*w*r*c*t+4*l^2*t^2)/l^2/t^2)^(1/2)*(w^2*r^2/t^2/l^4)^(1/2)*l*(8*c*t^2+4*w*r*t))/(-2*w*r*c*t-w^2*r^2+((w^2*r^2+4*c^2*t^2+4*w*r*c*t+4*l^2*t^2)/l^2/t^2)^(1/2)*(w^2*r^2/t^2/l^4)^(1/2)*l^3*t^2)*l^2*t^2)/t/r/w-2*t/E/A*c/r)*sign(uz);