using TyMath
using TyPlot

eta_0=[0,0.383,0.707,0.921,1]

CLD_Lamda0P25=[1.4,1.185,0.83,0.46,0]
CLD_Lamda0P5 = [1.305,1.15,0.89,0.54,0]
CLD_Lamda0P4 = zeros(5,1)
k=zeros(5,1)
b=zeros(5,1)
for m=1:4
    k[m]=(CLD_Lamda0P5[m]-CLD_Lamda0P25[m])/(0.5-0.25)
    b[m]=CLD_Lamda0P5[m]
    CLD_Lamda0P4[m]=k[m]*(0.4-0.5)+b[m]
end

p =  polyfit(eta_0,CLD_Lamda0P4,3)