clc;
clear;
syms x y t1 t2 t3 t4 real
E1=sym('- (t1*(x - 1/10))/(abs(x - 1/10)^2 + abs(y - 3873/10000)^2 + 7396/15625)^(1/2) - (t2*(x - 1/10))/(abs(x - 1/10)^2 + abs(y + 3873/10000)^2 + 7396/15625)^(1/2) - (t4*(x + 23/250))/(abs(x + 23/250)^2 + abs(y - 23/100)^2 + 444889/1000000)^(1/2) - (t3*(x + 23/250))/(abs(x + 23/250)^2 + abs(y + 111/500)^2 + 444889/1000000)^(1/2)=0')
E2=sym('- (t1*(y - 3873/10000))/(abs(x - 1/10)^2 + abs(y - 3873/10000)^2 + 7396/15625)^(1/2) - (t2*(y + 3873/10000))/(abs(x - 1/10)^2 + abs(y + 3873/10000)^2 + 7396/15625)^(1/2) - (t4*(y - 23/100))/(abs(x + 23/250)^2 + abs(y - 23/100)^2 + 444889/1000000)^(1/2) - (t3*(y + 111/500))/(abs(x + 23/250)^2 + abs(y + 111/500)^2 + 444889/1000000)^(1/2)=0')
E3=sym('294/25 - (86*t2)/(125*(abs(x - 1/10)^2 + abs(y + 3873/10000)^2 + 7396/15625)^(1/2)) - (667*t4)/(1000*(abs(x + 23/250)^2 + abs(y - 23/100)^2 + 444889/1000000)^(1/2)) - (667*t3)/(1000*(abs(x + 23/250)^2 + abs(y + 111/500)^2 + 444889/1000000)^(1/2)) - (86*t1)/(125*(abs(x - 1/10)^2 + abs(y - 3873/10000)^2 + 7396/15625)^(1/2))=0')
E4=sym('- (t4*(y/1000 + 12949/250000))/(abs(x + 23/250)^2 + abs(y - 23/100)^2 + 444889/1000000)^(1/2) - (t3*(y/1000 - 12951/250000))/(abs(x + 23/250)^2 + abs(y + 111/500)^2 + 444889/1000000)^(1/2) - (t1*(y/100 + 75923/5000000))/(abs(x - 1/10)^2 + abs(y - 3873/10000)^2 + 7396/15625)^(1/2) - (t2*(y/100 - 75923/5000000))/(abs(x - 1/10)^2 + abs(y + 3873/10000)^2 + 7396/15625)^(1/2)=0')
E5=sym('(t1*(x/100 + 339/5000))/(abs(x - 1/10)^2 + abs(y - 3873/10000)^2 + 7396/15625)^(1/2) + (t2*(x/100 + 339/5000))/(abs(x - 1/10)^2 + abs(y + 3873/10000)^2 + 7396/15625)^(1/2) + (t4*(x/1000 - 34661/250000))/(abs(x + 23/250)^2 + abs(y - 23/100)^2 + 444889/1000000)^(1/2) + (t3*(x/1000 - 34661/250000))/(abs(x + 23/250)^2 + abs(y + 111/500)^2 + 444889/1000000)^(1/2)=0')
E6=sym('(t1*((277*x)/10000 - y/10 + 899/25000))/(abs(x - 1/10)^2 + abs(y - 3873/10000)^2 + 7396/15625)^(1/2) - (t2*((277*x)/10000 + y/10 + 899/25000))/(abs(x - 1/10)^2 + abs(y + 3873/10000)^2 + 7396/15625)^(1/2) + (t3*((26*y)/125 - (39*x)/500 + 39/1000))/(abs(x + 23/250)^2 + abs(y + 111/500)^2 + 444889/1000000)^(1/2) + (t4*((39*x)/500 + (26*y)/125 - 5083/125000))/(abs(x + 23/250)^2 + abs(y - 23/100)^2 + 444889/1000000)^(1/2)=0')
solve(E1,E2,E3,E4,E5,E6,'x','y','t1','t2','t3','4')