[t,y]=heun('satelite',0,10,[0.4,0,0.1,2],0.01);  
plot(y(:,1),y(:,3),'k') 
xlabel('x');ylabel('y'); 
title('Solucion particular del problema de dos cuerpos por el metodo de Heun.')