J(1,1) = ((x-ksi)*exp(-((x-ksi)**2/w**2)/2.0d+0)*(erf((alpha*(x-ksi))/(sqrt(2)*w))+1))/(sqrt(2)*sqrt(%pi)*w**3)-(alpha*exp((-((alpha**2*(x-ksi)**2)/w**2)/2.0d+0)-((x-ksi)**2/w**2)/2.0d+0))/(%pi*w**2)
J(1,2) = ((x-ksi)**2*exp(-((x-ksi)**2/w**2)/2.0d+0)*(erf((alpha*(x-ksi))/(sqrt(2)*w))+1))/(sqrt(2)*sqrt(%pi)*w**4)-(exp(-((x-ksi)**2/w**2)/2.0d+0)*(erf((alpha*(x-ksi))/(sqrt(2)*w))+1))/(sqrt(2)*sqrt(%pi)*w**2)-(alpha*(x-ksi)*exp((-((alpha**2*(x-ksi)**2)/w**2)/2.0d+0)-((x-ksi)**2/w**2)/2.0d+0))/(%pi*w**3)
J(1,3) = ((x-ksi)*exp((-((alpha**2*(x-ksi)**2)/w**2)/2.0d+0)-((x-ksi)**2/w**2)/2.0d+0))/(%pi*w**2)