ymat(imin) = ymat(imin)*bet(imin)
ymat(i) = ( ymat(i)
& - amat(i)*ymat(i-1)
ymat(i) = ymat(i)
& - cmat(i)*bet(i)*ymat(i+1)
cg_Uin (i,j,1,1) = ymat(i)
ymat(i)=0.0
yMat(i) = ymat(i)+cg_Bu(i,j+2,1,1)