function [c,q,c0,q0,f0,t0,a] =chuanre_13
%chuanre Summary of this function goes here
%  t1 是1點的溫度
%  t2 是2點的溫度
%  t3 是3點的溫度
%  c  是t1點的比熱
%  q  是t1點的導(dǎo)熱系數(shù)
%  c0  是試件比熱
%  q0  是導(dǎo)熱系數(shù)
%  f0 是t0點的傅立葉準(zhǔn)則數(shù)
%  t0 是表面溫度
%  a  是換熱系數(shù)
[t1,t2,t3]=textread('t.txt','%f%f%f','delimiter','');
T=[t1,t3];
m=size(T);
t0=zeros(m(1)-1,1);
f0=zeros(m(1)-1,1);
a=zeros(m(1)-2,1);
q=11.77034+0.00833*t1+5.13709*10^(-6)*t1.^2-1.03044*10^(-8)*t1.^3+9.51184*10^(-12)*t1.^4;
c=474.73669+0.07218*t1+8.40877*10^(-5)*t1.^2-2.17759*10^(-7)*t1.^3+2.12546*10^(-10)*t1.^4;
for i=1:(m(1)-1)
    t0(i)=0.3144*c(i)*1.5*(t1(i+1)-t1(i))/q(i)+0.5*(3*t1(i)-t3(i));
end 
q0=11.77034+0.00833*t0+5.13709*10^(-6)*t0.^2-1.03044*10^(-8)*t0.^3+9.51184*10^(-12)*t0.^4;
c0=474.73669+0.07218*t0+8.40877*10^(-5)*t0.^2-2.17759*10^(-7)*t0.^3+2.12546*10^(-10)*t0.^4;
for i=1:(m(1)-1)
    f0(i)=2.5*10^4*q0(i)/(7.86*10^3*c0(i));
end
for i=1:(m(1)-2)
    a(i)=157.2*(0.5*(t0(i+1)-t0(i))+f0(i)*(t0(i)-t1(i)))*c(i)/(10-t0(i));
end