program test
      real der(3)
      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-2
     &4*pi+16)*alpha**5+(10*pi**2-24*pi)*alpha**3+4*pi**2*alpha)*abs(alp
     &ha)*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*alp
     &ha)*abs(alpha)*ksi*w*exp(((alpha**2*ksi**2)/w**2)/2.0d+0+(2*pi)/ab
     &s(alpha)))*exp(((alpha**2*x**2)/w**2)/2.0d+0+(ksi*x)/w**2)*erf((al
     &pha*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)*alph
     &a**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)*al
     &pha**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.0
     &d+0)*pi)*alpha**6+(3*2**(5.0d+0/2.0d+0)*pi**3-2**(1.10000000000000
     &01d+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)*p
     &i**3-2**(5.0d+0/2.0d+0)*pi**2+2**(5.0d+0/2.0d+0)*pi)*alpha**7+(3*s
     &qrt(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(((a
     &lpha**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)*alp
     &ha**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)*p
     &i)*alpha**5+(2**(5.0d+0/2.0d+0)*pi**2-3*sqrt(2)*pi**3)*alpha**3-sq
     &rt(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)*al
     &pha**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)*ks
     &i)*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(al
     &pha)*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.0
     &d+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(a
     &lpha)*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)*al
     &pha**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)/a
     &bs(alpha))*x**2+((((2*pi**4-12*pi**3+24*pi**2-16*pi)*alpha**8+(8*p
     &i**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(((al
     &pha**2*ksi**2)/w**2)/2.0d+0)*x+((((-4*pi**2)+16*pi-16)*alpha**7+(1
     &6*pi-8*pi**2)*alpha**5-4*pi**2*alpha**3)*exp((2*pi)/abs(alpha))*ab
     &s(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(al
     &pha)+((-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*al
     &pha**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*p
     &i**2-8*pi+8)*alpha**8+(4*pi**2-8*pi)*alpha**6+2*pi**2*alpha**4)*si
     &gma*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)*alph
     &a**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)*k
     &si**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)*ab
     &s(alpha)*ksi*sigma*w*exp(((alpha**2*ksi**2)/w**2)/2.0d+0+(2*pi)/ab
     &s(alpha)))*exp(((alpha**2*x**2)/w**2)/2.0d+0+(ksi*x)/w**2)*erf((al
     &pha*x-alpha*ksi)/(sqrt(2)*w))+((((-5*sqrt(2)*pi**3)+2**(9.0d+0/2.0
     &d+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.0
     &d+0)*pi**2-2**(7.0d+0/2.0d+0)*pi)*alpha**8+(3*2**(5.0d+0/2.0d+0)*p
     &i**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.0
     &d+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-7
     &2*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)/ab
     &s(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(alp
     &ha))*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)*al
     &pha**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)*a
     &bs(alpha)*w**4)
      end