
figure(1);
xx = -5:1:5;

[X,Y] = meshgrid( xx, xx );

R3 = ( X.^2 + Y.^2 ).^(3/2);

Ex = X./R3;
Ey = Y./R3;


quiver( X, Y, Ex, Ey );

doPageFormat( [4,4] );
ylabel( 'Y (m)' );
xlabel( 'X (m)' );

print -depsc ~/Desktop/Efield.eps

%%
figure(2);

xx = -5:.01:5;

Ex = xx./abs(xx).^(3/2);

figure(2);
clf;

plot( xx, Ex );
hold on;
plot( xx, 0*Ex, ':b' );
plot( [0,0],[-10,10], ':b' );

doPageFormat( [5,3] );
ylabel( 'E_x (AU)' );
xlabel( 'x (m)' );

ylim( [-10,10] );

print -depsc ~/Desktop/Efield2.eps




