%messung

function messung(fun,N)

[A,q,x] = heat(fun,N);

tic;
  y = A\q;
zeit = toc

speicher = whos('A')

plot(x,y)
title(fun);
xlabel(N);
ylabel('Temperatur');

