| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 835 人关注过本帖
标题:为什么手算的和编程算,不一样了,请帮助看看,程序哪里错了,谢谢
只看楼主 加入收藏
sunstudent
Rank: 1
等 级:新手上路
帖 子:110
专家分:0
注 册:2007-10-14
结帖率:100%
收藏
 问题点数:0 回复次数:3 
为什么手算的和编程算,不一样了,请帮助看看,程序哪里错了,谢谢
#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我是用高斯列主元素消去法解的,觉的没什么错,但就是不对,请帮助看看,谢谢
搜索更多相关主题的帖子: include 
2007-12-15 12:33
sunstudent
Rank: 1
等 级:新手上路
帖 子:110
专家分:0
注 册:2007-10-14
收藏
得分:0 
谢谢了,不用了
我知哪里不对了
加号变乘号了
谢谢
2007-12-15 13:10
StarWing83
Rank: 8Rank: 8
来 自:仙女座大星云
等 级:贵宾
威 望:19
帖 子:3951
专家分:748
注 册:2007-11-16
收藏
得分:0 
顺便说说,也可以使用矩阵求值的方法计算,可能会简单一点……
还有,楼主要是解决了的话麻烦编辑一下标题,加个“(已解决)”

专心编程………
飞燕算法初级群:3996098
我的Blog
2007-12-16 06:23
sunstudent
Rank: 1
等 级:新手上路
帖 子:110
专家分:0
注 册:2007-10-14
收藏
得分:0 
过了时间
好像就不能改了
不好意思
2007-12-21 15:50
快速回复:为什么手算的和编程算,不一样了,请帮助看看,程序哪里错了,谢谢
数据加载中...
 
   



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

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