F(1,1) = 1
F(1,2) = -((-((alpha*exp(-(2*%pi)/abs(alpha)))/abs(alpha))/2.0d+0)-(%pi**((-3.0d+0)/2.0d+0)*(4-%pi)*alpha**3*(alpha**2+1)**((-3.0d+0)/2.0d+0))/(sqrt(2)*(1-(2*alpha**2)/(%pi*(alpha**2+1)))**2.0d+0)+(sqrt(2)*alpha)/(sqrt(%pi)*sqrt(alpha**2+1)))/sqrt(1-(2*alpha**2)/(%pi*(alpha**2+1)))
F(1,3) = (((4*alpha**3)/(%pi*(alpha**2+1)**2)-(4*alpha)/(%pi*(alpha**2+1)))*(1-(2*alpha**2)/(%pi*(alpha**2+1)))**((-3.0d+0)/2.0d+0)*((-((alpha*exp(-(2*%pi)/abs(alpha)))/abs(alpha))/2.0d+0)-(%pi**((-3.0d+0)/2.0d+0)*(4-%pi)*alpha**3*(alpha**2+1)**((-3.0d+0)/2.0d+0))/(sqrt(2)*(1-(2*alpha**2)/(%pi*(alpha**2+1)))**2.0d+0)+(sqrt(2)*alpha)/(sqrt(%pi)*sqrt(alpha**2+1)))*sigma)/2.0d+0-(((-(%pi*exp(-(2*%pi)/abs(alpha)))/alpha**2)-(3*%pi**((-3.0d+0)/2.0d+0)*(4-%pi)*alpha**2*(alpha**2+1)**((-3.0d+0)/2.0d+0))/(sqrt(2)*(1-(2*alpha**2)/(%pi*(alpha**2+1)))**2.0d+0)+(3*%pi**((-3.0d+0)/2.0d+0)*(4-%pi)*alpha**4*(alpha**2+1)**((-5.0d+0)/2.0d+0))/(sqrt(2)*(1-(2*alpha**2)/(%pi*(alpha**2+1)))**2.0d+0)+(2.0d+0*%pi**((-3.0d+0)/2.0d+0)*(4-%pi)*alpha**3*(alpha**2+1)**((-3.0d+0)/2.0d+0)*((4*alpha**3)/(%pi*(alpha**2+1)**2)-(4*alpha)/(%pi*(alpha**2+1))))/(sqrt(2)*(1-(2*alpha**2)/(%pi*(alpha**2+1)))**3.0d+0)+sqrt(2)/(sqrt(%pi)*sqrt(alpha**2+1))-(sqrt(2)*alpha**2*(alpha**2+1)**((-3.0d+0)/2.0d+0))/sqrt(%pi))*sigma)/sqrt(1-(2*alpha**2)/(%pi*(alpha**2+1)))
F(2,1) = 0
F(2,2) = 1/sqrt(1-(2*alpha**2)/(%pi*(alpha**2+1)))
F(2,3) = -(((4*alpha**3)/(%pi*(alpha**2+1)**2)-(4*alpha)/(%pi*(alpha**2+1)))*(1-(2*alpha**2)/(%pi*(alpha**2+1)))**((-3.0d+0)/2.0d+0)*sigma)/2.0d+0
F(3,1) = 0
F(3,2) = 0
F(3,3) = 1