pom=pom;
Remove["Global`*"];
SetDirectory[NotebookDirectory[]];
$HistoryLength=3;
mm=0.001;
Tk=273.15;
σ=5.67*(10^(-8));

Dchar=22.2mm;
o=Pi*Dchar;
ϵ=0.95;
xka=mm*{0,42,85.3,133.2,181.9,230.51,269.7,314.2};
xmax=xka[[-1]];
S=o*xmax;
ui={22.65,1.01};
P=Times@@ui;
Tcka={53,63,79,88,91,93,88,76};
Ta=19;
Tsteny=Mean[{17,17,20,16}];
xT=Thread[{xka,Tcka}];
T=Interpolation[xT,InterpolationOrder->1];
Plot[T[x],{x,0,xmax},AxesOrigin->{0,0}];
Prad=o*σ*ϵ*NIntegrate[(((T[x]+Tk))^4)-(((Tsteny+Tk))^4),{x,0,xmax}];
Pkonv=P-Prad;
Tstr=Mean[Tcka];
αstr=(Pkonv/(S*(Tstr-Ta)))

ui={22.65,1.01};
P=Times@@ui;
Tcka={54,62,77,86,92,91,86,76};
Ta=19.4;
Tsteny=Mean[{16.4,15,20.5,16}];
xT=Thread[{xka,Tcka}];
T=Interpolation[xT,InterpolationOrder->1];
Plot[T[x],{x,0,xmax},AxesOrigin->{0,0}];
Prad=o*σ*ϵ*NIntegrate[(((T[x]+Tk))^4)-(((Tsteny+Tk))^4),{x,0,xmax}];
Pkonv=P-Prad;
Tstr=Mean[Tcka];
αstr=(Pkonv/(S*(Tstr-Ta)))

ui={32.2,1.407};
P=Times@@ui;
Tcka={76,93,114,132,142,144,135,114};
Ta=19.8;
Tsteny=Mean[{24,17,16.4,16}];
xT=Thread[{xka,Tcka}];
T=Interpolation[xT,InterpolationOrder->1];
Plot[T[x],{x,0,xmax},AxesOrigin->{0,0}];
Prad=o*σ*ϵ*NIntegrate[(((T[x]+Tk))^4)-(((Tsteny+Tk))^4),{x,0,xmax}];
Pkonv=P-Prad;
Tstr=Mean[Tcka];

αstr=(Pkonv/(S*(Tstr-Ta)))

ui={11.77,0.528};
P=Times@@ui;
Tcka={32,34,38,41.5,43,43.2,41.5,38.5};
Ta=18.4;
Tsteny=Mean[{16,19.8,16,15}];
xT=Thread[{xka,Tcka}];
T=Interpolation[xT,InterpolationOrder->1];
Plot[T[x],{x,0,xmax},AxesOrigin->{0,0}];
Prad=o*σ*ϵ*NIntegrate[(((T[x]+Tk))^4)-(((Tsteny+Tk))^4),{x,0,xmax}];
Pkonv=P-Prad;
Tstr=Mean[Tcka];

αstr=(Pkonv/(S*(Tstr-Ta)))