function SensiticityCompare()

load('TPMS_CHsample.mat', 'CH','cc','num_c')
i=1;j=1;
CHij = zeros(num_c,1);
for ic=1:num_c
    CHij(ic)=CH{ic}(i,j);
end
P=polyfit(cc.',CHij,9);
% plot(cc,CHij)
Sfit = polyval(P(1:9),cc);



load('TPMS_CHdc_sample.mat', 'dDHdc')
CHdcij = zeros(num_c,1);
for ic=1:num_c
    CHdcij(ic)=dDHdc{ic}(i,j);
end

plot(cc.',[CHij,Sfit.',CHdcij])
legend('CHij','fitting','Analytical')

end