function [eps1 eps2]=Cole_Cole(freq,epss,epsinf,fr,beta,sigmas);
eps0=8.8541878176e-12;
eps1=epsinf+(epss-epsinf)*(1+((freq/fr)^(1-beta))*sin(beta*pi/2))/(1+(freq/fr)^(2*(1-beta))+2*((freq/fr)^(1-beta))*sin(beta*pi/2));
eps2=((epss-epsinf)*((freq/fr)^(1-beta))*cos(beta*pi/2))/(1+(freq/fr)^(2*(1-beta))+2*((freq/fr)^(1-beta))*sin(beta*pi/2))+sigmas/(2*pi*freq*eps0);
