% Fil temperatur.m
% Programmet beregner temperaturprofilen i kolonnen
% Input er sammensetningsvektoren fra et trinn i kolonnen, hentes fra
% utvektoren til programmet nyco7.m
% Programmet bruker funksjonen ycomp.m og sumy.m

% Programmet er laget av Per Fur og Ketil Eik
% K=[x(:,64) x(:,65) x(:,66)], x er utvektoren fra nyco7.m 
% Definer Tkglob=337*ones(1,length(K))
% Kall i matlab: I=temperatur(K)
% plot(t,I)

function Tk=temperatur(x)


global Xglob Tkglob 



tk=Tkglob;
Tkut=[];

for i=1:length(x);    

Xglob=[x(i,1); x(i,2); x(i,3) ];
Tk=fzero('sumy',tk(i));
Tkut(i)=Tk;
y(:,:,i)=ycomp(Tk);
Tkut=[Tkut];
end

Tk=Tkut;
