eps=ones(1,n_windings)*0.01;
T_WATERback=ones(1,n_windings)*T_water_hp;
T_WATERfront=T_WATER;
T_WIREback=ones(1,n_windings)*(180+273.15);
T_WIREfront=T_mean;
iteration=0;
% i split the condition as two parts
cond_1 = sum(abs(T_WATERfront-T_WATERback)>eps) == length(T_WATERfront)
cond_2 = sum(abs(T_WIREfront-T_WIREback)>eps) == length(T_WIREfront);
while (cond_1 || cond_2)
end