| 查看: 1148 | 回复: 5 | ||
| 【悬赏金币】回答本帖问题，作者夕阳西下将赠送您 20 个金币 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
夕阳西下金虫 (著名写手)
|
[求助]
帮忙写一个差分程序
|
|
如图所示的 差分方程，用FORTRAN 程序代码如何写 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有5人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有6人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有6人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有6人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有7人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有7人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 急急急，请兄弟们帮忙看一个溶解性很差的配合物该如何去培养单晶！！
已经有8人回复
- 求助：能不能帮忙写一个用matlab求一个函数的傅里叶变换的程序
已经有16人回复
- 我最近在做一个实验，错了点差错，请求师兄，师姐们帮忙！师妹非常感激！
已经有13人回复
- 两个物质，一个沸点170，一个沸点210，浓度差不多，应该怎么设置程序升温
已经有10人回复
- 写了一个fortran90的小程序，编译通不过，请大侠帮忙
已经有59人回复
- 【求助】求助哪位大侠帮我分析下电子密度差分图
已经有16人回复
- 【求助】请教虫友关于matlab差分方程的求解和画图程序
已经有4人回复
- 【求助】帮忙改写展宽程序
已经有11人回复
- 【求助】望做过差热的大虫帮忙看下
已经有9人回复
- 【求助】求高手帮忙改写数据实时动态显示程序~
已经有8人回复
- 急求帮我分析一个关于卡方的引起的误差的问题？？？
已经有2人回复
夕阳西下
金虫 (著名写手)
- 应助: 138 (高中生)
- 贵宾: 0.041
- 金币: 9247.1
- 散金: 73
- 红花: 28
- 沙发: 1
- 帖子: 1268
- 在线: 1164.4小时
- 虫号: 220750
- 注册: 2006-03-20
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 物理
3楼2012-01-31 21:48:46
2楼2012-01-31 20:40:45
夕阳西下
金虫 (著名写手)
- 应助: 138 (高中生)
- 贵宾: 0.041
- 金币: 9247.1
- 散金: 73
- 红花: 28
- 沙发: 1
- 帖子: 1268
- 在线: 1164.4小时
- 虫号: 220750
- 注册: 2006-03-20
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 物理
4楼2012-02-01 11:50:19
夕阳西下
金虫 (著名写手)
- 应助: 138 (高中生)
- 贵宾: 0.041
- 金币: 9247.1
- 散金: 73
- 红花: 28
- 沙发: 1
- 帖子: 1268
- 在线: 1164.4小时
- 虫号: 220750
- 注册: 2006-03-20
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 物理
|
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
10