Numeryczne rozwiazanie rownania Friedmana

Program niestabilny !!!!
proba numerycznego rozwiazania rownania Friedmana (dominacja pylu ): %MATHMODE{(\frac{da}{dt})=\frac{8\pi }{3}\frac{1}{a}-k}% przy zalozeniach G=1 q0=1;

problem z pierwiastkiem przy k>0
zielony kolor rozwiazanie scisle 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))}%

mozliwosc ustalenia
k- krzywizna
dt - krok
y0, t0 - warunki poczatkowe

Java Applet