function yprim = blad(t,y)
yprim = [y(2);...% Första ekvationen
        -y(1)^3-0.05*y(2)];%Andra ekvationen

end