function [C11, XiC,Steo] = cost_T_fit(K,Sobs, KB, iKfit, Tdis, q, dof, Sn,D,Snfact) if nargin<7 Sn = zeros(size(Sobs)); end if nargin<6 dof = 6; end iKB = find(K<=KB,1,'last'); iKnoi = find(Sobs