| 查看: 10740 | 回复: 157 | |||
| 本帖产生 1 个 程序强帖 ，点击这里进行查看 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
maomao1210金虫 (正式写手)
|
[交流]
【交流】Fortran语言答疑专帖已有46人参与
|
||
[ Last edited by nono2009 on 2009-11-18 at 10:34 ] |
回复此楼» 本帖已获得的红花（最新10朵）
ynn112» 猜你喜欢
孩子确诊有中度注意力缺陷
已经有14人回复
-
三甲基碘化亚砜的氧化反应
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
-
论文投稿，期刊推荐
已经有4人回复
-
请问2026国家基金面上项目会启动申2停1吗
已经有5人回复
lirenshun
铁虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 389.8
- 帖子: 33
- 在线: 43.1小时
- 虫号: 539592
- 注册: 2008-04-04
- 性别: GG
- 专业: 核技术及其应用
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
小木虫(金币+0.5):给个红包，谢谢回帖交流
|
程序如下，是计算辛普森积分的一个程序 module INTEGRAL implicit none real, parameter :: t=3000 contains ! 产生数列 subroutine GenerateData(datas, width, f) real datas(  , widthreal, external :: f real r integer i,n n = size(datas,1) width = t/(n-1) r = 0 do i=1,n datas(i) = f(r) r = r+width end do end subroutine real function f(t) implicit none real,parameter::x1=1e-7 real,parameter::d=1e-10 real,parameter::pi=3.1415926 integer::t f=(x1/(60*t*((4*pi*d*t)**0.5)))*exp((-x1)/(4*d*t)) end function real function Simpson_Integral(datas, width) IMPLICIT NONE real datas( , widthreal sum integer i,n n = size(datas,1) if ( mod(n,2)==0 ) then write(*,*) "要有奇数条数据" stop end if sum = datas(1) + datas(n) ! 先算出头尾的和 do i=2,n-1 if ( mod(i,2)==0 ) then sum = sum + 4*datas(i) ! 把4*f(x)的部分累加起来 else sum = sum + 2*datas(i) ! 把2*f(x)的部分累加起来 end if end do Simpson_Integral = sum * width/3.0 ! SUM再乘上H/3 就好了 return end function end module ! SIMPSON法积分范例 program main use integral implicit none integer,parameter::n=1001 real,external::f real datas(n),width call generatedata(datas,width,f) write(*,*)simpson_integral(datas,width) stop end |
53楼2009-11-18 09:17:01
2楼2009-06-01 12:32:17
maomao1210
金虫 (正式写手)
- 程序强帖: 5
- 应助: 2 (幼儿园)
- 金币: 1431.3
- 散金: 242
- 红花: 16
- 沙发: 1
- 帖子: 991
- 在线: 441.5小时
- 虫号: 253215
- 注册: 2006-05-20
- 性别: MM
- 专业: 考古理论
★ ★ ★ ★ ★
wangen994(金币+2,VIP+0):鼓励一下，哈哈 6-4 11:55
gwdavid(金币+3,VIP+0):辛苦了!答疑帖加大奖励力度!;) 6-7 10:32
wangen994(金币+0,VIP+0):请你讲九月份十月份的参与应助的帖子整理附在http://emuch.net/bbs/viewthread.php?tid=1358729&fpage=1后面，以便发放津贴 11-9 19:00
wangen994(金币+2,VIP+0):鼓励一下，哈哈 6-4 11:55
gwdavid(金币+3,VIP+0):辛苦了!答疑帖加大奖励力度!;) 6-7 10:32
wangen994(金币+0,VIP+0):请你讲九月份十月份的参与应助的帖子整理附在http://emuch.net/bbs/viewthread.php?tid=1358729&fpage=1后面，以便发放津贴 11-9 19:00
|
哦。我给你举个例子吧，比如对称矩阵 A[3,3] A(1,1) A(1,2) A(1,3) A(2,1) A(2,2) A(2,3) A(3,1) A(3,2) A(3,3) 只要存储成一个一维数组即可： B(1)=A(1,1),B(2)=A(2,1),B(3)=A(2,2), B(4)=A(3,1),B(5)=A(3,2),B(6)=A(3,3). 还要记住这个： 行和列与存储该值的关系为：II=MAX(I,J)*(MAX(I,J)+1)/2+MIN(I,J), 那么B(II)==A(I,J). 不知道我讲的能听懂与否。如果不懂，继续发问。 |
3楼2009-06-01 12:41:49
anbb1009
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 555.6
- 散金: 100
- 帖子: 375
- 在线: 233小时
- 虫号: 644698
- 注册: 2008-11-03
- 专业: 凝聚态物性 II ：电子结构
4楼2009-06-05 14:35:16