Ошибки при компиляции

struct XYZ
{
  double x, y, z;
}
 
double sqr(double x) {return x * x; }
 
void CalcNormals(XYZ a, XYZ b, XYZ c, XYZ& n)
{
  double wrki;
  XYZ v1, v2;
 
  v1.x = a.x - b.x;
  v1.y = a.y - b.y;
  v1.z = a.z - b.z;
 
  v2.x = b.x - c.x;
  v2.y = b.y - c.y;
  v2.z = b.z - c.z;
 
  wrki = sqrt(sqr(v1.y*v2.z - v1.z * v2.y) + sqr(v1.z * v2.x - v1.x * v2.z) + sqr(v1.x * v2.y - v1.y * v2.x));
  n.x = (v1.y * v2.z - v1.z * v2.y) / wrki;
  n.y = (v1.z * v2.x - v1.x * v2.z) / wrki;
  n.z = (v1.x * v2.y - v1.y * v2.x) / wrki;
}

struct XYZ;

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

class Matrix
{
 
    int n,m;
    double **array;
 
public:
    Matrix(Matrix & A);//--конструктор копирования
    ~Matrix(void);//---деструктор
    void OutMatrix(void);//---вывод матрицы на экран
    void EnterMatrix(void);//---ввод матрицы
    void SerchAkstrem(void); //---поиск эксремумов и их строки
    void Polindrom(void);  //---поиск чисел полиндромов
    //-----перегрузка операторов----
    Matrix&operator=(Matrix &a); //---оператор =
    Matrix operator+(Matrix &a);//---оператор +
    Matrix operator-(Matrix &a);//---оператор -
    Matrix operator*(double a);//---оператор *a
    Matrix operator/(double a);//---оператор /a
    Matrix operator*(Matrix & a);//---оператор Matrix*Matrix
    Matrix Transponiring(void);//---транспонирование матрицы
    double SkPrOnVek(Matrix b);//---сколярное произведение 2х векторов doublel Matrix1(Matrix2)
    Matrix operator-(void);//---внесение минуса в матрицу (*(-1))
    double Determinant3x3max(void);//---детерминант матрицы. 3х3-максимальный рамер матрицы
};

Matrix q;
    q.EnterMatrix(); // здесь вызывается  метод , в частности ввода матрицы.

XYZ CalcNormals(XYZ a, XYZ b, XYZ c)
{
  double wrki;
  XYZ v1, v2;
 
  v1.x = a.x - b.x;
  v1.y = a.y - b.y;
  v1.z = a.z - b.z;
 
  v2.x = b.x - c.x;
  v2.y = b.y - c.y;
  v2.z = b.z - c.z;
 XYZ n;
  wrki = sqrt(sqr(v1.y*v2.z - v1.z * v2.y) + sqr(v1.z * v2.x - v1.x * v2.z) + sqr(v1.x * v2.y - v1.y * v2.x));
  n.x = (v1.y * v2.z - v1.z * v2.y) / wrki;
  n.y = (v1.z * v2.x - v1.x * v2.z) / wrki;
  n.z = (v1.x * v2.y - v1.y * v2.x) / wrki;
return n;
}

void main()
{
    XYZ x1,x2,x3;
x1=CalcNormals(/*тут 3 точки*/);
x2=CalcNormals(/*тут 3 точки*/);
x3=CalcNormals(/*тут 3 точки*/);
 
}

x3=CalcNormals(/*тут 3 точки*/);