程序没问题，但是就是不能循环，请指教，那问题了

小白12306 · 发表于 2015-9-24 19:56:29

program main
implicit none
real x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x14,x15,x16,f1,f2,f3,f4,f5,f6,f7,f8,a,m,b,c,i
integer x(1:6)
print*,'请输入x1,x2,x3,x4,x5,x6,b:'
read*,b
do i=1,6
      read*,x(i)
f1=x1-x2+2.0*x1*x1+2.0*x1*x2+x2*x2
f2=x3-x4+2.0*x3*x3+2.0*x3*x4+x4*x4
f3=x5-x6+2.0*x5*x5+2.0*x5*x6+x6*x6
f4=max(f1,f2,f3)
if (f4==f1) then
      x7=x1
      x8=x2
      x11=(x3+x5)/2.0
x12=(x4+x6)/2.0
else if (f4==f2)then
      x7=x3
      x8=x4
      x11=(x1+x5)/2.0
x12=(x2+x6)/2.0
else
      x7=x5
      x8=x6
      x11=(x3+x1)/2.0
x12=(x4+x2)/2.0
end if
f5=min(f1,f2,f3)
if (f5==f1) then
      x9=x1
      x10=x2
else if (f5==f2) then
      x9=x3
      x10=x4
else
      x9=x5
      x10=x6
end if
f6=x11-x12+2.0*x11*x12+2.0*x11*x11+x12*x12
x13=2.0*x11-x7
x14=2.0*x12-x8
f7=x13-x14+2.0*x13*x13+2.0*x13*x14+x14*x14
if (f7<f6) then
      x15=2.0*x13-x11
      x16=2.0*x14-x12
else
      x15=0.5*x7+0.5*x11
      x16=0.5*x8+0.5*x12
end if
      f8=x15-x16+2.0*x15*x15+2.0*x15*x16+x16*x16
      if (f8<f5) then
      x1=x1
      x2=x2
      x3=x15
      x4=x16
      x5=x5
      x6=x6
      else
         x1=x1
         x2=x2
         x3=x13
         x4=x14
         x5=x5
         x6=x6
      end if
      if (f8>f5.and.f8<f4) then
         x1=x1
         x2=x2
         x3=x3
         x4=x4
         x5=x15
         x6=x16
         end if
      m=(f1-f6)**2/3.0+(f8-f6)**2/3.0+(f3-f6)**2/3.0
      a=sqrt(m)
      if (a<b) exit
      end do
      print*,a
      end

小白12306 · 发表于 2015-9-24 19:59:07

我这个循环结构，要循环好几次，但是我不知道那个do while循环从哪开始

csutjf · 发表于 2015-9-24 20:24:04

1) x(1)是数组元素，x1是变量。程序中x1~x16未赋值，缺省可能是0；
2) if (a<b) exit 这一句退出了循环

小白12306 · 发表于 2015-9-24 20:37:43

csutjf 发表于 2015-9-24 20:24
1) x(1)是数组元素，x1是变量。程序中x1~x16未赋值，缺省可能是0；
2) if (a

能具体点吗？、我就是想把这个程序循环一下，今天在图书馆调整了一天，都没法正常

[求助] 程序没问题，但是就是不能循环，请指教，那问题了