syms x y
eq = sqrt(y) == x;
syms t real
ysol = solve(subs(eq,x,t),y,'returnconditions',true)
The constraint tells you that t must be positive. This is because the original equation was written as sqrt(y), which presumes the positive branch of the square root function.
