function [cp,lambda,M,rho,eta] = updtProp(T,p) global p_i; global hfTP; global dT; cp = cp_cpt(T,p); lambda = lambda_cpt(T,p); [M,rho] = molarM_cpt(T,p); eta = eta_cpt(T,p); %eta*=10; if 1 cp = 1005*ones(80,1); lambda = 0.025*ones(80,1); M = 29*ones(80,1); rho = p/8314.*M./T; eta = ones(80,1)*1e-5; endif endfunction if 0 T = (300:10:50000)'; p=100*100000*ones(rows(T),1); [cp,lambda,M,rho,eta] = updtProp(T,p); logrho = log(rho); logT = log(T); plot(logT,logrho) endif