clear all
prop2 = zeros(3,3,3,3);
for i = 1:3
  for j = i:3
    for k = j:3
      for l = k:3
        prop2(i,j,k,l) = get_interp_stiff(0.012165696685514,i,j,k,l,1,1);
      end
    end
  end
end
[Cijkl,C]=cijkl_to_cij(prop2,zeros(6,6),1);
C