clear;

%-- To Generate plots for complex permittivity at two different
%temperatures
f = 0:1:50; %frequency in GHz (input)

t1= 0; % temperature in degree C

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

t2 = 20;
[epsr2, epsi2]= RelDielConst_PureWater(t2,f);

subplot(2,1,1)
loglog(f, epsr1, f, epsi1)
grid

subplot(2,1,2)
loglog(f,epsr1, f, epsi2)
grid