As researchers focus on the field of quantum machine learning, there are more and more quantum algorithms for data fitting. However, the majority of quantum algorithms for data fitting are only applicable to sparse matrices, which bring a large number of limitations to their application. In this paper, to solve this problem, we propose a quantum algorithm for data fitting for general matrices, not only sparse and dense matrices. The quantum algorithm is mainly composed of three subroutines which are generating target state, quantum state tomography and estimating fit quality. The time complexity of the quantum algorithm is O κ 4 ‖ A ‖ F polylog (M N) (M ′ 2 log (M ′) 2 / ɛ 3 + 1 / (ɛ ϵ 2)) , where A ∈ R M × N , κ is the conditional number of matrix A † A ∈ R N × N , | | A | | F is the Frobenius norm of A , M ′ is the maxima number of fit functions in allowing fit, ɛ is the fitting error, and ϵ is the error of fit quality estimation. Moreover, our quantum algorithm is independent of the sparsity of the matrix, so our algorithm has wide application in dense matrices and poses a new state of the art for solving dense data fitting. • Data fitting is a mathematical method to construct the least squares fitting problem. • The proposed quantum algorithm is independent of the sparsity of the matrix. • The algorithm achieves polynomial-scale speedup compared to other algorithms. [ABSTRACT FROM AUTHOR]