function [A1, A2, A3] = ch16ex9;

% Exercise 16.9 (Fall 2011)
% experiment design

D = 4;

% source positions
x1=[-1.5; 0];  x2=[-0.5; 0]; x3 = [0.5; 0]; x4 = [1.5; 0];

% sensors positions for 1st configuration
y1=[-1.5; D];  y2=[-0.5; D]; y3 = [0.5; D]; y4 = [1.5; D];
A1 = [1/norm(x1-y1)^2  1/norm(x2-y1)^2 1/norm(x3-y1)^2 1/norm(x4-y1)^2;
    1/norm(x1-y2)^2  1/norm(x2-y2)^2 1/norm(x3-y2)^2 1/norm(x4-y2)^2;
    1/norm(x1-y3)^2  1/norm(x2-y3)^2 1/norm(x3-y3)^2 1/norm(x4-y3)^2;
    1/norm(x1-y4)^2  1/norm(x2-y4)^2 1/norm(x3-y4)^2 1/norm(x4-y4)^2];

% sensor positions for 2nd configuration
y1=[-1.5; D];  y2=[1.5; D]; y3 = [-1.5; -D]; y4 = [1.5; -D];
A2 = [1/norm(x1-y1)^2  1/norm(x2-y1)^2 1/norm(x3-y1)^2 1/norm(x4-y1)^2;
    1/norm(x1-y2)^2  1/norm(x2-y2)^2 1/norm(x3-y2)^2 1/norm(x4-y2)^2;
    1/norm(x1-y3)^2  1/norm(x2-y3)^2 1/norm(x3-y3)^2 1/norm(x4-y3)^2;
    1/norm(x1-y4)^2  1/norm(x2-y4)^2 1/norm(x3-y4)^2 1/norm(x4-y4)^2];

% sensor positions for 3rd configuration
y1=[-1.5; D];  y2=[1.5; D]; y3 = [-1.5-D; 0]; y4 = [1.5+D; 0];
A3 = [1/norm(x1-y1)^2  1/norm(x2-y1)^2 1/norm(x3-y1)^2 1/norm(x4-y1)^2;
    1/norm(x1-y2)^2  1/norm(x2-y2)^2 1/norm(x3-y2)^2 1/norm(x4-y2)^2;
    1/norm(x1-y3)^2  1/norm(x2-y3)^2 1/norm(x3-y3)^2 1/norm(x4-y3)^2;
    1/norm(x1-y4)^2  1/norm(x2-y4)^2 1/norm(x3-y4)^2 1/norm(x4-y4)^2];