Tim Peters wrote: > > for (j=KK+KK-2;j>KK-LL;j-=2) > ul[KK+KK-1-j]=0.0,u[KK+KK-1-j]=u[j]-ul[j]; y(u[KK])! -- Greg Ewing, Computer Science Dept, University of Canterbury, Christchurch, New Zealand To get my email address, please visit my web page: http://www.cosc.canterbury.ac.nz/~greg