Hilbert 希尔伯特矩阵

duan

用vc6 去求Hilbert矩阵，给不同的阶数 n

当 n 小于等于 27 时    A[0][0]=1;

n 大于 27 时 就开始变化
为什么？？？？？？？？？？？？？？？？？？？

1. #include<stdio.h>
   
2. int main()
   
3. {
   
4.         const int n=40;
   
5.         int i,j;
   
6.         double**A;
   
7.         double*b;
   
8.         A=new double*[n];
   
9.         for(i=0;i<n;i++)
   
10.                 A[i]=new double[n];
    
11.         b=new double[n];
    
12.         
    
13.         for(i=0;i<n;i++)
    
14.                 for(j=0;j<n;j++)
    
15.                         A[i][j]=1.0/(i+1+j+1-1);
    
16.         for(i=0;i<n;i++)
    
17.         {        
    
18.                 b[i]=0.0;
    
19.                 for(j=0;j<n;j++)
    
20.                         b[i]+=1.0/(i+1+j+1-1);               
    
21.         }
    
22.         //打印A
    
23.         printf("A=\n");
    
24.         for(i=0;i<n;i++)
    
25.         {
    
26.                 for(j=0;j<n;j++)
    
27.                         printf("%f\t",A[i][j]);
    
28. 
    
29.         }
    
30.         //打印b
    
31.         printf("b=\n");
    
32.         for(i=0;i<n;i++)
    
33.                         printf("%f\t",b[i]);
    
34.         return 0;
    
35. }

复制代码

vanentu

:sad鼓励一下楼主

vank

就是来顶 支持

逆战时代

顶。。。。。。。。。。。

dps521

一起研究一下呵呵

鱼C工作室.YCGZS

感谢分享

孤心傲

就是来顶 支持

孤心傲

就是来顶 支持