Create a palette which fades nicely eg: black->blue->cyan->white->black then plot your pixels at x, y palette index (x*x)+(y*y) it will create a load of circle patterns.

What you wrote is just the equation for the circle:

radius = (x-cx)² + (y-cy)²