Реферат: Разработка библиотечных средств

if (n==1)

return (*this)(0,0);

for(int i=0; i<m; i++)

if ((*this)(i,0))

{

static Matrix<VARTYPE> M;

M <<= *this;

VARTYPE d=M(i,0)*(i%2?-1:1);

return d*M.Gauss(i,0).minor(i,0).determ();

}

return 0.0;

}

template <class VARTYPE>

VARTYPE Matrix<VARTYPE>::vmodul() // Модуль вектора

{

VARTYPE d=0;

if (n!=1) ERROR_MATRIX(9);

static Matrix<VARTYPE> M;

M <<= *this;

for(int i=0; i<m; i++)

d=d+M(i,0)*M(i,0);

return sqrt(d);

}

template <class VARTYPE>

Matrix<VARTYPE>& Matrix<VARTYPE>::Gauss(dim M, dim N)

{

Matrix<VARTYPE>& A=*this;

if (!A(M,N)) ERROR_MATRIX(5);

for(int i=0,j;i<m;i++)