We investigate the relationship between the optical variability of galaxies and their distances from the centre of the Virgo Cluster using Palomar Transient Factory data. We define the ratio between the standard deviation of the galaxy brightness and the mean value of the standard deviation as a measure of a galaxy's optical variability. A sample of 814 Virgo galaxies with 230263 observations shows a monotonically decreasing trend of optical variability with increasing clustercentric distance. The variability level inside the cluster is 3.2$\sigma$ higher than the level outside. We fit the variability with a linear function and find that the data reject a distance-independent model. We examine 217 background galaxies for comparison and find no significant trend in galaxy variability. We assess the relation with Monte Carlo simulation by rebuilding the brightness of each galaxy. The simulation shows a monotonically decreasing relation for member galaxy variability and a distance-independent relation for background galaxies. Our result is consistent with the theory that the cold gas flowing inwards the cluster centre fuels AGN activity. This work is a new implementation of the method using optical variability to investigate the relation between galaxies evolution and their environment.

Binaries play key roles in determining stellar parameters and exploring stellar evolution models. We build a catalog of 88 eclipsing binaries with spectroscopic information, taking advantage of observations from both the Large Sky Area Multi-Object fiber Spectroscopic Telescope (LAMOST) and the Palomar Transient Factory (PTF) surveys. A software pipeline is constructed to identify binary candidates by examining their light curves. The orbital periods of binaries are derived from the Lomb-Scargle method. The key distinguishing features of eclipsing binaries are recognized by a new filter \textit{Flat Test}. We classify the eclipsing binaries by applying Fourier analysis on the light curves. Among all the binary stars, 13 binaries are identified as eclipsing binaries for the first time. The catalog contains information: position, primary eclipsing magnitude and time, eclipsing depth, the number of photometry and radial velocity observations, largest radial velocity difference, binary type, the effective temperature of observable star $T_{\rm eff}$, and surface gravity of observable star log \emph{g}. The false-positive probability is calculated by using both a Monte Carlo simulation and real data from the SDSS Stripe 82 Standard Catalog. The binaries in the catalog are mostly with a period of less than one day. The period distribution shows a 0.22-day cut-off which is consistent with the low probability of an eclipsing binary rotating with such a period.

We apply the vertical Jeans equation to the kinematics of Milky Way stars in the solar neighbourhood to measure the local dark matter density. More than 90,000 G- and K-type dwarf stars are selected from the cross-matched sample of LAMOST DR5 and Gaia DR2 for our analyses. The mass models applied consist of a single exponential stellar disc, a razor thin gas disc and a constant dark matter density. We first consider the simplified vertical Jeans equation which ignores the tilt term and assumes a flat rotation curve. Under a Gaussian prior on the total stellar surface density, the local dark matter density inferred from Markov Chain Monte Carlo simulations is $0.0133_{-0.0022}^{+0.0024}\ {\rm M}_{\odot}\,{\rm pc}^{-3}$. The local dark matter densities for subsamples in an azimuthal angle range of $-10^{\circ} < \phi < 5^{\circ}$ are consistent within their 1$\sigma$ errors. However, the northern and southern subsamples show a large discrepancy due to plateaux in the northern and southern vertical velocity dispersion profiles. These plateaux may be the cause of the different estimates of the dark matter density between the north and south. Taking the tilt term into account has little effect on the parameter estimations and does not explain the north and south asymmetry. Taking half of the difference of $\sigma_{z}$ profiles as unknown systematic errors, we then obtain consistent measurements for the northern and southern subsamples. We discuss the influence of the vertical data range, the scale height of the tracer population, the vertical distribution of stars and the sample size on the uncertainty of the determination of the local dark matter density.

