*** empty log message ***
to do:form c matrix and form u(i,j)
GammaN is missing
omega, lambda, LAMDA are missing...s-sigma, sigman+1 and m+1 to be added
small change