| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 617 人关注过本帖
标题:不知我的这个程序哪里不对,算的结果和手算不一样了,请看看,谢谢
只看楼主 加入收藏
sunstudent
Rank: 1
等 级:新手上路
帖 子:110
专家分:0
注 册:2007-10-14
结帖率:100%
收藏
 问题点数:0 回复次数:1 
不知我的这个程序哪里不对,算的结果和手算不一样了,请看看,谢谢
#include "stdio.h"
#include "math.h"
#include <malloc.h>
float *gaoshi(float*c,int n)
{
    int i,j,k,t;
    float *x,p;
    x=(float*)malloc(n*sizeof(float));
    for (i=0;i<=n-2;i++)
    {
            k=i;
        for (j=i+1;j<=n-1;j++)

            if(fabs(*(c+j*(n+1)+i))>fabs(*(c+k*(n+1)+i)))
            {    k=j;
            if(k!=j)
                for (j=i;j<=n;j++)
                {
                    p=*(c+i*(n+1)+j);
*(c+i*(n+1)+j)=*(c+k*(n+1)+j);
*(c+k*(n+1)+j)=p;
                }
                }
            else
            {
            
    for (j=i+1;j<=n-1;j++)
    {
        p=(*(c+j*(n+1)+i))/(*(c+i*(n+1)+i));
        for(t=i;t<=n;t++)
            *(c+j*(n+1)+t)-=p*(*(c+i*(n+1)+t));
        
    }
            }
    }
    
    for (i=n-1;i>=0;i--)
    {
        for(j=n-1;j>=i+1;j--)
        {    (*(c+i*(n+1)+n))-=x[j]*(*(c+i*(n+1)+j));}
        x[i]=*(c+i*(n+1)+n)/(*(c+i*(n+1)*i));
        
    }
    return x;
}

void main()
{
    float a[4][5]={1,1,0,3,4,
                   2,1,-1,1,1,
                   3,-1,-1,3,-3,
                   -1,2,3,-1,4};
    
float *x;
//x=a[0];
int i;
x=gaoshi(a[0],4);
for (i=0;i<=3;i++)
printf(" x[%d]=%f\n",i+1,x[i]);
    }我要算的是个方程
x1+x2+3*x4=4
2*x1+x2-x3+x4=1
3*x1-x2-x3+3*x4=-3
-x1+2*x2+3*x3-x4=4我是用高斯列主元素消去法解的,觉的没什么错,但就是不对,请帮助看看,谢谢
搜索更多相关主题的帖子: float include 结果 
2007-12-15 11:21
sunstudent
Rank: 1
等 级:新手上路
帖 子:110
专家分:0
注 册:2007-10-14
收藏
得分:0 
谢谢大家了,不用了
我知哪里不对了
2007-12-15 13:09
快速回复:不知我的这个程序哪里不对,算的结果和手算不一样了,请看看,谢谢
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.016080 second(s), 7 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved