I compiled the svdcmp() code copied from Numerical Recipes in C book(On-Line) using g++ at a UNIX machine. The two bugs reported have been fixed. The fifth singular value I got was sometimes negative, and it gave me "bus error" or "fragmentation fault" when I tried to printf matrix u or even my pre-saved matrix a, yet there is no such problem when printing matrix v. Could somebody please tell me what's wrong with the code? Thanks.