In this article, an algorithm to identify LPV State Space models is proposed. The LPV State Space system is in the companion reachable canonical form. Both the state matrix and the output vector coefficients are linear combinations of a set of nonlinear basis functions dependent on the scheduling signal. This model structure, although simple, can describe accurately the behaviour of many nonlinear systems by an adequate choice of the scheduling signal. The identification algorithm minimises a quadratic criterion of the output error. Since this error is a linear function of the output vector parameters, a separable nonlinear least squares approach is used to minimise the criterion function by a gradient method. The derivatives required by the algorithm are the states of LPV systems that need to be simulated at every iteration. The effectiveness of the algorithm is assessed by two simulated examples.