Skärningspunkt mellan två cirklar MATLAB

Hej!

Jag har ett problem med en uppgift i skolan som jag inte kan lyckas komma runt! Har suttit i flera timmar och letat utan framgång :/

vad jag vill finna är skärningspunkterna i x och y-led med hjälp av funktionen solve. Vad som ska göras är alltså att använda cirkelns ekvation för att räkna fram x och y koordinaterna. Hoppas det finns någon hjälpsam som kan detta då uppgiften ska vara inne innan dagens slut :/

x=linspace (-6,6)
centerX1=13;
centerY1=3;
radie1=6;
centerX2=13;
centerY2=2;
radie2=5;
axis equal
hold on
th = 0:pi/50:2*pi;
xunit1 = radie1 * cos(th) + centerX1;
yunit1 = radie1 * sin(th) + centerY1;
h = plot(xunit1, yunit1);
hold on
xunit2 = radie1 * cos(th) + centerX2;
yunit2 = radie1 * sin(th) + centerY2;
h = plot(xunit2, yunit2);
clearvars x
syms x

x1=solve(x^2+26*x+193.25==0) Det är här jag inte får rätt på ekvationen så det är egentligen bara denna som behövs hjälp med!

Uppskattar all hjälp även om den bara ger mig lite insikt på vägen så jag förhoppningsvis kan klara av detta!