[X,Y]=solve('x^2+y^2=10^2*(1.568-1.538)','x*2=arctan(y/x)')X =-0.60651959867130434289594767388317Y =1.6223852737335851185433543673061