added tests/checks for Mishra/Neuman implementation #1 (fully penetration ac=ak), fixed incorrect point in documentation about input for l and d being dimensionless. adding comments to bwalya's matlab routine in an attempt to better understand it and compare to fortran code (still not same as original or FD version)