| 查看: 1000 | 回复: 5 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
hym5701木虫 (小有名气)
|
[交流]
【求助】求fortran 求函数导数和复数积分的程序
|
||
|
那位能提供一下计算任意函数的导数和复数积分的fortran小程序啊，由于时间紧，没时间编程了，而且对fortran也不是很内行，希望大家帮帮忙谢谢了。提供地址或发到邮箱hym97353@163.com，小程序测试通过的话，20金币送上。^_^ [ Last edited by hym5701 on 2008-12-22 at 23:46 ] |
回复此楼» 猜你喜欢
环境工程 085701，267求调剂
已经有5人回复
-
275求调剂
已经有17人回复
-
334分 一志愿武理 材料求调剂
已经有6人回复
-
一志愿北京化工大学材料与化工（085600）296求调剂
已经有15人回复
-
085600 286分 材料求调剂
已经有9人回复
-
材料与化工272求调剂
已经有10人回复
-
南京大学化学调剂
已经有9人回复
-
085601材料工程找调剂
已经有13人回复
-
考研调剂
已经有6人回复
-
085600 材料与化工 329分求调剂
已经有13人回复
★ ★ ★ ★ ★ ★ ★
hym5701(金币+5,VIP+0):谢谢提供，不过想问下，这里如何体现复数
lei0736(金币+2,VIP+0):谢谢
hym5701(金币+5,VIP+0):谢谢提供，不过想问下，这里如何体现复数
lei0736(金币+2,VIP+0):谢谢
|
复数积分(来自Fortran常用算法程序集-徐士良) 变步长梯形求积法 SUBROUTINE FFFTS(A,B,F,EPS,T) DOUBLE PRECISION A,B,F,T,FA,FB,H,T1,S,X FA=F(A) FB=F(B) N=1 H=B-A T1=H*(FA+FB)/2.0 5 S=0.0 DO 10 K=0,N-1 X=A+(K+0.5)*H S=S+F(X) 10 CONTINUE T=(T1+H*S)/2.0 IF (ABS(T1-T).GE.EPS) THEN T1=T N=N+N H=H/2.0 GOTO 5 END IF RETURN END |
3楼2008-12-23 11:50:50
woshilsh
荣誉版主 (职业作家)
- 应助: 43 (小学生)
- 贵宾: 10.33
- 金币: 21115.1
- 散金: 737
- 红花: 22
- 帖子: 3777
- 在线: 639.7小时
- 虫号: 293038
- 注册: 2006-11-04
- 专业: 计算机应用技术
- 管辖: 仿真模拟
2楼2008-12-23 11:04:18
★ ★ ★
lei0736(金币+3,VIP+0):谢谢
lei0736(金币+3,VIP+0):谢谢
|
求非线性方程实根的对分法(来自Fortran常用算法程序集-徐士良) SUBROUTINE DDHRT(A,B,H,EPS,X,N,M,F) DIMENSION X(N) M=0 Z=A Y=F(Z) 10 IF ((Z.GT.B+H/2.0).OR.(M.EQ.N)) RETURN IF (ABS(Y).LT.EPS) THEN M=M+1 X(M)=Z Z=Z+H/2.0 Y=F(Z) GOTO 10 END IF Z1=Z+H Y1=F(Z1) IF (ABS(Y1).LT.EPS) THEN M=M+1 X(M)=Z1 Z=Z1+H/2.0 Y=F(Z) GOTO 10 END IF IF (Y*Y1.GT.0.0) THEN Y=Y1 Z=Z1 GOTO 10 END IF 20 IF (ABS(Z1-Z).LT.EPS) THEN M=M+1 X(M)=(Z1+Z)/2.0 Z=Z1+H/2.0 Y=F(Z) GOTO 10 END IF Z0=(Z1+Z)/2.0 Y0=F(Z0) IF (ABS(Y0).LT.EPS) THEN M=M+1 X(M)=Z0 Z=Z0+H/2.0 Y=F(Z) GOTO 10 END IF IF (Y*Y0.LT.0.0) THEN Z1=Z0 Y1=Y0 ELSE Z=Z0 Y=Y0 END IF GOTO 20 END |
4楼2008-12-23 11:59:47
king_zhc
木虫 (著名写手)
- 应助: 4 (幼儿园)
- 金币: 3739.2
- 红花: 5
- 帖子: 1018
- 在线: 132.1小时
- 虫号: 9
- 注册: 2003-03-01
- 性别: GG
- 专业: 工程热力学
5楼2009-02-27 23:54:32
