矩阵的相关操作
if(kk%n==0) cout<<endl;
}
cout<<endl;
return 0;
}
void temp(double aa[],double bb[],int n)
{ int i;
double temp1;
for(i=0;i<n;i++)
{temp1=aa[i];aa[i]=bb[i];bb[i]=temp1; }
}
double fun(double array[n][n])
{ int ii,jj,k,u;
double det1=1,yin;
for(ii=0;ii<n;ii++)
{ if(array[ii][ii]==0)
for(jj=ii;jj<n;jj++)
{ if(array[jj][ii]!=0)
temp(array[ii],array[jj],n);
}
for(k=ii+1;k<n;k++)
{ yin=-1*array[k][ii]/array[ii][ii];
for(u=0;u<n;u++)
{ array[k][u]=array[k][u]+array[ii][u]*yin;
}
}
}
for(ii=0;ii<n;ii++)
det1=det1*array[ii][ii];
return (det1);
}
6.矩阵求逆的程序代码(此时矩阵也必须为方阵):
#include<iostream>
#include<iomanip>
using namespace std; //注意在程序运行过程中那些量是变化的。
int const n=3;
int main()
{
void temp(double aa[],double bb[],int n);
百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读,免费范文网,提供经典小说教育文库矩阵代码大全(矩阵的逆、乘法、加法、行列式)(c++程序)(4)在线全文阅读。
相关推荐: