24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

brqhl_ing

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 262.9
帖子: 115
在线: 19.5小时
虫号: 854693

[交流] 【求助】这个带有分段函数偏微分方程怎么用for程序去解决

方程如下：

其中：

如何去解出这个偏微分方程望各位大侠帮帮忙啊给点力嘻嘻

回复此楼

» 猜你喜欢

急招9月入学博士，要有4级、最晚7月硕士毕业。精密电机驱控课题；学位材料已经有5人回复
又一批高校组建人工智能学院师资行吗不是骗人吗已经有7人回复
国自然面上和省基金B类撒花已经有18人回复
有没有学校收留已经有3人回复
312求调剂已经有3人回复
华师大读博已经有5人回复
急需审稿人！！！已经有3人回复
申博/考博已经有8人回复
295分求调剂已经有6人回复
085600材料与化工调剂已经有6人回复

» 本主题相关价值贴推荐，对您同样有帮助:

Max-Planck-Institute for Coal Research的电面通知，求指点~~ 已经有19人回复
fortran能求带字母矩阵的本征值吗？（矩阵元素带有字母，能得到形式解吗）已经有5人回复
求大虾指导一个关于atoi函数的问题已经有3人回复
利用概率密度函数生成随机变量过程中出现超几何分布函数问题已经有4人回复
fortran中到底该如何比较两个实数是否相等已经有8人回复
origin进行非线性拟合如何选择函数或者自定义合适的函数已经有19人回复
6-31G基组添加极化和弥散函数已经有14人回复
怎能让fortran调用c的fopen和fileio函数已经有3人回复
【交流】合流超比函数的 fortran程序（自己搞定了，不过也希望你来交流）已经有7人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2010-12-10 16:55:36

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

brqhl_ing

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 262.9
帖子: 115
在线: 19.5小时
虫号: 854693

引用回帖:

Originally posted by snail594 at 2010-12-11 13:46:24:
for 嵌套 if else语句，细节自己看了

看我这个我写的外部函数可有问题：
      external f
dimension  y(3),d(3),z(3,5000000),b(4,3),e(3),s(3),g(3)

   open(26,file='variable.dat')

h=0.001d0
   m=5000000
n=3
eps=0.00001d0


y(1)=1.0d0
y(2)=2.0d0
y(3)=3.0d0

   call adams(a,h,m,y,n,z,f,d,eps,b,e,s,g)
do 20 i=1,m
   t=(i-1)*h
   if(t.ge.1d2) write(26,'(4f15.6)')t,z(1,i),z(3,i)
20    continue

   close(26)

   end

subroutine f(t,y,n,d)
c    ----------------------------------------------------------------------
   dimension  y(n),d(n)
   real*8 mm

c    if (abs(y(1)).gt.2.9d0)then
  mm=0.366d0
else
mm=0.061d0
  endif
   d(1)=(1.15d0-mm*y(1)-sin(y(2))-y(3))/2.6d0
   d(2)=y(1)
   d(3)=(y(1)-y(3))/0.707d0
return
end
这才是  最后一部分重要的  上面以及以前都是模块化的程序