public static void LinFit2(double[] x, double[] y, double[] sig, double[] a, out double chisq, LfitFunc funcs) { double[,] covar; if (sig == null) { sig = new double[x.Length]; for (int i = 0; i < sig.Length; i++) { sig[i] = 1E-2; } } NumericalRecipes.Lfit(x, y, sig, a, out covar, out chisq, funcs); }