An N-body code containing live stellar evolution through combination of the software packages nbody6 and stars is presented. Operational details of the two codes are outlined and the changes that have been made to combine them discussed. We have computed the evolution of clusters of 10⁴ stars using the combined code and we compare the results with those obtained using nbody6 and the synthetic stellar evolution code sse. We find that, providing the physics package within stars is set up correctly to match the parameters of the models used to construct sse, the results are very similar. This provides a good indication that the new code is working well. We also demonstrate how this physics can be changed simply in the new code with convective overshooting as an example. Similar changes in sse would require considerable reworking of the model fits. We conclude by outlining proposed future development of the code to include more complete models of single stars and binary star systems.