Hilbert 希尔伯特矩阵

duan · 发表于 2015-4-9 20:59:09

马上注册，结交更多好友，享用更多功能^_^

您需要登录才可以下载或查看，没有账号？立即注册

x

本帖最后由 duan 于 2015-4-9 21:01 编辑

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

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

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

#include<stdio.h>
int main()
{
        const int n=40;
        int i,j;
        double**A;
        double*b;
        A=new double*[n];
        for(i=0;i<n;i++)
                A[i]=new double[n];
        b=new double[n];
        
        for(i=0;i<n;i++)
                for(j=0;j<n;j++)
                        A[i][j]=1.0/(i+1+j+1-1);
        for(i=0;i<n;i++)
        {        
                b[i]=0.0;
                for(j=0;j<n;j++)
                        b[i]+=1.0/(i+1+j+1-1);                
        }
        //打印A
        printf("A=\n");
        for(i=0;i<n;i++)
        {
                for(j=0;j<n;j++)
                        printf("%f\t",A[i][j]);

        }
        //打印b
        printf("b=\n");
        for(i=0;i<n;i++)
                        printf("%f\t",b[i]);
        return 0;
}

vanentu · 发表于 2015-5-26 01:29:37

:sad鼓励一下楼主

vank · 发表于 2015-5-27 05:23:58

就是来顶支持

逆战时代 · 发表于 2015-6-5 15:20:30

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

dps521 · 发表于 2015-11-21 16:06:00

一起研究一下呵呵

鱼C工作室.YCGZS · 发表于 2015-11-24 16:49:02

感谢分享

孤心傲 · 发表于 2015-11-30 12:58:51

就是来顶支持

孤心傲 · 发表于 2015-11-30 12:59:26

就是来顶支持

账号		自动登录	找回密码
密码			立即注册