clc; clear all; r=0.12e-9;//atomic radius of Se e0=8.854e-12; x=4*%pi*e0*(r^3);//the electronic polarisability of Se disp('F m^2',x,'the electronic polarisability of Se is:')