Предыдущая страница реферата | 2  3  4  5  6  7  8  9  10  11  12 | Следующая страница реферата



A(j,i)=(*this)(i,j);

return A;

}

template <class VARTYPE>

VARTYPE Matrix<VARTYPE>::determ() // рекурсивно находит определитель матрицы

{

if (n!=m) ERROR_MATRIX(4);

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++)

for(j=0;j<n;j++)

if (i!=M && j!=N)

A(i,j)-=A(M,j)*A(i,N)/A(M,N);

for(j=0;j<n;j++)

if (j!=N)

A(M,j)/=A(M,N);

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

A(i,N)=0;

A(M,N)=1;

Рекомендуем скачать другие рефераты по теме: инновационная деятельность, решебник по английскому языку.

Предыдущая страница реферата | 2  3  4  5  6  7  8  9  10  11  12 | Следующая страница реферата

Поделитесь этой записью или добавьте в закладки