%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% CALCULA LA ÓRBITA EN EL SISTEMA SOLAR
% RESOLVIENDO LA EDO A PARTIR DEL PERIHELIO
% Y AFELIO
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


% DATOS 
% (unidades en el sistema internacional)
% Mercurio(1), Venus(2), Tierra(3), Marte(4)
% Ceres(5), Júpiter(6), Halley(7) 
%
% SELECCIONA LOS QUE QUIERAS DIBUJAR
ind = [5,6,7];

% número de veces que pinta el planeta
N = 6;
% número de pasos en el método numérico
nst = 1200;


p = 1e10*[4.6001, 10.748, 14.710, 20.67, 38.262, 74.057, 8.7664];
a = 1e10*[6.9817, 10.894, 15.193, 24.92, 44.541, 81.652, 525.09];



for k=1:size(ind,2)
	peri = p(ind(k));
	afel = a(ind(k));

	[t,y] = gravi_edo2(peri,afel, nst);

	figure(1)
	plot(y(:,1).*cos(y(:,3)), y(:,1).*sin(y(:,3)),'LineWidth',2)
	hold on
	% Dibuja el planeta en N posiciones
	if N~=0
		v = [1:floor(nst/N):nst];
		plot(y(v,1).*cos(y(v,3)), y(v,1).*sin(y(v,3)), '.r', 'MarkerSize',30)
	end
	axis equal

end

% Dibuja el Sol
rectangle('Position',0.5*p(1)*[-1,-1,2,2],'Curvature',[1,1],'FaceColor','r')
axis equal;

hold off