function Temperaturverteilung

t=[288.15 216.65 216.65 228.65 270.65 270.65 252.65 180.65 180.65];
h=[0 11000 20000 32000 47000 52000 61000 79000 88000];

x=min(t):20:max(t);
y=interp1(t,h,x,'spline')

plot(x,y,'-')



end