颗粒倍数放大是出错！

黎亮

这是源程序！求高手指点！
*new
set random
def qiangneikeli
n_stiff=10e8
s_stiff=0
;；；；生成四面墙
command
         wall id 1 ks=s_stiff kn=n_stiff nodes (0,0) (5,0)
         wall id 2 ks=s_stiff kn=n_stiff nodes (5,0) (5,10)
         wall id 3 ks=s_stiff kn=n_stiff nodes (5, 10) (0,10)
         wall id 4 ks=s_stiff kn=n_stiff nodes (0,10) (0,0)
    end_command
command
          gen id=240,246 rad=0.2,0.8 x=0.2,4.8 y=02,9.8
          prop dens=1000 ks=0 kn=10e8
   end_command
   command
          gen id=247,257 rad=0.025,0.05 x=0.2,4.8 y=02,9.8
          prop dens=1000 ks=0 kn=10e8
   end_command
   command
          gen id=258,810 rad=0.0075,0.025 x=0.2,4.8 y=0.2,9.8
          prop dens=1000 ks=0 kn=10e8
   end_command
command
          gen id=811,1530 rad=0.0005,0.0075 x=0.2,4.8 y=0.2,9.8
          prop dens=1000 ks=0 kn=10e8
   end_command

    end
def get_poros
         sum = 0.0
         count=0
         bp = ball_head
     loop while bp # null
           sum = sum + pi * b_rad(bp)^2
           bp = b_next(bp)
           count=count+1
     end_loop
          pmeas = 1.0 - sum / 50;(width * height)
     end
   
qiangneikeli

get_poros
group app range x=(0.2,4.8) y=(0.2,9.8)
print count sum pmeas
plot ball wall
print pmeas ;就是这peams 值是0.8
mult = sqrt((1.0 - 0.12) / (1.0 - pmears));提示 bad command (mult)
print mult ;为什么这的mult  不可以算出等于1！
pause
     ini rad mul mult
         cycle 1000
         prop fric 0.2
           cycle 2000
         
      
print count sum pmeas
plot ball wall

ccb137731

编程需细心。pmeas pmears
一般参数需在函数内部定义才能赋值，且赋值一般都在函数内部调用。举个例子：
def aa
     b=1
end
c=b
出错，因为c未定义
再举个例
def aa
  bb=1
end
aa
bb=2
print bb
此时的结果 bb为1
就不解释了，你自己好好看看manual吧！

黎亮

谢谢了！