We study the mass distribution and kinematics of the giant elliptical galaxy M87 (NGC 4486) using discrete chemo-dynamical, axisymmetric Jeans equation modelling. Our catalogue comprises 894 globular clusters (GCs) extending to a projected radius of $\sim 430$ kpc with line-of-sight velocities and colours, and Multi Unit Spectroscopic Explorer (MUSE) integral field unit data within the central $2.4$ kpc of the main galaxy. The gravitational potential for our models is a combination of a luminous matter potential with a varying mass-to-light ratio for the main galaxy, a supermassive black hole and a dark matter (DM) potential with a cusped or cored DM halo. The best-fitting models with either a cusped or a cored DM halo show no significant differences and both are acceptable. We obtain a total mass of $(2.16 \pm 0.38) \times 10^{13} M_{\odot}$ within $\sim$ 400 kpc. By including the stellar mass-to-light ratio gradient, the DM fraction increases from $\sim$ 26 percent (with no gradient) to $\sim$ 73 percent within $1\,R_e^{\rm maj}$ (major axis of half-light isophote, 14.2 kpc), and from $\sim$ 84 percent to $\sim$ 94 percent within $5\,R_e^{\rm maj}$ (71.2 kpc). Red GCs have moderate rotation with $V_{\rm max}/\sigma \sim$ 0.4, and blue GCs have weak rotation with $V_{\rm max}/\sigma \sim$ 0.1. Red GCs have tangential velocity dispersion anisotropy, while blue GCs are consistent with being nearly isotropic. Our results suggest that red GCs are more likely to be born in-situ, while blue GCs are more likely to be accreted.

In our earlier 2019 paper, we evaluated the reliability of Schwarzschild's orbit-superposition dynamical modelling method in estimating the internal mass distribution, intrinsic stellar shapes and orbit distributions of early-type galaxies (ETGs) taken from the Illustris cosmological simulation. We now apply the same techniques to galaxies taken from the integral-field survey Mapping Nearby Galaxies with APO (MaNGA), using a sample of 149 ETGs in the mass range of $10^{9.90}\sim10^{11.80} M_{\odot}$ and made up of 105 central and 44 satellite galaxies. We find that low-mass ETGs with $\log(M_*/M_{\odot})<11.1$ have an average dark matter fraction of $\sim0.2$ within one effective radius $R_{\rm e}$, tend to be oblate-like, and are dominated by rotation about their minor axis. High-mass ETGs with $\log(M_*/M_{\odot})>11.1$ have an average dark matter fraction of $\sim0.4$ within one effective radius $R_{\rm e}$, tend to be prolate-like, and are dominated by rotation about their major axis and by centrophilic orbits. The changes of internal structures within one $R_{\rm e}$ are dominated by the total stellar mass of the individual galaxies. We find no differences of internal structures between central and satellite ETGs for the same stellar masses. However, for similar stellar mass and colour distributions, we find that ETGs more prolate-like, or with more hot orbits, tend to have higher close neighbour counts at $r_p\sim40$ kpc.

We evaluate the capabilities of Schwarzschild's orbit-superposition method by applying it to galaxies from the large scale, high resolution Illustris simulation. Nine early-type galaxies with a range of triaxiality are selected, and we create mock integral field unit data for five line-of-sight projections of each galaxy. Each of the 45 mock data sets is taken as an independent observed galaxy. Using van den Bosch's 2008 triaxial Schwarzschild implementation, we assess model estimates of various galaxy properties, covering mass profiles, intrinsic shapes, stellar orbit distributions and velocity anisotropies. Total mass within $\overline{R_{\rm e}}$ is recovered well with average deviations within $\pm15$ percent. Stellar mass is underestimated by $\sim24$ percent and dark matter overestimated by $\sim38$ percent (assuming an NFW dark matter profile and allowing for degeneracy between stellar mass and dark matter mass). Using a gNFW profile, these values improve to $\sim13$ percent for stellar mass and $\sim18$ percent for dark matter. Axis ratio estimates show a moderate bias of $\Delta (b/a)=0.07$ and $\Delta (c/a)=0.14$ ($a\ge b\ge c$). Distributions of the orbit circularities $\lambda_z$ and $\lambda_x$, representing rotation about the minor and major axes, are well reconstructed. Separating orbits into thermal categories, our models match the average fractions of these categories to within $10$ percent. Velocity anisotropy is well estimated with values matching in the inner regions but becoming slightly radially biased in the outer regions. Overall, the galaxy property estimates we obtained using Schwarzschild modelling are not implausible and are representative of the simulated galaxies we modelled.

We use a sample of ~3000 galaxies from the MaNGA MPL-7 internal data release to study the alpha abundance distribution within low-redshift galaxies. We use the Lick index ratio Mgb/<Fe> as an alpha abundance indicator to study relationships between the alpha abundance distribution and galaxy properties such as effective stellar velocity dispersion within 0.3 effective radii (sigma_*), galaxy environment, and dark matter halo formation time (z_f). We find that (1) all galaxies show a tight correlation between Mgb/<Fe> and sigma_*; (2) `old' (H_beta < 3) low-sigma_* galaxies in high local density environment and inner regions within galaxy groups are enhanced in Mgb/<Fe>, while `young' (H_beta>3) galaxies and high-mass galaxies show no or less environmental dependence; (3) `old' galaxies with high-z_f show enhanced Mgb/<Fe> over low- and medium-z_f; (4) Mgb/<Fe> gradients are close to zero and show dependence on sigma_* but no obvious dependence on the environment or z_f. Our study indicates that stellar velocity dispersion or galaxy mass is the main parameter driving the Mgb/<Fe> enhancement, although environments appear to have modest effects, particularly for low- and medium-mass galaxies.

