(%i1)

k: (x-m)^2+(y-1)^2=16;

(%i2)

p: x-3*y+6=0;

(%i3)

solve([k,p],[x,y]);

(%i4)

solve(p,y);

(%i5)

rhs(%[1]);

(%i6)

subst(%,y,k);

(%i7)

expand(%);

(%i8)

%-16;

(%i9)

lhs(%);

(%i10)

ratcoef(%,x,2);

(%i11)

let(a,%o10);

(%i12)

ratcoef(%o9,x,1);

(%i13)

let(b,%o12);

(%i14)

ratcoef(%o9,x,0);

(%i15)

let(c,%o14);

(%i16)

letsimp(b^2-4*a*c);

(%i17)

%=0;

(%i18)

solve(%);

(%i19)

rhs(%[1]);

(%i20)

subst(%,m,k);

(%i21)

solve([%,p],[x,y]);

(%i22)

rhs(%o18[2]);

(%i23)

subst(%,m,k);

(%i24)

solve([%,p],[x,y]);

---

Created with wxMaxima.