function cp=libr_birerong(t,k)
a=[4.07;-5.123;2.297];
b=[9.92e-4;6.29e-3;-9.38e-3];
c=[-1.1988e-5;4.3855e-6;1.2567e-5];
for i=1:3
    p(i)=a(i)+b(i)*t+c(i)*t^2;
end
v=rot90(p,2);
cp=polyval(v,k);