function[]=MakeRotatingOrbitFigures
	
	standard_gravitational_parameter = 398600.4418e9;
	sidereal_day = 86164.1;
	earth_radius = 6378100.0;
	omega = 2 * pi / sidereal_day;
	
	
	r0 = 343000 + 6378100;			% initial radius (also the minimum radius)
	u0 = 7700;						% initial velocity
	inclination = 51.6;				% inclination of the orbit
	duration = 3600*4;
		
	[time, x, y, z, period] = ExactOrbitSolution( r0, u0, inclination, duration );
	figure('Position',[0 0 1024 1024]);	
	DrawSphere
	hold on;
	plot3( x, y, z, 'Color','black','LineWidth',6 );
	hold off;
	axis equal
	view(120,12)
	zoom(1.2)
	export_fig RotatingOrbitExample1.tif -r200
	
	r0 = 6.6*6378100;			% initial radius (also the minimum radius)
	u0 = 1.0*omega*r0;			% initial velocity
	inclination = 3;			% inclination of the orbit
	duration = 86400*4;
	
	[time, x, y, z, period] = ExactOrbitSolution( r0, u0, inclination, duration );
	figure('Position',[0 0 1024 1024]);	
	DrawSphere
	hold on;
	plot3( x, y, z, 'Color','black','LineWidth',6 );
	hold off;
	axis equal
	view(98,0)
	%view(90,12)
	zoom(2.5)
	axis equal tight
	export_fig RotatingOrbitExample2.tif -r200
	
	r0 = 6.6108*6378100;
	u0 = 1.00*omega*r0;			% initial velocity
	inclination = 10;			% inclination of the orbit
	duration = 86400*1;
	
	[time, x, y, z, period] = ExactOrbitSolution( r0, u0, inclination, duration );
	figure('Position',[0 0 1024 1024]);	
	DrawSphere
	hold on;
	plot3( x, y, z, 'Color','black','LineWidth',6 );
	hold off;
	axis equal
	view(90,0)
	%view(90,12)
	zoom(3)
	export_fig RotatingOrbitExample3.tif -r200
end

function[]=DrawSphere
	[xx,yy,zz] = sphere(18);
	scale = 1.01;

	cla reset
	axis square off
	
	props.AmbientStrength = 0.8;
	props.DiffuseStrength = 1;
	props.SpecularColorReflectance = .2;
	props.SpecularExponent = 20;
	props.SpecularStrength = 1;
	props.FaceColor= 'interp';
	props.EdgeColor = [.3 .3 .3];
	props.LineWidth = 1.5;
	props.FaceLighting = 'phong';
	surface(xx,yy,zz,props);
	
	colormap(gray(512));

	light('position',[1 1 -0.1]);
end