子程序使用do while循环后计算失败

jujudefuren

之前用的是for循环（能计算成功），现因需要计算更精确，用do while循环来控制计算精度，但改用do while循环后总计算失败，并且出现代码为1073741819的错误。子程序应该没有语法错误，我查过fortran语言中do while循环的用法，应该没有错。do while循环语句如下：

1.      do while(epsilon.gt.1e-3)
   
2.           V=V1-V2+p(i)*V1/K
   
3.           V3=V/n
   
4.           i=i+1
   
5.           p(i)=R*TEMP1/(V3-b)-a*arf/(V3**2+2*b*V3-b**2)
   
6.           epsilon=abs(p(i)-p(i-1))
   
7.       end do

复制代码

求大神指教，应该怎么修改？是子程序的问题吗？如果不是的话，应该怎样做才能使计算顺利进行？请指教，万分感谢！