| 查看: 1742 | 回复: 9 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】为什么我这个程序编译都没有问题，就是运行不出来？【已完结】 已有5人参与
|
|||
|
#include #include #include using namespace std; class sdet { private: int n; double **a,det; public: sdet(int nn) { int i; n=nn; a=new double * [n]; for(i=0;i } void input(); void a_det(); void output(); sdet() { int i; for(i=0;i delete[]a; } }; void sdet::input() { int i,j; char str1[20]; cout<<"\n 输入文件名:"; cin>>str1; ifstream fin(str1); if(!fin) {cout<<"\n不能打开这个文件"< fin.close(); } void sdet::a_det() { int i,j,k,is,js; double f,q,d; f=1.0;det=1.0; for(k=0;k<=n-2;k++) { q=0.0; for(i=k;i<=n-1;i++) for(j=k;j<=n-1;j++) { d=fabs(a[j]); if(d>q){q=d;is=i;js=j;} } if(q+1.0==1.0)det=0.0; else { if(is!=k) { f=-f; for(j=k;j<=n-1;j++) { d=a[k][j];a[k][j]=a[is][j];a[is][j]=d; } } if(js!=k) { f=-f; for(j=k;j<=n-1;j++) { d=a[js];a[js]=a[k];a[k]=d; } } det=det*a[k][k]; for(i=k+1;i<=n-1;i++) { d=a[k]/a[k][k]; for(j=k+1;j<=n-1;j++) { a[j]=a[j]-d*a[k][j]; } } } } det=f*det*a[n-1][n-1]; } void sdet:  utput(){ char str2[20]; cout<<"\n输入文件名:"; cin>>str2; ofstream fout (str2); if(!fout) {cout<<"\n不能打开这个文件"< } void main() { sdet d(4); d.input(); d.a_det(); d.output(); } 文件是这样的： 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 [ Last edited by nono2009 on 2010-6-28 at 12:56 ] |
回复此楼» 猜你喜欢
今天系统多次维护，明天很可能放榜！
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有74人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
zhangeil04
铁杆木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 7928
- 红花: 1
- 帖子: 1939
- 在线: 25.6小时
- 虫号: 729893
- 注册: 2009-03-24
- 专业: 物理学I
9楼2010-06-27 21:33:58
2楼2010-06-25 01:09:12
3楼2010-06-25 08:41:31
alexissp
金虫 (正式写手)
干活的
- 应助: 0 (幼儿园)
- 金币: 1428.3
- 帖子: 670
- 在线: 14.9小时
- 虫号: 478745
- 注册: 2007-12-14
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2010-06-25 09:53:56
25