We extend Schwarzschild's dynamical modelling method to model absorption line strength data as well as the more usual luminosity and kinematic data. Our approach draws on earlier published work by the first author with the Syer and Tremaine made-to-measure (M2M) dynamical modelling method and uses similar ideas to create a chemo-Schwarzschild method. We apply our extended Schwarzschild method to the same four early type galaxies (NGC 1248, NGC 3838, NGC 4452, NGC 4551) as the chemo-M2M work, and are able to recover successfully the 2D absorption line strength for the three lines we model (Hbeta, Fe5015, Mgb). We believe that this is the first time Schwarzschild's method has been used in this way. The techniques developed can be applied to modelling other aspects of galaxies, for example age and metallicity data coming from stellar population modelling, and are not specific to absorption line strength data.

The MaNGA project has obtained IFU data for several thousand nearby galaxies, including barred galaxies. With the two dimensional spectral and kinematic information provided by IFUs, we can measure the pattern speed of a barred galaxy, which determines the bar dynamics. We apply the non-parametric method proposed by Tremaine \& Weinberg to estimate the bar pattern speed for 53 barred galaxies, making this the largest sample studied so far in this way. Our sample is selected from the MaNGA first public data release as part of SDSS Data Release 13 according mainly to the axis ratio and position angle difference between the bar and disc, while kinematic data is from the later SDSS Data Release 14. We have used both the photometric position angle from the photometric image and the kinematic position angle from the stellar velocity map to derive the pattern speed. Combining three independent bar length measurements and the circular velocity from Jeans Anisotropic modelling (JAM), we also determine the dimensionless ratio $\cal{R}$ of the corotation radius to the bar length. We find that the galaxy's position angle is the main uncertainty in determining the bar pattern speed. The kinematic position angle leads to fewer ultrafast bars than the photometric position angle, and this could be due to the method of measuring the kinematic position angle. We study the dependence of $\cal{R}$ values on galaxy properties such as the dark matter fraction from JAM modelling and the stellar age and metallicity from stellar population synthesis (SPS). A positive correlation between the bar length and bar strength is found: the longer the bar, the stronger the bar. However, no other significant correlations are found. This may result from errors in deriving the $\cal{R}$ values or from the complex formation and slowdown processes of galactic bars.

By inverting the distributions of galaxies' apparent ellipticities and misalignment angles (measured around the projected half-light radius $R_{\rm e}$) between their photometric and kinematic axes, we study the intrinsic shape distribution of 189 slow rotator early-type galaxies with stellar masses $2\times 10^{11} M_{\odot}<M_\ast<2\times 10^{12} M_{\odot}$, extracted from a sample of about 2200 galaxies with integral-field stellar kinematics from the DR14 of the SDSS-IV MaNGA IFU survey. Thanks to the large sample of slow rotators, Graham+18 showed that there is clear structure in the misalignment angle distribution, with two peaks at both $0^{\circ}$ and $90^{\circ}$ misalignment (characteristic of oblate and prolate rotation respectively). Here we invert the observed distribution from Graham+18. The large sample allows us to go beyond the known fact that slow rotators are weakly triaxial and to place useful constraints on their intrinsic triaxiality distribution (around $1R_{\rm e}$) for the first time. The shape inversion is generally non-unique. However, we find that, for a wide set of model assumptions, the observed distribution clearly requires a dominant triaxial-oblate population. For some of our models, the data suggest a hint for a minor triaxial-prolate population, but a dominant prolate population is ruled out.

