program test
      real*8 der
      der(1) = -(pi**((-3.0d+0)/2.0d+0)*exp((-((alpha**2*x**2)/w**2)/2.0d+0)-(x**2/w**2)/2.0d+0-((alpha**2*ksi**2)/w**2)/2.0d+0-(ksi**2/w**2)/2.0d+0)*((pi*ksi*exp(((alpha**2*ksi**2)/w**2)/2.0d+0)-pi*exp(((alpha**2*ksi**2)/w**2)/2.0d+0)*x)*exp(((alpha**2*x**2)/w**2)/2.0d+0+(ksi*x)/w**2)*erf((alpha*x-alpha*ksi)/(sqrt(2)*w))+(pi*ksi*exp(((alpha**2*ksi**2)/w**2)/2.0d+0)-pi*exp(((alpha**2*ksi**2)/w**2)/2.0d+0)*x)*exp(((alpha**2*x**2)/w**2)/2.0d+0+(ksi*x)/w**2)+sqrt(2)*sqrt(pi)*alpha*w*exp((alpha**2*ksi*x)/w**2+(ksi*x)/w**2)))/(sqrt(2)*w**3)
      der(2) = -(sqrt((pi-2)*alpha**2+pi)*sqrt(pi*alpha**2+pi)*exp((-((alpha**2*x**2)/w**2)/2.0d+0)-(x**2/w**2)/2.0d+0-((alpha**2*ksi**2)/w**2)/2.0d+0-(ksi**2/w**2)/2.0d+0-(2*pi)/abs(alpha))*(sqrt(alpha**2+1)*(sqrt(pi)*((((6*pi**2-24*pi+16)*alpha**5+(10*pi**2-24*pi)*alpha**3+4*pi**2*alpha)*abs(alpha)*w*exp(((alpha**2*ksi**2)/w**2)/2.0d+0+(2*pi)/abs(alpha))*x+(((-6*pi**2)+24*pi-16)*alpha**5+(24*pi-10*pi**2)*alpha**3-4*pi**2*alpha)*abs(alpha)*ksi*w*exp(((alpha**2*ksi**2)/w**2)/2.0d+0+(2*pi)/abs(alpha)))*exp(((alpha**2*x**2)/w**2)/2.0d+0+(ksi*x)/w**2)*erf((alpha*x-alpha*ksi)/(sqrt(2)*w))+(((6*pi**2-24*pi+16)*alpha**5+(10*pi**2-24*pi)*alpha**3+4*pi**2*alpha)*abs(alpha)*w*exp(((alpha**2*ksi**2)/w**2)/2.0d+0+(2*pi)/abs(alpha))*x+(((-6*pi**2)+24*pi-16)*alpha**5+(24*pi-10*pi**2)*alpha**3-4*pi**2*alpha)*abs(alpha)*ksi*w*exp(((alpha**2*ksi**2)/w**2)/2.0d+0+(2*pi)/abs(alpha)))*exp(((alpha**2*x**2)/w**2)/2.0d+0+(ksi*x)/w**2))+(((-3*2**(3.0d+0/2.0d+0)*pi**2)+3*2**(7.0d+0/2.0d+0)*pi-2**(9.0d+0/2.0d+0))*alpha**6+(3*2**(7.0d+0/2.0d+0)*pi-5*2**(3.0d+0/2.0d+0)*pi**2)*alpha**4-2**(5.0d+0/2.0d+0)*pi**2*alpha**2)*abs(alpha)*w**2*exp((alpha**2*ksi*x)/w**2+(ksi*x)/w**2+(2*pi)/abs(alpha)))+((((-2**(3.0d+0/2.0d+0)*pi**3)+2**(7.0d+0/2.0d+0)*pi**2-2**(7.0d+0/2.0d+0)*pi)*alpha**6+((-3*2**(3.0d+0/2.0d+0)*pi**3)+2**(9.0d+0/2.0d+0)*pi**2-2**(7.0d+0/2.0d+0)*pi)*alpha**4+(2**(7.0d+0/2.0d+0)*pi**2-3*2**(3.0d+0/2.0d+0)*pi**3)*alpha**2-2**(3.0d+0/2.0d+0)*pi**3)*abs(alpha)*exp(((alpha**2*ksi**2)/w**2)/2.0d+0+(2*pi)/abs(alpha))*x**2+((((-sqrt(2)*pi**3)+2**(5.0d+0/2.0d+0)*pi**2-2**(5.0d+0/2.0d+0)*pi)*alpha**7+((-3*sqrt(2)*pi**3)+2**(7.0d+0/2.0d+0)*pi**2-2**(5.0d+0/2.0d+0)*pi)*alpha**5+(2**(5.0d+0/2.0d+0)*pi**2-3*sqrt(2)*pi**3)*alpha**3-sqrt(2)*pi**3*alpha)*w+((2**(5.0d+0/2.0d+0)*pi**3-2**(9.0d+0/2.0d+0)*pi**2+2**(9.0d+0/2.0d+0)*pi)*alpha**6+(3*2**(5.0d+0/2.0d+0)*pi**3-2**(1.1000000000000001d+1/2.0d+0)*pi**2+2**(9.0d+0/2.0d+0)*pi)*alpha**4+(3*2**(5.0d+0/2.0d+0)*pi**3-2**(9.0d+0/2.0d+0)*pi**2)*alpha**2+2**(5.0d+0/2.0d+0)*pi**3)*exp((2*pi)/abs(alpha))*abs(alpha)*ksi)*exp(((alpha**2*ksi**2)/w**2)/2.0d+0)*x+(((2**(3.0d+0/2.0d+0)*pi**3-2**(7.0d+0/2.0d+0)*pi**2+2**(7.0d+0/2.0d+0)*pi)*alpha**6+(3*2**(3.0d+0/2.0d+0)*pi**3-2**(9.0d+0/2.0d+0)*pi**2+2**(7.0d+0/2.0d+0)*pi)*alpha**4+(3*2**(3.0d+0/2.0d+0)*pi**3-2**(7.0d+0/2.0d+0)*pi**2)*alpha**2+2**(3.0d+0/2.0d+0)*pi**3)*exp((2*pi)/abs(alpha))*abs(alpha)*w**2+((sqrt(2)*pi**3-2**(5.0d+0/2.0d+0)*pi**2+2**(5.0d+0/2.0d+0)*pi)*alpha**7+(3*sqrt(2)*pi**3-2**(7.0d+0/2.0d+0)*pi**2+2**(5.0d+0/2.0d+0)*pi)*alpha**5+(3*sqrt(2)*pi**3-2**(5.0d+0/2.0d+0)*pi**2)*alpha**3+sqrt(2)*pi**3*alpha)*ksi*w+(((-2**(3.0d+0/2.0d+0)*pi**3)+2**(7.0d+0/2.0d+0)*pi**2-2**(7.0d+0/2.0d+0)*pi)*alpha**6+((-3*2**(3.0d+0/2.0d+0)*pi**3)+2**(9.0d+0/2.0d+0)*pi**2-2**(7.0d+0/2.0d+0)*pi)*alpha**4+(2**(7.0d+0/2.0d+0)*pi**2-3*2**(3.0d+0/2.0d+0)*pi**3)*alpha**2-2**(3.0d+0/2.0d+0)*pi**3)*exp((2*pi)/abs(alpha))*abs(alpha)*ksi**2)*exp(((alpha**2*ksi**2)/w**2)/2.0d+0))*exp(((alpha**2*x**2)/w**2)/2.0d+0+(ksi*x)/w**2)*erf((alpha*x-alpha*ksi)/(sqrt(2)*w))+((((-2**(3.0d+0/2.0d+0)*pi**3)+2**(7.0d+0/2.0d+0)*pi**2-2**(7.0d+0/2.0d+0)*pi)*alpha**6+((-3*2**(3.0d+0/2.0d+0)*pi**3)+2**(9.0d+0/2.0d+0)*pi**2-2**(7.0d+0/2.0d+0)*pi)*alpha**4+(2**(7.0d+0/2.0d+0)*pi**2-3*2**(3.0d+0/2.0d+0)*pi**3)*alpha**2-2**(3.0d+0/2.0d+0)*pi**3)*abs(alpha)*exp(((alpha**2*ksi**2)/w**2)/2.0d+0+(2*pi)/abs(alpha))*x**2+((((-sqrt(2)*pi**3)+2**(5.0d+0/2.0d+0)*pi**2-2**(5.0d+0/2.0d+0)*pi)*alpha**7+((-3*sqrt(2)*pi**3)+2**(7.0d+0/2.0d+0)*pi**2-2**(5.0d+0/2.0d+0)*pi)*alpha**5+(2**(5.0d+0/2.0d+0)*pi**2-3*sqrt(2)*pi**3)*alpha**3-sqrt(2)*pi**3*alpha)*w+((2**(5.0d+0/2.0d+0)*pi**3-2**(9.0d+0/2.0d+0)*pi**2+2**(9.0d+0/2.0d+0)*pi)*alpha**6+(3*2**(5.0d+0/2.0d+0)*pi**3-2**(1.1000000000000001d+1/2.0d+0)*pi**2+2**(9.0d+0/2.0d+0)*pi)*alpha**4+(3*2**(5.0d+0/2.0d+0)*pi**3-2**(9.0d+0/2.0d+0)*pi**2)*alpha**2+2**(5.0d+0/2.0d+0)*pi**3)*exp((2*pi)/abs(alpha))*abs(alpha)*ksi)*exp(((alpha**2*ksi**2)/w**2)/2.0d+0)*x+(((2**(3.0d+0/2.0d+0)*pi**3-2**(7.0d+0/2.0d+0)*pi**2+2**(7.0d+0/2.0d+0)*pi)*alpha**6+(3*2**(3.0d+0/2.0d+0)*pi**3-2**(9.0d+0/2.0d+0)*pi**2+2**(7.0d+0/2.0d+0)*pi)*alpha**4+(3*2**(3.0d+0/2.0d+0)*pi**3-2**(7.0d+0/2.0d+0)*pi**2)*alpha**2+2**(3.0d+0/2.0d+0)*pi**3)*exp((2*pi)/abs(alpha))*abs(alpha)*w**2+((sqrt(2)*pi**3-2**(5.0d+0/2.0d+0)*pi**2+2**(5.0d+0/2.0d+0)*pi)*alpha**7+(3*sqrt(2)*pi**3-2**(7.0d+0/2.0d+0)*pi**2+2**(5.0d+0/2.0d+0)*pi)*alpha**5+(3*sqrt(2)*pi**3-2**(5.0d+0/2.0d+0)*pi**2)*alpha**3+sqrt(2)*pi**3*alpha)*ksi*w+(((-2**(3.0d+0/2.0d+0)*pi**3)+2**(7.0d+0/2.0d+0)*pi**2-2**(7.0d+0/2.0d+0)*pi)*alpha**6+((-3*2**(3.0d+0/2.0d+0)*pi**3)+2**(9.0d+0/2.0d+0)*pi**2-2**(7.0d+0/2.0d+0)*pi)*alpha**4+(2**(7.0d+0/2.0d+0)*pi**2-3*2**(3.0d+0/2.0d+0)*pi**3)*alpha**2-2**(3.0d+0/2.0d+0)*pi**3)*exp((2*pi)/abs(alpha))*abs(alpha)*ksi**2)*exp(((alpha**2*ksi**2)/w**2)/2.0d+0))*exp(((alpha**2*x**2)/w**2)/2.0d+0+(ksi*x)/w**2)+sqrt(pi)*(((4*pi**2-16*pi+16)*alpha**7+(12*pi**2-32*pi+16)*alpha**5+(12*pi**2-16*pi)*alpha**3+4*pi**2*alpha)*exp((2*pi)/abs(alpha))*abs(alpha)*w*x+((2*pi**2-8*pi+8)*alpha**8+(6*pi**2-16*pi+8)*alpha**6+(6*pi**2-8*pi)*alpha**4+2*pi**2*alpha**2)*w**2+(((-4*pi**2)+16*pi-16)*alpha**7+((-12*pi**2)+32*pi-16)*alpha**5+(16*pi-12*pi**2)*alpha**3-4*pi**2*alpha)*exp((2*pi)/abs(alpha))*abs(alpha)*ksi*w)*exp((alpha**2*ksi*x)/w**2+(ksi*x)/w**2)))/(sqrt(pi)*((4*pi**4-24*pi**3+48*pi**2-32*pi)*alpha**8+(16*pi**4-72*pi**3+96*pi**2-32*pi)*alpha**6+(24*pi**4-72*pi**3+48*pi**2)*alpha**4+(16*pi**4-24*pi**3)*alpha**2+4*pi**4)*abs(alpha)*w**4)
      der(3) = (exp((-((alpha**2*x**2)/w**2)/2.0d+0)-(x**2/w**2)/2.0d+0-((alpha**2*ksi**2)/w**2)/2.0d+0-(ksi**2/w**2)/2.0d+0-(2*pi)/abs(alpha))*(sqrt(pi*alpha**2+pi)*(sqrt(alpha**2+1)*(sqrt(pi)*((((4*pi**2-16*pi+16)*alpha**7+(8*pi**2-16*pi)*alpha**5+4*pi**2*alpha**3)*abs(alpha)*sigma*exp(((alpha**2*ksi**2)/w**2)/2.0d+0+(2*pi)/abs(alpha))*x**2+((((2*pi**4-12*pi**3+24*pi**2-16*pi)*alpha**8+(8*pi**4-36*pi**3+48*pi**2-16*pi)*alpha**6+(12*pi**4-36*pi**3+24*pi**2)*alpha**4+(8*pi**4-12*pi**3)*alpha**2+2*pi**4)*abs(alpha)+(2*pi**2-8*pi+8)*alpha**8+(4*pi**2-8*pi)*alpha**6+2*pi**2*alpha**4)*sigma*w+(((-8*pi**2)+32*pi-32)*alpha**7+(32*pi-16*pi**2)*alpha**5-8*pi**2*alpha**3)*exp((2*pi)/abs(alpha))*abs(alpha)*ksi*sigma)*exp(((alpha**2*ksi**2)/w**2)/2.0d+0)*x+((((-4*pi**2)+16*pi-16)*alpha**7+(16*pi-8*pi**2)*alpha**5-4*pi**2*alpha**3)*exp((2*pi)/abs(alpha))*abs(alpha)*sigma*w**2+((((-2*pi**4)+12*pi**3-24*pi**2+16*pi)*alpha**8+((-8*pi**4)+36*pi**3-48*pi**2+16*pi)*alpha**6+((-12*pi**4)+36*pi**3-24*pi**2)*alpha**4+(12*pi**3-8*pi**4)*alpha**2-2*pi**4)*abs(alpha)+((-2*pi**2)+8*pi-8)*alpha**8+(8*pi-4*pi**2)*alpha**6-2*pi**2*alpha**4)*ksi*sigma*w+((4*pi**2-16*pi+16)*alpha**7+(8*pi**2-16*pi)*alpha**5+4*pi**2*alpha**3)*exp((2*pi)/abs(alpha))*abs(alpha)*ksi**2*sigma)*exp(((alpha**2*ksi**2)/w**2)/2.0d+0))*exp(((alpha**2*x**2)/w**2)/2.0d+0+(ksi*x)/w**2)*erf((alpha*x-alpha*ksi)/(sqrt(2)*w))+(((4*pi**2-16*pi+16)*alpha**7+(8*pi**2-16*pi)*alpha**5+4*pi**2*alpha**3)*abs(alpha)*sigma*exp(((alpha**2*ksi**2)/w**2)/2.0d+0+(2*pi)/abs(alpha))*x**2+((((2*pi**4-12*pi**3+24*pi**2-16*pi)*alpha**8+(8*pi**4-36*pi**3+48*pi**2-16*pi)*alpha**6+(12*pi**4-36*pi**3+24*pi**2)*alpha**4+(8*pi**4-12*pi**3)*alpha**2+2*pi**4)*abs(alpha)+(2*pi**2-8*pi+8)*alpha**8+(4*pi**2-8*pi)*alpha**6+2*pi**2*alpha**4)*sigma*w+(((-8*pi**2)+32*pi-32)*alpha**7+(32*pi-16*pi**2)*alpha**5-8*pi**2*alpha**3)*exp((2*pi)/abs(alpha))*abs(alpha)*ksi*sigma)*exp(((alpha**2*ksi**2)/w**2)/2.0d+0)*x+((((-4*pi**2)+16*pi-16)*alpha**7+(16*pi-8*pi**2)*alpha**5-4*pi**2*alpha**3)*exp((2*pi)/abs(alpha))*abs(alpha)*sigma*w**2+((((-2*pi**4)+12*pi**3-24*pi**2+16*pi)*alpha**8+((-8*pi**4)+36*pi**3-48*pi**2+16*pi)*alpha**6+((-12*pi**4)+36*pi**3-24*pi**2)*alpha**4+(12*pi**3-8*pi**4)*alpha**2-2*pi**4)*abs(alpha)+((-2*pi**2)+8*pi-8)*alpha**8+(8*pi-4*pi**2)*alpha**6-2*pi**2*alpha**4)*ksi*sigma*w+((4*pi**2-16*pi+16)*alpha**7+(8*pi**2-16*pi)*alpha**5+4*pi**2*alpha**3)*exp((2*pi)/abs(alpha))*abs(alpha)*ksi**2*sigma)*exp(((alpha**2*ksi**2)/w**2)/2.0d+0))*exp(((alpha**2*x**2)/w**2)/2.0d+0+(ksi*x)/w**2))+((((-2**(5.0d+0/2.0d+0)*pi**2)+2**(9.0d+0/2.0d+0)*pi-2**(9.0d+0/2.0d+0))*alpha**8+(2**(9.0d+0/2.0d+0)*pi-2**(7.0d+0/2.0d+0)*pi**2)*alpha**6-2**(5.0d+0/2.0d+0)*pi**2*alpha**4)*exp((2*pi)/abs(alpha))*abs(alpha)*sigma*w*x+((((-2**(3.0d+0/2.0d+0)*pi**4)+3*2**(5.0d+0/2.0d+0)*pi**3-3*2**(7.0d+0/2.0d+0)*pi**2+2**(9.0d+0/2.0d+0)*pi)*alpha**9+((-2**(7.0d+0/2.0d+0)*pi**4)+9*2**(5.0d+0/2.0d+0)*pi**3-3*2**(9.0d+0/2.0d+0)*pi**2+2**(9.0d+0/2.0d+0)*pi)*alpha**7+((-3*2**(5.0d+0/2.0d+0)*pi**4)+9*2**(5.0d+0/2.0d+0)*pi**3-3*2**(7.0d+0/2.0d+0)*pi**2)*alpha**5+(3*2**(5.0d+0/2.0d+0)*pi**3-2**(7.0d+0/2.0d+0)*pi**4)*alpha**3-2**(3.0d+0/2.0d+0)*pi**4*alpha)*abs(alpha)+((-2**(3.0d+0/2.0d+0)*pi**2)+2**(7.0d+0/2.0d+0)*pi-2**(7.0d+0/2.0d+0))*alpha**9+(2**(7.0d+0/2.0d+0)*pi-2**(5.0d+0/2.0d+0)*pi**2)*alpha**7-2**(3.0d+0/2.0d+0)*pi**2*alpha**5)*sigma*w**2+((2**(5.0d+0/2.0d+0)*pi**2-2**(9.0d+0/2.0d+0)*pi+2**(9.0d+0/2.0d+0))*alpha**8+(2**(7.0d+0/2.0d+0)*pi**2-2**(9.0d+0/2.0d+0)*pi)*alpha**6+2**(5.0d+0/2.0d+0)*pi**2*alpha**4)*exp((2*pi)/abs(alpha))*abs(alpha)*ksi*sigma*w)*exp((alpha**2*ksi*x)/w**2+(ksi*x)/w**2))+((((-5*sqrt(2)*pi**3)+2**(9.0d+0/2.0d+0)*pi**2+2**(7.0d+0/2.0d+0)*pi)*alpha**8+((-3*2**(5.0d+0/2.0d+0)*pi**3)+9*2**(5.0d+0/2.0d+0)*pi**2+2**(7.0d+0/2.0d+0)*pi)*alpha**6+(5*2**(5.0d+0/2.0d+0)*pi**2-9*sqrt(2)*pi**3)*alpha**4-2**(3.0d+0/2.0d+0)*pi**3*alpha**2)*abs(alpha)*sigma*w*exp(((alpha**2*ksi**2)/w**2)/2.0d+0+(2*pi)/abs(alpha))*x+((5*sqrt(2)*pi**3-2**(9.0d+0/2.0d+0)*pi**2-2**(7.0d+0/2.0d+0)*pi)*alpha**8+(3*2**(5.0d+0/2.0d+0)*pi**3-9*2**(5.0d+0/2.0d+0)*pi**2-2**(7.0d+0/2.0d+0)*pi)*alpha**6+(9*sqrt(2)*pi**3-5*2**(5.0d+0/2.0d+0)*pi**2)*alpha**4+2**(3.0d+0/2.0d+0)*pi**3*alpha**2)*abs(alpha)*ksi*sigma*w*exp(((alpha**2*ksi**2)/w**2)/2.0d+0+(2*pi)/abs(alpha)))*exp(((alpha**2*x**2)/w**2)/2.0d+0+(ksi*x)/w**2)*erf((alpha*x-alpha*ksi)/(sqrt(2)*w))+((((-5*sqrt(2)*pi**3)+2**(9.0d+0/2.0d+0)*pi**2+2**(7.0d+0/2.0d+0)*pi)*alpha**8+((-3*2**(5.0d+0/2.0d+0)*pi**3)+9*2**(5.0d+0/2.0d+0)*pi**2+2**(7.0d+0/2.0d+0)*pi)*alpha**6+(5*2**(5.0d+0/2.0d+0)*pi**2-9*sqrt(2)*pi**3)*alpha**4-2**(3.0d+0/2.0d+0)*pi**3*alpha**2)*abs(alpha)*sigma*w*exp(((alpha**2*ksi**2)/w**2)/2.0d+0+(2*pi)/abs(alpha))*x+((5*sqrt(2)*pi**3-2**(9.0d+0/2.0d+0)*pi**2-2**(7.0d+0/2.0d+0)*pi)*alpha**8+(3*2**(5.0d+0/2.0d+0)*pi**3-9*2**(5.0d+0/2.0d+0)*pi**2-2**(7.0d+0/2.0d+0)*pi)*alpha**6+(9*sqrt(2)*pi**3-5*2**(5.0d+0/2.0d+0)*pi**2)*alpha**4+2**(3.0d+0/2.0d+0)*pi**3*alpha**2)*abs(alpha)*ksi*sigma*w*exp(((alpha**2*ksi**2)/w**2)/2.0d+0+(2*pi)/abs(alpha)))*exp(((alpha**2*x**2)/w**2)/2.0d+0+(ksi*x)/w**2)+sqrt(pi)*((10*pi**2-32*pi-16)*alpha**9+(24*pi**2-72*pi-16)*alpha**7+(18*pi**2-40*pi)*alpha**5+4*pi**2*alpha**3)*abs(alpha)*sigma*w**2*exp((alpha**2*ksi*x)/w**2+(ksi*x)/w**2+(2*pi)/abs(alpha)))+sqrt(alpha**2+1)*sqrt((pi-2)*alpha**2+pi)*(((2**(3.0d+0/2.0d+0)*pi**3-3*2**(5.0d+0/2.0d+0)*pi**2+3*2**(7.0d+0/2.0d+0)*pi-2**(9.0d+0/2.0d+0))*alpha**10+(2**(7.0d+0/2.0d+0)*pi**3-9*2**(5.0d+0/2.0d+0)*pi**2+3*2**(9.0d+0/2.0d+0)*pi-2**(9.0d+0/2.0d+0))*alpha**8+(3*2**(5.0d+0/2.0d+0)*pi**3-9*2**(5.0d+0/2.0d+0)*pi**2+3*2**(7.0d+0/2.0d+0)*pi)*alpha**6+(2**(7.0d+0/2.0d+0)*pi**3-3*2**(5.0d+0/2.0d+0)*pi**2)*alpha**4+2**(3.0d+0/2.0d+0)*pi**3*alpha**2)*exp((2*pi)/abs(alpha))*abs(alpha)*w**2*x+(((-2**(3.0d+0/2.0d+0)*pi**3)+3*2**(5.0d+0/2.0d+0)*pi**2-3*2**(7.0d+0/2.0d+0)*pi+2**(9.0d+0/2.0d+0))*alpha**10+((-2**(7.0d+0/2.0d+0)*pi**3)+9*2**(5.0d+0/2.0d+0)*pi**2-3*2**(9.0d+0/2.0d+0)*pi+2**(9.0d+0/2.0d+0))*alpha**8+((-3*2**(5.0d+0/2.0d+0)*pi**3)+9*2**(5.0d+0/2.0d+0)*pi**2-3*2**(7.0d+0/2.0d+0)*pi)*alpha**6+(3*2**(5.0d+0/2.0d+0)*pi**2-2**(7.0d+0/2.0d+0)*pi**3)*alpha**4-2**(3.0d+0/2.0d+0)*pi**3*alpha**2)*exp((2*pi)/abs(alpha))*abs(alpha)*ksi*w**2)*exp((alpha**2*ksi*x)/w**2+(ksi*x)/w**2)))/(sqrt(alpha**2+1)*sqrt((pi-2)*alpha**2+pi)*((2**(3.0d+0/2.0d+0)*pi**4-3*2**(5.0d+0/2.0d+0)*pi**3+3*2**(7.0d+0/2.0d+0)*pi**2-2**(9.0d+0/2.0d+0)*pi)*alpha**10+(2**(7.0d+0/2.0d+0)*pi**4-9*2**(5.0d+0/2.0d+0)*pi**3+3*2**(9.0d+0/2.0d+0)*pi**2-2**(9.0d+0/2.0d+0)*pi)*alpha**8+(3*2**(5.0d+0/2.0d+0)*pi**4-9*2**(5.0d+0/2.0d+0)*pi**3+3*2**(7.0d+0/2.0d+0)*pi**2)*alpha**6+(2**(7.0d+0/2.0d+0)*pi**4-3*2**(5.0d+0/2.0d+0)*pi**3)*alpha**4+2**(3.0d+0/2.0d+0)*pi**4*alpha**2)*abs(alpha)*w**4)
     end