Fix an error in the equations of motion. Allow w != kappa

New equations of motion and initial conditions.
There was previously a term K tau'^2 missing from
the lambda'' -equation. The term is proportional to
n^2 tau'^2, so the error is expected to be small,
and to not change the results in 1312.5199 significantly.
Preliminarily, the numerics seem to also support such
a conclusion.

The new equations of motion also implement w != kappa,
allowing a wider range of potentials.
Add a check on Veff in tauhFromQuarkMass

Add a check that Veff > 0 before attempting anything
else in tauhFromQuarkMass. This has a very minor
negative effect on runtime when the check is
passed, and a significant positive effect when not.

In addition, this is necessary in Mathematica 9,
which would attempt to evaluate bignums for very
large values of tauh, effectively choking the
entire computation.
Add an alternative method for numerical integration

Add an alternative method to NIntegralFunction,
which evaluates the input function on a specified
grid, forms an interpolation, and integrates that.
This is useful in some cases when the integrand
is itself a combination of interpolating functions
specified on the same grid.
Enforce order of sum in NIntegralFunction

Add the possibility, on by default, to enforce the
order of the sum to be from the lower limit to the
upper limit, even when their numerical order is
reversed, when using the IntegrationGrid -method
in NIntegralFunction. This allows some integrals
to remain numerically  stable even when they go
through several orders of magnitude change.
Fix a bug which broke FindZeroCrossings
Fix a dangling &, which broke lambdaendnt0
A number of fixes to ComputeParametricRawThermo
Add functions for retrieving data from thermo files.
Add a new method for FindFunctionRoots

Add an option to FindFunctionRoots, which gives
an explicit list of points for looking for sign
changes. This could especially be useful when
the function is defined using an interpolatingfunction.
Fix broken lahpointsFromData, ntPointsFromData

A typo completely broke lahpointFromData and 
ntPointsFromData, both would return tauhPointsFromData.
Fixed.
Add errors when non-numeric values are generated.
Fix old information in a comment
Fix kappa, should be omega, in qhFromThermoData.
Add the optimal potential according to 1309.2286

Add the optimal potential as DefineAIJKOptPot. Also
implemented the changes needed to support w != kappa,
in the potential definitions,  which required
minor changes in all the potential definitions and
in PotFromAnsatzAndBeta.
Add range control for tachyonic nt, and fix bug

Add options and limits to control the nt range used
for calculation. Also fix a bug related to finding
the computation boxes. This makes it possible to do 
computations for potentials where ntcrit -> infinity
as tauh, lah -> infinity.
Fix a mistake in the optimal potential.
Fix the previous fix
Fix some bugs related to options for initial conditions.

Fix a bug that prevented options from being passed
to VQCDBHInitialConditions. In addition, fix the
initial conditions with absolute epsilon, which
were broken and inconsistent with the automatically
adjusted ic's.
Fix a bug in tauhFromQuarkMass that caused spurious warnings
Fix limit finding in ComputeThermoCurve

ComputeThermoCurve would fail to find the numeric
limits, if the computation failed already when
finding the value of tauh. Fixed.