We perform full spectrum fitting stellar population analysis and Jeans Anisotropic modelling (JAM) of the stellar kinematics for about 2000 early-type galaxies (ETGs) and spiral galaxies from the MaNGA DR14 sample. Galaxies with different morphologies are found to be located on a remarkably tight mass plane which is close to the prediction of the virial theorem, extending previous results for ETGs. By examining an inclined projection (`the mass-size' plane), we find that spiral and early-type galaxies occupy different regions on the plane, and their stellar population properties (i.e. age, metallicity and stellar mass-to-light ratio) vary systematically along roughly the direction of velocity dispersion, which is a proxy for the bulge fraction. Galaxies with higher velocity dispersions have typically older ages, larger stellar mass-to-light ratios and are more metal rich, which indicates that galaxies increase their bulge fractions as their stellar populations age and become enriched chemically. The age and stellar mass-to-light ratio gradients for low-mass galaxies in our sample tend to be positive ($\rm centre<outer$), while the gradients for most massive galaxies are negative. The metallicity gradients show a clear peak around velocity dispersion $\log_{10} \sigma_{\rm e}\approx 2.0$, which corresponds to the critical mass $\sim 3\times 10^{10}M_{\odot}$ of the break in the mass-size relation. Spiral galaxies with large mass and size have the steepest gradients, while the most massive ETGs, especially above the critical mass $M_{\rm crit}\ge 2\times 10^{11} M_{\odot}$, where slow rotator ETGs start dominating, have much flatter gradients. This may be due to differences in their evolution histories, e.g. mergers.

With a sample of 48,161 K giant stars selected from the LAMOST DR 2 catalogue, we construct torus models in a large volume extending, for the first time, from the solar vicinity to a Galactocentric distance of $\sim 20$ kpc, reaching the outskirts of the Galactic disc. We show that the kinematics of the K giant stars match conventional models, e.g. as created by Binney in 2012, in the Solar vicinity. However such two-disc models fail if they are extended to the outer regions, even if an additional disc component is utilised. If we loosen constraints in the Sun's vicinity, we find that an effective thick disc model could explain the anti-centre of the MW. The LAMOST data imply that the sizes of the Galactic discs are much larger, and that the outer disc is much thicker, than previously thought, or alternatively that the outer structure is not a conventional disc at all. However, the velocity dispersion $\sigma_{0z}$ of the kinematically thick disc in the best-fitting model is about 80 km s$^{-1}$ and has a scale parameter $R_{\sigma}$ for an exponential distribution function of $\sim 19$ kpc. Such a height $\sigma_{0z}$ is strongly rejected by current measurements in the solar neighbourhood, and thus a model beyond quasi-thermal, two or three thin or thick discs is required.

We perform Jeans anisotropic modeling (JAM) on elliptical and spiral galaxies from the MaNGA DR13 sample. By comparing the stellar mass-to-light ratios estimated from stellar population synthesis (SPS) and from JAM, we find a similar systematic variation of the initial mass function (IMF) as in the earlier $\rm ATLAS^{3D}$ results. Early type galaxies (elliptical and lenticular) with lower velocity dispersions within one effective radius are consistent with a Chabrier-like IMF while galaxies with higher velocity dispersions are consistent with a more bottom heavy IMF such as the Salpeter IMF. Spiral galaxies have similar systematic IMF variations, but with slightly different slopes and larger scatters, due to the uncertainties caused by higher gas fractions and extinctions for these galaxies. Furthermore, we examine the effects of stellar mass-to-light ratio gradients on our JAM modeling, and find that the trends from our results becomes stronger after considering the gradients.

We present a study on the stellar age and metallicity distributions for 1105 galaxies using the STARLIGHT software on MaNGA integral field spectra. We derive age and metallicity gradients by fitting straight lines to the radial profiles, and explore their correlations with total stellar mass M*, NUV-r colour and environments, as identified by both the large scale structure (LSS) type and the local density. We find that the mean age and metallicity gradients are close to zero but slightly negative, which is consistent with the inside-out formation scenario. Within our sample, we find that both the age and metallicity gradients show weak or no correlation with either the LSS type or local density environment. In addition, we also study the environmental dependence of age and metallicity values at the effective radii. The age and metallicity values are highly correlated with M* and NUV-r and are also dependent on LSS type as well as local density. Low-mass galaxies tend to be younger and have lower metallicity in low-density environments while high-mass galaxies are less affected by environment.

We enhance the Syer & Tremaine made-to-measure (M2M) particle method of stellar dynamical modelling to model simultaneously both kinematic data and absorption line strength data thus creating a `chemo-M2M' modelling scheme. We apply the enhanced method to four galaxies (NGC 1248, NGC 3838, NGC 4452, NGC 4551) observed using the SAURON integral-field spectrograph as part of the ATLAS3D programme. We are able to reproduce successfully the 2D line strength data achieving mean chi^2 per bin values of ~1 with >95\% of particles having converged weights. Because M2M uses a 3D particle system, we are also able to examine the underlying 3D line strength distributions. The extent to which these distributions are plausible representations of real galaxies requires further consideration. Overall we consider the modelling exercise to be a promising first step in developing a `chemo-M2M' modelling system and in understanding some of the issues to be addressed. Whilst the made-to-measure techniques developed have been applied to absorption line strength data, they are in fact general and may be of value in modelling other aspects of galaxies.

We construct a suite of discrete chemo-dynamical models of the giant elliptical galaxy NGC 5846. These models are a powerful tool to constrain both the mass distribution and internal dynamics of multiple tracer populations. We use Jeans models to simultaneously fit stellar kinematics within the effective radius $R_{\rm e}$, planetary nebula (PN) radial velocities out to $3\, R_{\rm e}$, and globular cluster (GC) radial velocities and colours out to $6\,R_{\rm e}$. The best-fitting model is a cored DM halo which contributes $\sim 10\%$ of the total mass within $1\,R_{\rm e}$, and $67\% \pm 10\%$ within $6\,R_{\rm e}$, although a cusped DM halo is also acceptable. The red GCs exhibit mild rotation with $v_{\rm max}/\sigma_0 \sim 0.3$ in the region $R > \,R_{\rm e}$, aligned with but counter-rotating to the stars in the inner parts, while the blue GCs and PNe kinematics are consistent with no rotation. The red GCs are tangentially anisotropic, the blue GCs are mildly radially anisotropic, and the PNe vary from radially to tangentially anisotropic from the inner to the outer region. This is confirmed by general made-to-measure models. The tangential anisotropy of the red GCs in the inner regions could stem from the preferential destruction of red GCs on more radial orbits, while their outer tangential anisotropy -- similar to the PNe in this region -- has no good explanation. The mild radial anisotropy of the blue GCs is consistent with an accretion scenario.

We assess the effectiveness of the Jeans-Anisotropic-MGE (JAM) technique with a state-of-the-art cosmological hydrodynamic simulation, the Illustris project. We perform JAM modelling on 1413 simulated galaxies with stellar mass M^* > 10^{10}M_{sun}, and construct an axisymmetric dynamical model for each galaxy. Combined with a Markov Chain Monte Carlo (MCMC) simulation, we recover the projected root-mean-square velocity (V_rms) field of the stellar component, and investigate constraints on the stellar mass-to-light ratio, M^*/L, and the fraction of dark matter f_{DM} within 2.5 effective radii (R_e). We find that the enclosed total mass within 2.5 R_e is well constrained to within 10%. However, there is a degeneracy between the dark matter and stellar components with correspondingly larger individual errors. The 1 sigma scatter in the recovered M^*/L is 30-40% of the true value. The accuracy of the recovery of M^*/L depends on the triaxial shape of a galaxy. There is no significant bias for oblate galaxies, while for prolate galaxies the JAM-recovered stellar mass is on average 18% higher than the input values. We also find that higher image resolutions alleviate the dark matter and stellar mass degeneracy and yield systematically better parameter recovery.

Measurement of the local dark matter density plays an important role in both Galactic dynamics and dark matter direct detection experiments. However, the estimated values from previous works are far from agreeing with each other. In this work, we provide a well-defined observed sample with 1427 G \& K type main-sequence stars from the LAMOST spectroscopic survey, taking into account selection effects, volume completeness, and the stellar populations. We apply a vertical Jeans equation method containing a single exponential stellar disk, a razor thin gas disk, and a constant dark matter density distribution to the sample, and obtain a total surface mass density of $\rm {78.7 ^{+3.9}_{-4.7}\ M_{\odot}\ pc^{-2}}$ up to 1 kpc and a local dark matter density of $0.0159^{+0.0047}_{-0.0057}\,\rm M_{\odot}\,\rm pc^{-3}$. We find that the sampling density (i.e. number of stars per unit volume) of the spectroscopic data contributes to about two-thirds of the uncertainty in the estimated values. We discuss the effect of the tilt term in the Jeans equation and find it has little impact on our measurement. Other issues, such as a non-equilibrium component due to perturbations and contamination by the thick disk population, are also discussed.

We study the dynamics of the giant elliptical galaxy M87 from the central to the outermost regions with the made-to-measure (M2M) method. We use a new catalogue of 922 globular cluster line-of- sight velocities extending to a projected radius of 180 kpc (equivalent to 25 M87 effective radii), and SAURON integral field unit data within the central 2.4 kpc. 263 globular clusters, mainly located beyond 40 kpc, are newly observed by the Next Generation Virgo Survey (NGVS). For the M2M modelling, the gravitational potential is taken as a combination of a luminous matter potential with a constant stellar mass-to-light ratio and a dark matter potential modelled as a logarithmic potential. Our best dynamical model returns a stellar mass-to-light ratio in the I band of M/LI = 6.0(+ -0.3) M_sun/L_sun with a dark matter potential scale velocity of 591(+ -50) km/s and scale radius of 42(+ -10) kpc. We determine the total mass of M87 within 180 kpc to be (1.5 + - 0.2) 10^13 M_sun. The mass within 40 kpc is smaller than previous estimates determined using globular cluster kinematics that did not extend beyond 45 kpc. With our new globular cluster velocities at much larger radii, we see that globular clusters around 40 kpc show an anomalously large velocity dispersion which affected previous results. The mass we derived is in good agreement with that inferred from ROSAT X-ray observation out to 180 kpc. Within 30 kpc our mass is also consistent with that inferred from Chandra and XMM-Newton X-ray observations, while within 120 kpc it is about 20% smaller. The model velocity dispersion anisotropy beta parameter for the globular clusters in M87 is small, varying from -0.2 at the centre to 0.2 at 40 kpc, and gradually decreasing to zero at 120 kpc.

Using the potential from N-body simulations, we construct the Galactic bar models with the Schwarzschild method. By varying the pattern speed and the position angle of the bar, we find that the best-fit bar model has pattern speed $\Omega_{\rm p}=40\ \rm{km\ s^{-1}\ kpc^{-1}}$, and bar angle $\theta_{\rm bar}=45^{\circ}$. $N$-body simulations show that the best-fit model is stable for more than 1.5 Gyrs. Combined with the results in Wang et al. (2012), we find that the bar angle and/or the pattern speed are not well constrained by BRAVA data in our Schwarzschild models. The proper motions predicted from our model are slightly larger than those observed in four fields. In the future, more kinematic data from the ground and space-based observations will enable us to refine our model of the Milky Way bar.

We demonstrate how the Syer & Tremaine made-to-measure method of stellar dynamical modelling can be adapted to model a rotating galactic bar. We validate our made-to-measure changes using observations constructed from the existing Shen et al. (2010) N-body model of the Milky Way bar, together with kinematic observations of the Milky Way bulge and bar taken by the Bulge Radial Velocity Assay (BRAVA). Our results for a combined determination of the bar angle and bar pattern speed (~30 degrees and ~40 km/s/kpc) are consistent with those determined by the N-body model. Whilst the made-to-measure techniques we have developed are applied using a particular N-body model and observational data set, they are in fact general and could be applied to other Milky Way modelling scenarios utilising different N-body models and data sets. Additionally, we use the exercise as a vehicle for illustrating how N-body and made-to-measure methods might be combined into a more effective method.

We take a sample of 24 elliptical and lenticular galaxies previously analysed by the SAURON project using three-integral dynamical models created with Schwarzschild's method, and re-analyse them using the made-to-measure (M2M) method of dynamical modelling. We obtain good agreement between the two methods in determining the dynamical mass-to-light (M/L) ratios for the galaxies with over 80% of ratios differing by < 10% and over 95% differing by < 20%. We show that (M/L)_M2M is approximately equal to (M/L)_Sch. For the global velocity dispersion anisotropy parameter delta, we find similar values but with fewer of the made-to-measure models tangentially anisotropic by comparison with their SAURON Schwarzschild counterparts. Our investigation is the largest comparative application of the made-to-measure method to date.

We re-derive the made-to-measure method of Syer and Tremaine 1996, for modelling stellar systems and individual galaxies, and demonstrate how extensions to the made-to-measure method may be implemented and used. We illustrate the enhanced made-to-measure method by determining the mass-to-light ratio of a galaxy modelled as a Plummer sphere. From the standard galactic observables of surface brightness and line-of-sight velocity dispersion together with the h_4 Gauss-Hermite coefficient of the line-of-sight velocity distribution, we successfully recover the true mass-to-light ratio of our toy galaxy. Using kinematic data from Kleyna et al 2002, we then estimate the mass-to-light ratio of the dwarf spheroidal galaxy Draco achieving a V-band value of 539 \pm 136 M_{\odot} / L_{\odot}. We describe the main aspects of creating a made-to-measure galaxy model and show how the key modelling parameters may be determined.