A:[xA,yA];

B:[xB,yB];

"Milieu de [A,B] :" ; I:(A+B)/2;

"Vecteur AB :" ; u:B-A;

"Fonction affine :"; f(x):=a*x+b;

"Pour xA<>xB : ";show(f(xA)=yA,f(xB)=yB);

t:linsolve([f(xA)=yA, f(xB)=yB], [a,b]);

a:t[1];

b:t[2];

f(x);

xA:1;yA:-3;xB:3;yB:1;

f(x);

A;

ev(A, nouns);

ev(B);

ev(a);

ev(b);

f(x);

ev(%, nouns);