function sol = Heat_Transfer_Coef(CAS_lok1,CAS_lok2,T,P)
%HEAT_TRANSFER_COEF Summary of this function goes here
%   Detailed explanation goes here

%HX_Shell_Tube;

%Media in shell


%Vapour Pressure for Shell and for Tubes
vap_pres_CAS1=@(T) evalf('-C:\Users\Admin\Documents\MATLAB\DB_CHEMENG\_DBparts\_DB_vapour_pressure',@vap_pres,CAS_lok1,T);
vap_pres_CAS2=@(T) evalf('-C:\Users\Admin\Documents\MATLAB\DB_CHEMENG\_DBparts\_DB_vapour_pressure',@vap_pres,CAS_lok2,T);
%Compressibility Factor for Gas
compress_CAS1=@(T) evalf('-C:\Users\Admin\Documents\MATLAB\DB_CHEMENG\_DBparts\_DB_compressibility_factor',@compress_gas,CAS_lok1,T,P);
compress_CAS2=@(T) evalf('-C:\Users\Admin\Documents\MATLAB\DB_CHEMENG\_DBparts\_DB_compressibility_factor',@compress_gas,CAS_lok2,T,P);

fun=@(T) ((compress_CAS1(T)-1)/P);
fugacity=@(T) integral(fun,0,P);
fun2=@(T) fugacity(T)*vap_pres_CAS1(T)-P;
[N1,M]=size(CAS_lok1);
[N2,M]=size(CAS_lok2);
for i=0:N1;

   M(i)=     fsolve(@(T) fun2(T),T,optimoptions('fsolve','Display','off'));
   %Y=fsolve(@(Y) fun(Y),Y0,optimoptions('fsolve','Display','off'));

end






%Media in Shell
%alfa1= evalf('-C:\Users\Admin\Documents\MATLAB\DB_CHEMENG\_DBcorrs\_CO_GNIELINSKI',@Nu_Gnielinski_Gas,Re,Pr);
%Media in Tube
%alfa2= evalf('-C:\Users\Admin\Documents\MATLAB\DB_CHEMENG\_DBcorrs\_CO_',@Nu_Hobler,Re,Pr,ncols);
%Exchanger material
%HX_cond=HX_Shell_Tube.lambda;
%Exchanger thickness
%HX_thickness=HX_Shell_Tube.delta;

%sol=1/(1/alfa1+HX_thickness/HX_cond+1/alfa2);
sol=M;
%so
end