clc clear disp('ooooooooooooooooooooooooo') disp('ooo Coulomb Force Law ooo') disp('ooooooooooooooooooooooooo') disp('Asc.Prof W.Chainarong') n = input('\n? Number of Charge, n ='); R = zeros(n,3); aR = zeros(n-1,3); Q = zeros(1, n); for i =1:n fprintf('P%1.0f(x,y,z)\n',i) R(i,1) = input(': x = '); R(i,2) = input(': y = '); R(i,3) = input(': z = '); Q(i) = input(': Q = '); end j = input('? Force on Ponit: P'); k = 1; for i =1:n if(i ~= j) aR(k,1) = R(j,1)-R(i,1); aR(k,2) = R(j,2)-R(i,2); aR(k,3) = R(j,3)-R(i,3); k = k+1; end end for i =1:n-1 aR(i,1:3) = aR(i,1:3)/(norm(aR(i,1:3)))^3; end F = zeros(n-1,3); K = 1/4/pi/8.8542e-12; k = 1; for i =1:n if(i ~= j) F(k,1:3) = Q(j)*Q(i)*K*aR(k,1:3); k = k+1; end end Ft = zeros(1, 3); for i =1:n-1 Ft = Ft + F(i,1:3); end display(F) display(Ft)