这是代码
CLEAR
cTALK=SET("TALK")
SET TALK OFF
cChoice="Y"
DO WHILE UPPER(cChoice)="Y"
  CLEAR
  INPUT"孔隙率="TO b
  INPUT"钻孔半径="TO m
  INPUT"湿润半径=" TO w
  INPUT"直径="TO d
  INPUT"液体压力="TO n
  INPUT"注水压力="TO v
  INPUT"粘性系数="TO k
  ?a=1+0.41*LOG10(1/b)
  ?c=d/4*[SQRT(2b/(1-b))+SQRT(2*PI()/SQRT((3)*(1-b)))-2]
  ?e=d*sprt(3/PI()-1/2)
  ?h=(3-logb)/LOG10(e/c)
  ?j=(1-b)/(b*PI()*h*c^2)/4*(2-h)
  ?t=(64*a^2*k*(4-h)*j*LOG10(m))/(PI()*h*c^4*(n-v))*[w^2*LOG10(m)-w^2*LOG10(w)+1/2*w^2]
  ?"迂曲度=",a
  ?"最大孔隙直径=",c
  ?"最小孔隙直径=",e
  ?"分形维数=",h
  ?"单位元=",j
  ?"注水时间=",t
  WAIT "继续运算吗?(Y/N)" TO cChoice
   ENDDO
  SET TALK &cTalk
  RETURN