function(price, lambda, l = 35, d = 3)
{
	options(warn=1)
	matrixX <- computeX(lambda, l, d)
	n <- length(price)
	regressionmatrix <- matrix(, n - l - d + 1, l + d)
	for(t in 1:(n - l - d + 1)) {
		print(t)
		regressionmatrix[t,  ] <- l1fit(matrixX, c(price[t:(t + l + d - 1)], rep(0, (l + d - 1))), int = F)[["coefficients"]]
	}
	return(regressionmatrix)
}
