高斯消元程序为何无法运行帮忙看下,谢谢

高斯消元程序为何无法运行帮忙看下,谢谢,第1张

主元,要考虑系数绝对值的大小,而不只是数值大小。

试试这个:

/ ==========================================

GaoSi_XiaoQuc

==========================================/

#include <stdioh>

#include <stdlibh>

#include <mathh>

float ColPivot(float c,int n);

void main()

{

int i,j;

float x;

int n=4;

float a[4][5]={4,-1,0,2,-10,-1,4,-1,0,-70,0,-1,4,-1,90,2,0,-1,40,00};

x=ColPivot(a[0],n);

for (i=0;i<n;i++) printf("x[%d]=%f\n",i,x[i]);

getch();

}

float ColPivot(float c,int n)

{

int i,j,t,k;

float x,p;

x=(float )malloc(nsizeof(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!=i)

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;

}

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;

}

SAP2000是国内外使用最广泛的通用有限元结构分析程序。SAP2000程序有别于其它一般结构有限元程序的最大特点就在于它的强大的分析功能。SAP2000中使用许多不同类型的分析,它基本上集成了现有结构分析中经常遇到方法,如时程分析、地震动输入、动力分析以及Push-over分析等等。另外还包括:静力分析、用特征向量或Ritz向量进行振动模式的模态分析、对地震反应的反应谱分析等等。这些不同类型的分析可在程序的同一次运行中进行,并把结果综合起来输出。从简单的二维框架静力分析到复杂的三维非线性动力分析,SAP2000能为所有结构分析和设计提供了解决方案。

为了让广大网友掌握SAP2000这种强大的分析工具,也为了推广这一著名的通用结构分析软件,我们在这里为广大网友推荐--如何学习SAP2000系列书籍,希望各位朋友在结构成长的道路上少走弯路,快速进步,早日成为合格的结构工程师!

取证、依据、陈述、决定、处罚。石臼湖管理站行政处罚二千元程序是1、调查取证,2、告知处罚事实、理由、依据和有关权利,3、听取陈述、申辩或者举行听证会,4、作出行政处罚决定,5、作出行政处罚。

以上就是关于高斯消元程序为何无法运行帮忙看下,谢谢全部的内容,包括:高斯消元程序为何无法运行帮忙看下,谢谢、如何学习SAP2000系列、石臼湖管理站行政处罚二千元程序是什么等相关内容解答,如果想了解更多相关内容,可以关注我们,你们的支持是我们更新的动力!

欢迎分享,转载请注明来源:内存溢出

原文地址: http://outofmemory.cn/zz/9321881.html

(0)
打赏 微信扫一扫 微信扫一扫 支付宝扫一扫 支付宝扫一扫
上一篇 2023-04-27
下一篇 2023-04-27

发表评论

登录后才能评论

评论列表(0条)

保存