Numeryczne rozwiazanie rownania Friedmana

Program niestabilny !!!!
próba numerycznego rozwiązania równania Friedmana (dominacja pylu ): %MATHMODE{(\frac{da}{dt})=\frac{8\pi }{3}\frac{1}{a}-k}% przy założeniach G=1 q0=1;

problem z pierwiastkiem przy k>0
zielony kolor rozwiązanie ścisłe dla k>0
%MATHMODE{a(\theta)=\frac{4\pi}{3k}(1-cos(\theta))}% %MATHMODE{t(\theta)=\frac{4\pi}{3k^{3/2}}(\theta-sin(\theta))}%

możliwość ustalenia
k- krzywizna
dt - krok
y0, t0 - warunki początkowe

Java Applet