kill(all)$
numer: true$

/* Position der Empfänger A, B, C */
/* (gleichseitiges Dreieck) */
xa: +1$  
ya: 0$  
xb: -1$  
yb: 0$  
xc: 0$  
yc: sqrt(3)$

/* Position des Senders S */
xs: 0.3$
ys: 0.9$

/* Exponent im Abstandsgesetz P(r) ~ 1/r^n */
n: 2$ 

/* Entfernungen S-A, S-B, S-C */
dSA: sqrt((xs - xa)^2 + (ys - ya)^2);
dSB: sqrt((xs - xb)^2 + (ys - yb)^2);
dSC: sqrt((xs - xc)^2 + (ys - yc)^2);

/* Empfangsstärken (f = Senderleistung) */
f: 14$   
PA: f/dSA^n;
PB: f/dSB^n;
PC: f/dSC^n;

s1: 1/PA^(1/n)$
s2: 1/PB^(1/n)$
s3: 1/PC^(1/n)$
alpha: xb$
beta:  yc$

/*------------------------------------*/
hx: (s1^2 - s2^2)/(4*alpha)$

delta: (s1^2-s3^2-2*alpha*hx)/(2*beta)$
gamma: (beta^2-alpha^2)/(2*beta)$

ga4: hx^2 + delta^2$
ga2: 2*hx*alpha + 2*delta*gamma - s1^2$
ga0: alpha^2 + gamma^2$

gb2: ga4/ga4$
gb1: ga2/ga4$
gb0: ga0/ga4$

p: gb1$
q: gb0$

ht1: -p/2 + sqrt(p^2/4-q)$
ht2: -p/2 - sqrt(p^2/4-q)$

hs1: sqrt(ht1)$
hs2: sqrt(ht2)$

x1: hs1^2*hx$
y1: (hs1^2*(s1^2-s3^2-2*alpha*hx)-alpha^2+beta^2)/(2*beta)$

x2: hs2^2*hx$
y2: (hs2^2*(s1^2-s3^2-2*alpha*hx)-alpha^2+beta^2)/(2*beta)$
/*------------------------------------*/

/* Ergebnis anzeigen */
x1;
y1;
x2; 
y2;