function [a,b]=data(Tfilm)
A=[3.13,1.94,1.26,.86,.59,.42,.312,.180,.111,.072,.049,.0346,.0251,.0187];
A=10^6.*A;
k=[.0131,0.0143,.0157,.0167,.0181,.0192,.0203,.0225,.0246,.0265,.0284,.0303,.0320,.0337];


if (Tfilm < 300)
     lower=300;
     if rem(Tfilm,50)==0
       num=(Tfilm-50)/50+2;
       a=A(num);
       b=k(num);
      else
        while (Tfilm < lower)
           lower = lower-50;      
        end       
        num1=(lower-50)/50+2;
        num2=num1+1; 
        a=-((Tfilm-lower)/50*(A(num1)-A(num2)))+A(num1);
        b=-((Tfilm-lower)/50*(k(num1)-k(num2)))+k(num1);
     end   
  else   
     lower=1000;
     if rem(Tfilm,100)==0
       num=(Tfilm-100)/100+5;
       a=A(num);
       b=k(num);
     else
        while (Tfilm < lower)
           lower = lower-100;      
        end       
        num1=(lower-100)/100+5;
        num2=num1+1; 
        a=-((Tfilm-lower)/100*(A(num1)-A(num2)))+A(num1);
        b=-((Tfilm-lower)/100*(k(num1)-k(num2)))+k(num1);
     end   
end
   
   
   
