#!/usr/bin/env python
from math import *

E=1E40 ## energie
r=5E10 ## polomer 
dlr=1E7  ## zmena polomeru
ro=1E-10  ## hustota
G=6.6732E-11
cs=200. ## rychlost zvuku
i=1 ## krok

   
while i < 1000:
	i +=1
	ddr=(ro*cs*cs+E/(2*pi*r*r*r)+(0.5*ro-1)*dr*dr+2*G*ro*ro*r*r/9)/(ro*r/3)*dlr
	dr=(((4*pi*r*r*(ro*cs*cs+0.5*E/(pi*r*r*r)+2*G*pi*ro*ro*r*r/9)/(1-2*pi*r*r*ro))-4*pi*ro*r*r*r*ddr/3)**0.5)*dlr
	r +=dr

   print r, dr, ddr 
