| 查看: 464 | 回复: 0 | ||
| 【悬赏金币】回答本帖问题，作者夏沫有约将赠送您 20 个金币 | ||
夏沫有约新虫 (著名写手)
|
[求助]
程序运行出错求助
|
|
|
program main dimension in(2),x(2) data ntrap /20/ integer i real*8 in,c,d in(1)=-sqrt(2.0)/2.0 in(2)=sqrt(2.0)/2.0 c=0.0 d=-1.0 do i=1,ntrap,1 call external f() c=x(1) d=x(2) in(1)=in(1)-2.0*c/1.0*vdotn(in(1),x(1),1) ！新方向的x坐标表达式 in(2)=in(2)-2.0*d/1.0*vdotn(in(1),x(1),1) ！新方向的y坐标表达式 write(*,*) 'x(1)=',x(1),'x(2)=',x(2) write(*,*) 'in(1)=',in(1),'in(2)=',in(2) end do stop end subroutine external f() dimension x(2),y(2) data x/2*0.0/ b=2.0 n=2 m=10 eps=1.0e-5 call dnmtc(x,n,b,m,eps,y) write(*,*) do 10 i=1,n 10 write(*,100) i,x(i) write(*,*) 100 format(5x,'x(',i2,1x,')=',e13.6) end subroutine function f(x,n) dimension x(n),in(n) f1=x(1)*x(1)+x(2)*x(2)-1.0 f2=in(2)*x(1)-in(2)*c-in(1)*x(2)+in(1)*d f=sqrt(f1*f1+f2*f2) end subroutine dnmtc(x,n,b,m,eps,y) dimension x(n),y(n) double precision r real nrnd1 a=b k=1 r=1.0d0 z=f(x,n) 10 if (a.ge.eps) then l=l+1 do 20 i=1,n 20 y(i)=-a+2.0*a*nrnd1(r)+x(i) z1=f(y,n) k=k+1 if (z1.ge.z) then if (k.gt.m) then k=1 a=a/2.0 end if goto 10 else k=1 do 30 i=1,n 30 x(i)=y(i) z=z1 if (z.ge.eps) goto 10 end if end if return end subroutine real function nrnd1(r) double precision s,u,v,r s=65536.0 u=2053.0 v=13849.0 m=r/s r=r-m*s r=u*r+v m=r/s r=r-m*s nrnd1=r/s return end 我写了一个代码，主要是由一个初始点和初始方向确定的直线和一个圆相交，得到的下一个点继续作为初始点，又有一个新的方向作为新方向，继续和圆相交，一直循环。循环的次数为ntrap.我运行不知道为什么没有执行循环，输出in(1),in(2)也没有执行。 x(1)为x，x(2)为y，c为x初坐标，d为y初坐标，in(1)为直线斜率的x坐标,in(2)为直线斜率的y坐标。直线方程为：(x(1)-c)/in(1)=(x(2)-c)/in(2) 新人，还不太懂，代码可能写的有问题，没有报错，但结果不对，还请大神们指点！感谢感谢！ |
回复此楼» 猜你喜欢
博三一直没文章怎么办
已经有12人回复
-
蹲一个男朋友
已经有24人回复
-
又见奇葩学生，可咋办
已经有7人回复
-
河北省基金
已经有17人回复
-
今天才知道硕士考进大专比博士进本科工资高得多
已经有9人回复
-
博士去高校就是为了有寒暑假吗？
已经有20人回复
-
一个审稿人已返回，也没邀请审稿人，也没人审稿中，编辑要闹啥？
已经有12人回复
-
有博导缺学生吗，我可以当牛做马
已经有56人回复
-
申请基金代表性成果
已经有10人回复
-
来用亲身体会，一起说说年轻老师的辛苦
已经有19人回复
80