function [muf]=muf(T) 
  %fluid dynamic viscosity [Ns/m^2] at T [celsius] 
 
 Tsat=[0:10:370,374.15];
  %saturation temperature in celsius

 muliq=1e-6*[1786,1304,1002,798.3,653.9,547.8,467.3,404.8,355.4,315.6,283.1,254.8,231.0,210.9,194.1,179.8,167.7,157.4,148.5,140.7,133.9,127.9,122.4,117.5,112.9,108.7,104.8,101.1,97.5,94.1,90.7,87.2,83.5,79.5,75.4,69.4,62.1,51.8,41.4];

muf=interp1(Tsat,muliq,T);