disp u
coor r,o
shap %1 %2
gaus %3
mate ek q 0.158;0.0;
$i
$c6 r0=0.0
$c6 do 1001 n=1,nnode
$c6 r0=r0+coorr(1,n)
$c0 1001 continue
$c6 r0=r0*1.0e-3/nnode
$c6 do 1002 n=1,nnode
$c6 if (coorr(1,n).lt.r0) coorr(1,n)=r0
$c0 1002 continue
stif
dist=[u/r;u/r]*r*ek+[u/o;u/o]/r*ek
load=[u]*r*q
end