We propose a method for obtaining the average values of physical quantities at temperatures different from the assumed value from a simulation calculation of the constant temperature molecular dynamics (MD) with a fixed temperature. The method is developed on the basis of Nosé's canonical MD. The present method is numerically tested in a simulation on a 13-atom Lennard-Jones cluster.