IDL の plot procedure で、sympol に塗りつぶした円を使う方法。

A = FINDGEN(17) * (!PI*2/16.)
USERSYM, COS(A), SIN(A), /FILL
plot,indgen(5),psym=8

IDL help の filling → plotting symbols に書いてありました。