return (z*z*cos(phi/2)*cos(phi/2)-z*dy*sin(phi)+(dy*dy-y*y)*sin(phi/2)*sin(phi/2))*exp(-1/(2*a)*(sqrt((x+dx)*(x+dx)+(y+dy)*(y+dy)+z*z)+sqrt((x-dx)*(x-dx)+(y-dy)*(y-dy)+z*z)));
Переписал вашу "Подинтегральную функцию".
Это она?
1. Проверьте это верно записано?
2. И аргументы под радикалами,
корни у вас должны возвращаться толко положительные или плюс/минус?
PS.Проверить значит - не в программе, а в физической формуле.
Например, в программе машина вычисляет корень только положительный.
А в физической формуле как надо ?
Надеюсь вы меня поняли.