Функция Грина

double grinn(vector_t const& x, vector_t const& y) {
   return ( 1. / ( 4. * pi * abs( x - y ) ) )
        - ( 1. / ( 4. * pi * abs( vector_t( x.x() - y.x(), x.y() - y.y(), x.z() + y.z() ) ) ) );
}

double Grin(double* x, double* y){
    double* yy;
           double* xy, *xyy;
    yy[2] = -y[2];
           
double t, z;
    for(int i = 0; i < 3; i++){
        t += pow(xy[i], 2);
        z += pow(xyy[i], 2);
    }
double mod1 = sqrt(t);
double mod2 = sqrt(z);
 
    return pow(4 * M_PI * mod1, -1) - pow(4 * M_PI * mod2, -1);
    
}

#include <vector>

struct vector_t {
 
   double x, y, z;
};

#define _USE_MATH_DEFINES
#include <iostream>
#include <cmath>
 
const int size = 3;
 
double vmod( double* v )
{
    double res = 0;
    for( int i = 0; i<size; ++i )
    {
        res += v[i]*v[i];
    }
    return sqrt( res );
}
 
void vsub( double* r, double* a, double* b )
{
    for( int i=0; i<size; ++i)
    {
        r[i] = a[i] - b[i];
    }
}
 
double grin( double* x, double* y, double* ys )
{
    double temp[size];
 
    vsub( temp, x, y );
    double a1 = 1/vmod( temp );
 
    vsub( temp, x, ys );
    double a2 = 1/vmod( temp );
 
    return (a1-a2)/(4*M_PI);
}
 
int main()
{
    double x[size] = { 1,2,3 };
    double y[size] = { 4,5,6 };
    double ys[size] = { y[0], y[1], -y[2] };
 
    std::cout << grin( x, y, ys ) << std::endl;
 
    std::system( "pause" );
 
    return 0;
}