clear all, clc;
h0=[0 3 10];
tidssteg=0.1;
for h=h0
    deltah=inf;%Blev onödig!!!
    H=h;
    while abs(dhdt(h))>1e-3
       h= h+tidssteg*dhdt(h);
       H=[H h];
    end
    t=0:tidssteg:tidssteg*(length(H)-1);
    plot(t,H)
    hold on
end