return new MultivariateMatrixFunction() {
public double[][] value(double[] inclination) {
try {
// Create new stellar source with input inclination angle
StellarSource source0 = getStellarSource(
star, massPlanet, period, timePeriastron, eccentricity, omega2, nodeAngle, inclination[0] );
// Run GOG for this source
runGog( source0 );
// Read epoch astrometry vector
double[] vector0 = readEpochAstrometryArray();
// Create new stellar source with input inclination angle + epsilon
StellarSource source1 = getStellarSource(
star, massPlanet, period, timePeriastron, eccentricity, omega2, nodeAngle, inclination[0] + angleIncrement );
// Run GOG for this source
runGog( source1 );
// Read epoch astrometry vector
double[] vector1 = readEpochAstrometryArray();