| 查看: 725 | 回复: 1 | ||
[求助]
求帮助！改Function中的可变系数进行循环计算的问题
|
|
没有金币了啊！！！希望发扬雷锋精神啊！求大神帮忙解决问题啊！ 我想用Fortran 求解非线性方程组，程序是Monte Carlo法 程序由主程序，Function, Subroutine三部分组成。 Function给出了非线性方程组，通过subroutine进行计算非线性方程组，主程序给出了初始计算值X1,X2,X3。 问题是， Function 中有个系数M是变化的， 在subroutine（或者主程序）中，先设定了初始值M0，每次计算结束了X后，由X计算新的M=f(x). 如果新Mc 和旧M差别在1.0E-05内就结束计算， 否则重新赋予Function 系数M, 重新调用子程序循环计算X. DIMENSION X(N), Y( N) REAL NRND1 External F Double precision R, M, M0, Mc Double precision R=1.0E-4; d=30.0E-3; k=383.0; * h=40000.0; T=303.0 A= B K= 1 R= 1.0E0 M0=450.0 40 M=M0 Z= F(X, N, M) 10 IF (A. GE. EPS) THEN L= L+ 1 DO 20 I= 1, N 20 Y(I) = -A+ 2.0* A* NRND1(R) + X(I) Z1= F(Y, N, M) K= K+ 1 IF ( Z1. GE. Z) THEN IF (K .GT. M) THEN K= 1 A= A/2.0 ENDIF GOTO 10 ELSE K= 1 DO 30 I= 1, N 30 X(I) = Y( I) Z= Z1 IF ( Z. GE. EPS) GOTO 10 ENDIF ENDIF Mc=T+(1000.0*X(1)-M)*(d/k+1/h)/R IF(ABS(M-Mc) .GE. EPS) THEN M0=(M+Mc)/2 GOTO 40 ENDIF END 上面程序是对subroutine的修改。倒数第6行直接调用X1对新的Mc计算，并赋值给M0，然后返回行40。 但是貌似根本就没有把新的M传递给方程F(X,N,M)进行循环计算。最后只有第一次的结果。 难道计算得到的X1不能直接调用，用于方程计算么？(尝试过把X1直接调出主程序计算倒数第6行的方程，方程结果根本就不对) 此子程序有否错误呢？ 谢谢！ |
回复此楼» 猜你喜欢
大佬帮着分析今年的评审意见？该课题连续申请了三年了，明年要不要换个方向。
已经有11人回复
-
两块石头
已经有10人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有6人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】怎么用spss计算两个变量间的相关系数
已经有5人回复
- 利用matlab求一阶变系数微分方程组的数值解。
已经有7人回复
- 求助！！审稿人提出的关于检测限的问题，不懂是什么意思，求帮助。。。
已经有14人回复
- 变系数偏微分方程matlab求解
已经有6人回复
- 求回归方程参数的最优化算法matlab代码
已经有17人回复
- 用1stOpt进行复数曲线拟合时，结果出错。急求教！！！
已经有5人回复
- 1stopt 拟合隐函数的四参数
已经有7人回复
- 麻烦大家帮忙看看这个审稿人什么意思（关于数据分析和”化学计量系数？“）
已经有5人回复
- 哪位虫友帮帮忙，我的comsol计算时出现如下错误
已经有9人回复
- 两条不同的氮气吸附脱附曲线，求帮助~~
已经有18人回复
- 关于用化学吸附仪做化学脉冲实验计算金属分散度的一些问题！求帮助！！
已经有11人回复
- 【求助】FLUENT UDF 气液两相流 汽蚀 汽化潜热(问题详见内容)
已经有21人回复
- 请问，参比制剂溶出曲线的变异系数偏大，是怎么回事啊？
已经有9人回复
- matlab中系数是矩阵的积分
已经有5人回复
- 如何计算变位齿轮的齿廓重叠干涉系数
已经有4人回复
- 【求助】迎风角变化下对长方体的风阻系数的关系，请求帮助！
已经有7人回复
- 【求助/交流】生物统计学问题，变异系数
已经有11人回复
- 【求助】MUSIC的MD模块计算甲烷在IRMOF-1上的扩散系数
已经有11人回复
- 【求助】MATLAB计算NRTL活度系数
已经有7人回复
- 【求助】常微分方程组求解中系数与某变量值关联的问题
已经有12人回复
- 【求助】如何从粒度分布曲线上求算变异系数CV值？
已经有4人回复
学虫dj
金虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 892.4
- 散金: 63
- 帖子: 688
- 在线: 89.9小时
- 虫号: 1829966
- 注册: 2012-05-22
- 性别: GG
- 专业: 应用地球物理学
2楼2014-05-26 16:40:35
100