diff --git a/CHANGES b/CHANGES
--- a/CHANGES
+++ b/CHANGES
@@ -54,3 +54,6 @@
 
 0.1.2.11:
 		fixed indentation that prevented compilation with ghc-7.2
+
+0.1.3:
+		Joachim Fasting fixed Numeric.GSL.Fitting.Linear.multifit
diff --git a/hmatrix-gsl-stats.cabal b/hmatrix-gsl-stats.cabal
--- a/hmatrix-gsl-stats.cabal
+++ b/hmatrix-gsl-stats.cabal
@@ -1,5 +1,5 @@
 Name:               hmatrix-gsl-stats
-Version:            0.1.2.12
+Version:            0.1.3
 License:            GPL
 License-file:       LICENSE
 Copyright:          (c) A.V.H. McPhail 2010, 2011
diff --git a/lib/Numeric/GSL/Fitting/Linear.hs b/lib/Numeric/GSL/Fitting/Linear.hs
--- a/lib/Numeric/GSL/Fitting/Linear.hs
+++ b/lib/Numeric/GSL/Fitting/Linear.hs
@@ -128,13 +128,13 @@
          -> Vector Double          -- ^ observations
          -> (Vector Double,Matrix Double,Double)   -- ^ (coefficients,covariance,chi_sq)
 multifit x y = unsafePerformIO $ do
-               let ys = dim y 
-               cov <- createMatrix RowMajor ys ys 
-               p <- createVector ys 
+               let p = cols x
+               cov <- createMatrix RowMajor p p
+               c <- createVector p
                alloca$ \chi_sq -> do
-                   app4 (fitting_multifit chi_sq) mat x vec y vec p mat cov "multifit"
+                   app4 (fitting_multifit chi_sq) mat x vec y vec c mat cov "multifit"
                    chi_sq' <- peek chi_sq
-                   return (p,cov,chi_sq')
+                   return (c,cov,chi_sq')
 
 -----------------------------------------------------------------------------
 
@@ -148,13 +148,13 @@
            -> Vector Double          -- ^ observations
            -> (Vector Double,Matrix Double,Double)   -- ^ (coefficients,covariance,chi_sq)
 multifit_w x w y = unsafePerformIO $ do
-                   let ys = dim y 
-                   cov <- createMatrix RowMajor ys ys 
-                   p <- createVector ys 
+                   let p = cols x
+                   cov <- createMatrix RowMajor p p
+                   c <- createVector p
                    alloca$ \chi_sq -> do
-                       app5 (fitting_multifit_w chi_sq) mat x vec w vec y vec p mat cov "multifit"
+                       app5 (fitting_multifit_w chi_sq) mat x vec w vec y vec c mat cov "multifit"
                        chi_sq' <- peek chi_sq
-                       return (p,cov,chi_sq')
+                       return (c,cov,chi_sq')
 
 -----------------------------------------------------------------------------