Function finiteDiffJacobian

Computes the Jacobian matrix using central difference method.

The Jacobian J has dimensions (residualCount × parameterCount) where:

• Each row corresponds to a residual component
• Each column corresponds to a parameter
• J[i][j] = ∂r_i / ∂p_j

Central difference is used for each partial derivative.