Fixed problem in bcd_mult.
Unable to trace error in original function. Bcd_mult now rewritten
to bcd_mult2. Initial testing with previous fail cases now works.
Mainly editorial.
Rearranged the text in bcd.h to make it easier to find functions
Working version of bcd_ident.
Set of test functions for bcd library based on (u+v)(u-v)=u^2 - v^2
Many changes including code to clear spurious digits in big array.