how to generate random data from equations of electromagnetism?

how to call matlab to run k means on the data and

how to show which circle do the row of data belong to?

r := sqrt(x^2+y^2+z^2);
divE := diff(e*x/(4*Pi*r^3), x)+diff(e*x/(4*Pi*r^3), y)+diff(e*x/(4*Pi*r^3),z);
simplify(divE);

1.   can not  define  equation  when not known  E

how to demonstrate the first quantization and second quantization of maxwell equation with physics package