计算程序
1个回答
2020-01-19 · 技术研发知识服务融合发展。
中地数媒
中地数媒（北京）科技文化有限责任公司奉行创新高效、以人为本的企业文化，坚持内容融合技术，创新驱动发展的经营方针，以高端培训、技术研发和知识服务为发展方向，旨在完成出版转型、媒体融合的重要使命
向TA提问
关注
![]()
展开全部
c name dfc3.for，无界承压含水层非稳定流抽水试验雅各布公式求T、S值
c 计算公式：s＝Q/（4πT）Ln（2.25Tt/（Srr）
dimension t0（500），s1（500），t1（500）
open（1，file＝’nqr432.txt’）
read（1，*）n，q，r
close（1）
open（1，file＝’st432.txt’）
read（1，*）（t0（i），s1（i），i＝1，n）
close（1）
y＝0
x＝0
a1＝0
b1＝0
do 10 i＝1，n
t1（i）＝alog（t0（i））
y＝y+s1（i）/n
x＝x+t1（i）/n
10 continue
do 20 i＝1，n
a1＝a1+（t1（i）-x）*（s1（i）-y）
b1＝b1+（t1（i）-x）**2
20 continue
b＝a1/b1
a＝y-b*x
t＝q/b/4/3.1415926
s＝（2.25*t*exp（-4*3.1415926*t*a/q））/r/r
40 write（*，*）’t＝’，t，’（m2/d）s＝’，s
open（1，file＝’ts432.txt’）
write（1，*）’Q（m3/d）r（m）T（m2/d）S’
write（1，60）q，r，t，s
close（1）
60 format（2f10.2，f10.2，f10.5）
open（1，file＝’st0432.txt’）
do 80 i＝1，n，3
write（1，90）t0（i），s1（i），t0（i+1），s1（i+1），t0（i+2），s1（i+2）
80 continue
close（1）
90 format（3（f10.4，f7.3，9x））
stop
end
-
官方服务
- 官方网站
已赞过
已踩过<
评论
收起
你对这个回答的评价是？
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
