行列式求值

来源:互联网 发布:电路模拟软件下载 编辑:程序博客网 时间:2024/04/30 15:04

#include <iostream>
using namespace std;
void main()
{
float **a,result;
int n,i;
float hls(float **a,int n);
cout<<"输入行列式的阶数:";
cin>>n;
a=new float *[n];
for(i=0;i<n;i++)
  a=new float[n];
cout<<"输入初始值:"<<endl;
for(i=0;i<n;i++)
  for(int j=0;j<n;j++)
   cin>>a[j];
result=hls(a,n);
cout<<"运算结果是:"<<result<<endl;
}
float hls(float **a,int n)
{
int i,j,k;
float result=1;
for(i=0;i<n-1;i++)
  for(j=i+1;j<n;j++)
   for(k=n;k>=0;k--)
    a[j][k]=a[j][k]-a[k]*a[j]/a;
for(i=0;i<n;i++)
  result=result*a;
return result;
}

原创粉丝点击