Write PhiDivergence class and tests

PhiDivergence class can return the value of the conjugate and its
derivative, stores the maximum value of s, and can calculate the second
derivative of phi.  Included types of phi:
LRO
Burg Entropy
Kullback-Leibler
Chi^2
Modified Chi^2
Hellinger

	new file:   @PhiDivergence/PhiDivergence.m
	new file:   TestPhi.m
Update Solution class for work with PhiDivergence

Have Solution calculate s values, check for feasibility.
Update TestSolutions.
PhiDivergence now stores divergence name, update TestPhi
Correct bug dealing with limiting s value of PhiDivergence

When called with a vector, Conjugate and ConjugateDerivative now return
the correct mix of real and Inf/NaN.  TestPhi updated to correspond.
Add method to calculate phi(t) in PhiDivergence

Update test with phi(1) = 0 for all except LRO.
Fix KL value function to work with vectors
Rewrite LRLP class to work with PhiDivergence, test script TestLRLP

Changed calculations associated with second stage problems, linear cost
of lambda in LRLP.  Added method Test() to LRLP, wich is called by the
test script TestLRLP, and allows for testing the private methods of LRLP
that needed to be changed.

	new file:   @LRLP/Test.m
	new file:   TestLRLP.m
Correct initializing mu for phi* with no upper bound on s

Now initialize mu to 0, and gracefully handle case when feasibility cuts
matrix is empty.  Add code to @LRLP/Test.m to check that
SolveMasterProblem functions without error.
Update SolveLRLP to correspond with using phi-divergences
Input rho < 0 sets automatic selection of rho by phi''(1)
Update TestTwoStageGraphical to work with PhiDivergence
Change formatting at end of TestTwoStage
Bug fix in LPModel.m
Bug fix in SolveLRLP.m
Update GatherLROData.m to work with PhiDivergence class
Change TestOneStageMatch.m to use assertElementsAlmostEqual

Necessary for test to pass on MORE lab machine copper.  The tests
changed only check the output of linprog, and thus do not need to be
exactly accurate.
Correct bug for mchi2 with s < -2

Previously, incorrectly used 'max' to find values for x < -2.  Changed
to using a logical system.
Change GatherLROData.m to use alpha to generate rho
Change upper bound limit on s in Solution to a hard inequality

This is correct for most phi-divergences, except variation distance.  If
variation distance is included later, this will have to be changed.
Change PhiDivergence Value method to Contribution method

Value() calculated the value at a given number, but this would produce
NaN when the denominator was zero.  Contribution() now takes a numerator
and denominator, and handles gracefully the case when the denominator is
zero.