clear;

%-- To Generate plots for complex permittivity at two different
%temperatures

f = 0.1:0.1:100; %frequency in GHz (input)
t1= -1 ; % temperature in degrees C

[epsr1, epsi1]= RelDielConst_PureIce(t1,f);


t2 = -15;
[epsr2, epsi2]= RelDielConst_PureIce(t2,f);

loglog(f, epsi1, f, epsi2)
grid

