The talk illustrates a user-written command that extends the official kdensity to estimate density functions by the kernel method. The extensions are of two types. Firstly, the new command allows the use of an 'adaptive kernel' approach with varying, rather than fixed, bandwidths. Secondly, estimates of pointwise variability bands around the estimated density functions are computed.