| 查看: 816 | 回复: 4 | ||
[求助]
各位帮我看下这个simpson积分怎么理解啊
|
|
SIMPSON INTEGRATION C NTOTH = 1000 C A = P(IC) C H = A/(TWO*DFLOAT(NTOTH)) C S = ONE/TWO*FUN(A,QSAT,PK,TN)C SIMPSONLOOP: DO I=1,NTOTH-1 C S = S+TWO*FUN((TWO*DFLOAT(I)-ONE)*H,QSAT,PK,TN) C + +FUN(TWO*DFLOAT(I)*H,QSAT,PK,TN)C S = S+TWO*FUN((TWO*DFLOAT(NTOTH)-ONE)*H,QSAT,PK,TN) C ENDDO SIMPSONLOOP C H = A/DFLOAT(NTOTH) C NTOTH1 = NTOTH + 1 红色字描述 P相当于x,那么第一个循环相当于x0,x1,x2,算得的值，我想问一下x0对应的值为什么是这个表达式，还有为什么会少了1/3啊 |
回复此楼» 猜你喜欢
药学105500专硕289分求调剂
已经有3人回复
-
0703总分331求调剂
已经有15人回复
-
一志愿北京2，材料与化工308求调剂
已经有13人回复
-
318求调剂
已经有14人回复
-
285求调剂
已经有10人回复
-
312求调剂
已经有11人回复
-
273求调剂
已经有6人回复
-
生物学308分求调剂（一志愿华东师大）
已经有8人回复
-
生物学求调剂 一志愿沪9，326分
已经有3人回复
-
080500求调剂
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 麻烦大家帮我看一下，想用软件预测2015年的垃圾组成，应该怎么做
已经有5人回复
- 我的质粒怎么回事，一直不好，请各位达人帮忙看一下
已经有21人回复
- 各位大侠帮忙看看我的这个液晶小分子的DSC和TGA图，求解答。
已经有14人回复
- 我编的Simpson积分法fortran程序给不出结果，大侠们看看哪里出了问题？
已经有4人回复
- 请各位大虫帮我看看啊，核磁做H1NMR的时候会不会出现少一个H的情况啊
已经有8人回复
- 请各位大侠帮忙看一下我的液相色谱图究竟怎么了？
已经有13人回复
- 【求助】请帮忙看一下这个偶合常数怎么计算
已经有10人回复
- 大家帮我看看这个图怎么做啊
已经有9人回复
- 投稿被拒了，不甘心想写信给编辑，大家帮我看下该怎么写？
已经有10人回复
- 大家帮我看看这个图是怎么画出来的
已经有12人回复
- 请各位投稿达人帮忙看一下这个审稿意见
已经有7人回复
- 【求助】请大家帮忙看一下这三个题（微积分方面）【已解决】
已经有11人回复
- 【求助】大家帮我看下这段话具体是什么意思，我每句能看懂，但整体不会理解
已经有14人回复
snoopyzhao
至尊木虫 (职业作家)
- 程序强帖: 16
- 应助: 157 (高中生)
- 贵宾: 0.02
- 金币: 18844.7
- 红花: 29
- 帖子: 3803
- 在线: 1422.4小时
- 虫号: 183750
- 注册: 2006-02-13
- 专业: 污染生态化学
2楼2011-12-15 18:34:28
3楼2011-12-15 20:36:14
snoopyzhao
至尊木虫 (职业作家)
- 程序强帖: 16
- 应助: 157 (高中生)
- 贵宾: 0.02
- 金币: 18844.7
- 红花: 29
- 帖子: 3803
- 在线: 1422.4小时
- 虫号: 183750
- 注册: 2006-02-13
- 专业: 污染生态化学
4楼2011-12-15 20:44:31
|
C ......... ISOTHERM CODE: C C C TOTH ISOTHERM CODE: C C THIRD ISOTHERM: NISO(IC)=3 C 3 PARAMETERS: NPISO(NISO(IC))= 3 C C PARAMETER STORAGE: C PARS(1,IC)= SATURATION LOADING QSAT (MOL/KG) OF IC C PARS(2,IC)= ADSORPION CONSTANT KADS (Pa^-1) OF IC C PARS(3,IC)= TN TOTH PARAMETER OF IC C QSAT = PARS(1,IC) PK = PARS(2,IC) TN = PARS(3,IC) DEN = (ONE + (PK*P(IC))**TN)**(ONE/TN) IF (IFLAG.EQ.0) THEN ! Calculate loading Qi at Pi Q(IC) = QSAT*PK*P(IC)/DEN ELSE IF (IFLAG .EQ. 1) THEN ! Calculate PSIi and dPSIi/dPi at Pi C... SIMPSON INTEGRATION C NTOTH = 1000 C A = P(IC) C H = A/(TWO*DFLOAT(NTOTH)) C S = ONE/TWO*FUN(A,QSAT,PK,TN) C SIMPSONLOOP: DO I=1,NTOTH-1 C S = S+TWO*FUN((TWO*DFLOAT(I)-ONE)*H,QSAT,PK,TN) C + +FUN(TWO*DFLOAT(I)*H,QSAT,PK,TN) C S = S+TWO*FUN((TWO*DFLOAT(NTOTH)-ONE)*H,QSAT,PK,TN) C ENDDO SIMPSONLOOP C H = A/DFLOAT(NTOTH) C NTOTH1 = NTOTH + 1 C C SIMPSONLOOP1: DO I = 1, NTOTH1 C YM(I) = FUN(DFLOAT(I-1)*H,QSAT,PK,TN) C ENDDO SIMPSONLOOP1 C SUM1 = ZERO C SUM2 = ZERO C SIMPSONLOOP2: DO I = 2, NTOTH, 2 C SUM1 = SUM1 + YM(I) C ENDDO SIMPSONLOOP2 C NTOTH2 = NTOTH - 1 C SIMPSONLOOP3: DO I = 3, NTOTH2, 2 C SUM2 = SUM2 + YM(I) C ENDDO SIMPSONLOOP3 C S = (YM(1)+4.D0*SUM1+2.D0*SUM2+YM(NTOTH1))/ C + DFLOAT(THREE*NTOTH) S1 = QSAT*PK*P(IC)/DEN S2 = S1/QSAT NTOTH = 0 SUM1 = 0.D0 300 CONTINUE NTOTH = NTOTH +1 C Equation from D. Do, Adsorption analysis: Equilibrium & Kinetics, C Imperial College Press: London, 1998, p.212. C TERM = S2**(1.D0+DFLOAT(NTOTH)*TN)/(1.D0+DFLOAT(NTOTH))/ C + (1.D0+DFLOAT(NTOTH)*TN) C From D. P.Valenzuela & A. L. Myers, Adsorption Equilibrium Data Handbook, C Prentice-Hall: New Jersey, 1989, p. 11. TERM = S2**(1.D0+DFLOAT(NTOTH)*TN)/(DFLOAT(NTOTH)*TN)/ + (1.D0+DFLOAT(NTOTH)*TN) SUM1 = SUM1 + TERM IF(TERM.GT.TOL) GO TO 300 C PSII(IC) = QSAT*(-(ONE/TN-ONE)*S2-S2**(ONE/TN)* C + DLOG(ONE-S2)/TN-(ONE/TN-ONE)*SUM1) PSII(IC) = QSAT*(S2-S2/TN*DLOG(ONE-S2**TN)-SUM1) DPSII(IC) = QSAT*PK/DEN ELSE IF (IFLAG .EQ. 2) THEN ! Calculate dQi/dPi at Pi DNDP(IC) = QSAT*PK/DEN/(ONE+(PK*P(IC))**TN) ELSE IF (IFLAG .EQ. 3) THEN ! Calculate Pi at Qi P(IC) = ONE/PK/((QSAT/Q(IC))**TN-ONE)**(ONE/TN) ENDIF ELSE IF (NISO(IC).EQ. 4) THEN ! Fourth ISOTHERM C |
5楼2011-12-15 21:07:19
5