disp("F*(a^(n+1))/n"); F=3.02; n=9; a=2.76; b=F*(a^(n+1))/n; Z1=1;Z2=-1;e=1.6;f=8.85;pi=3.14; printf('\n The value of b is %f*(10^-109) (Nm)^10',b); E1=(Z1*Z2*(e^2)/(4*pi*f*a)); e1=E1*10^3; E2=(b/(a^n)); E=e1+E2; printf('\n The net potentila energy of NaCl is %f*(10^-19)',E);