%PROGRAM TO CALCULATE FORCE BETWEEN TWO POINT CHARGES

Clear all;
clc;

Q1 = input('Enter the value of first point charge in coulombs = ');

Q2 = input('Enter the value of second point charge in coulombs = ');

epsilon_0 = input('Enter the value of free space permittivity =');

R = input('Enter the magnitude of distance between point charges =');

%Calculate Force Between Point Charges

disp('Force between two point charges = ');
F = (Q1*Q2)/(4*pi*epsilon_0*(R^2))

%If two point charges are joined together and then separated
%The new force between two charges 

q1 = (Q1+Q2)/2;

q2 = q1;

disp('The new force is = ')
Fnew = (q1*q2)/(4*pi*epsilon_0*(R^2))