%% PART A
function E = coulomb(q, x1, y1, z1, x2, y2, z2)
    k = 9*10^9;
    R = sqrt((x2-x1)^2+(y2-y1)^2+(z2-z1)^2);
    E = k*q/R^2;
end 