Sorry the parameter will not be the same
a = 1 ; b = 5 ; c = 30 ;
P = 4; Q = 10; R = 3;
syms x y z
% function F = root3d(x)
eq1 = ( x^2/(a^2 + P)) + ( y^2/(b^2 + P)) + ( z^2/(c^2 + P)) -1 == 0; % for P > -c^2
eq2 = ( x^2/(a^2 + Q)) + ( y^2/(b^2 + Q)) + ( z^2/(c^2 + Q)) -1 == 0; % for -c^2 < Q < - b^2
eq3 = ( x^2/(a^2 + R)) + ( y^2/(b^2 + R)) + ( z^2/(c^2 + R)) -1 == 0 ; % for -b^2 < R < - a^2