functions to do simulation without requiring positive definite matrix

git-svn-id: svn://svn.r-forge.r-project.org/svnroot/bmhyd@26 edb9625f-4e0d-4859-8d74-9fd3b1da38cb
was too rigorous with error checking, threw out clades we could handle

git-svn-id: svn://svn.r-forge.r-project.org/svnroot/bmhyd@27 edb9625f-4e0d-4859-8d74-9fd3b1da38cb
now with extra print statements for debugging

git-svn-id: svn://svn.r-forge.r-project.org/svnroot/bmhyd@28 edb9625f-4e0d-4859-8d74-9fd3b1da38cb
fixed bounds for simulating connection

git-svn-id: svn://svn.r-forge.r-project.org/svnroot/bmhyd@29 edb9625f-4e0d-4859-8d74-9fd3b1da38cb
got rid of extraneous plot

git-svn-id: svn://svn.r-forge.r-project.org/svnroot/bmhyd@30 edb9625f-4e0d-4859-8d74-9fd3b1da38cb