应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 帮忙写一个差分程序
61/1返回列表
查看: 1179  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者夕阳西下将赠送您 20 个金币

夕阳西下

金虫 (著名写手)

[求助] 帮忙写一个差分程序

如图所示的 差分方程，用FORTRAN 程序代码如何写
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-01-31 16:14:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

harryzs

铜虫 (正式写手)
★
jjdg(金币+1): 感谢参与 2012-02-01 00:04:34
公式都有了，还写不出来，只能说你不会用fortran，其实我也不会哈哈
赞一下(1人) 回复此楼
2楼2012-01-31 20:40:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

夕阳西下

金虫 (著名写手)
★
jjdg(金币+1): 感谢参与 2012-02-01 00:04:52
引用回帖:
2楼: Originally posted by harryzs at 2012-01-31 20:40:45:
公式都有了，还写不出来，只能说你不会用fortran，其实我也不会哈哈

有公式和写代码是两回事。
赞一下(1人) 回复此楼
3楼2012-01-31 21:48:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

夕阳西下

金虫 (著名写手)
★
ben_ladeng(金币+1): 2012-02-01 20:59:35
下面的语句用的中心差分，把一阶导数、二阶导数都写成了差分格式。  
DO I=2,N-1
   HX(I)=(H1(I+1)-H1(I-1))/(2*DX)
ENDDO

HX(1)=0
HX(N)=0

DO I=2,N-1
    H2X(I)=(H1(I+1)-2*H1(I)+H1(I-1))/(DX*DX)
ENDDO

H2X(1)=0
H2X(N)=0
赞一下(1人) 回复此楼
4楼2012-02-01 11:50:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

夕阳西下

金虫 (著名写手)
program Ne

integer :: I,J, N, Ny
real::  ne(J)
real::  Hx(I)
N=1000
Ny=1000
Dx=0.1
Dt=0.1
n0=5E+2
m=2.1
p=2
Te=5
q=0.6E-19
E=2
k=0.514E-9
v=n0*k
Hx(0)=sqrt(8*Te/pm)*ne*(1-0.25)
Hx(I)=-(q/m*v)*E(ne/(1-exp(q*E*Dx/Te)+ne/(1-exp(-q*E*Dx/Te))

open (unit=10,FILE='1.txt')
DO I=2,N-1

  DO J=2,Ny-1
    HX(I)=(H1(I+1)-H1(I-1))/(2*DX)
    ne(J)=(ne(J+1)-ne(J-1))/(2*Dt)
    ne(J)=Hx(I)+R
   end do
end do


write (10,"('N=?,Ny=?')" ne(I)
stop
end
这个程序
赞一下 回复此楼
5楼2012-02-05 16:13:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

夕阳西下

金虫 (著名写手)
上面的程序错误在那里？
赞一下 回复此楼
6楼2012-02-05 16:15:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 夕阳西下 的主题更新
61/1返回列表
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 277求调剂 +22 倪建设 2026-04-06 22/1100 2026-04-12 18:18 by 学zh
[考研] 一志愿双非085400电子信息344 求调剂，对材料和化学方向也感兴趣 +11 无情的小羊 2026-04-09 12/600 2026-04-12 17:41 by 686666
[考研] 297求调剂 +17 ORCHID1 2026-04-10 18/900 2026-04-12 09:59 by zhouwenxian
[考研] 280求调剂 +13 wzzz王 2026-04-09 13/650 2026-04-12 00:31 by 勇攀高峰0126
[考研] 266求调剂，一志愿哈工程电子信息，本科获多项国奖和省奖 +8 lumine1 2026-04-06 8/400 2026-04-11 18:35 by 逆水乘风
[考研] 283求调剂 086004考英二数二 +17 那个噜子 2026-04-10 18/900 2026-04-11 16:27 by 明月此时有
[考研] 机械专硕270求调剂，接受跨专业 +12 老师看看我吧aba 2026-04-09 14/700 2026-04-11 10:21 by laoshidan
[考研] 22408调剂求助 +7 毂12 2026-04-09 9/450 2026-04-11 09:23 by 哦哦123
[考研] 342电子信息专硕求调剂 +9 你让我怎么荔枝 2026-04-10 10/500 2026-04-11 08:33 by zhq0425
[考研] 材料与化工调剂 +12 否极泰来2026 2026-04-10 13/650 2026-04-11 00:28 by wangjihu
[论文投稿] mdpi小修rvr时间四五天了 20+3 哈哈high 2026-04-08 5/250 2026-04-10 16:02 by 北京莱茵润色
[考研] 298求调剂 +13 钉叮咚冬瓜 2026-04-09 13/650 2026-04-10 15:49 by jiajinhpu
[考研] 282，电气工程专业，求调剂，不挑专业 +9 jggshjkkm 2026-04-10 9/450 2026-04-10 14:55 by 逆水乘风
[考研] 296求调剂 +6 汪！？！ 2026-04-08 6/300 2026-04-10 11:02 by mattzhming
[考研] 求调剂 +8 吃口冰激凌 2026-04-07 8/400 2026-04-09 08:03 by 5268321
[考研] 281求调剂 +10 椰子蘑菇 2026-04-06 10/500 2026-04-08 11:43 by zzucheup
[考研] 287求调剂 +6 Fnhc 2026-04-07 6/300 2026-04-08 10:05 by xingguangj
[考研] 325 调剂 +6 QQ小虾 2026-04-07 6/300 2026-04-07 15:17 by Ccclqqq
[考研] 调剂 +4 mcbbc 2026-04-06 5/250 2026-04-07 12:33 by upczlm1989
[考研] 307求调剂 +3 Youth@@ 2026-04-07 3/150 2026-04-07 09:25 by 小黑不怕难
信息提示
关闭
请填处理意见
关闭