syms x 

[x,params,conds]=solve(sin(x)==1,'ReturnConditions', true) 

solx =pi/2+2*pi*k 

params =k 

conds =in(k,'integer')

syms a b c y x

[x,y]=solve([a*x^2+b*y+c==0,a*x+2*y==4],[x,y])

x =

 ((a*b)/4-(-(a*(- a*b^2+32*b +16*c))/16)^(1/2))/a

 ((a*b)/4+(-(a*(- a*b^2+32*b +16*c))/16)^(1/2))/a

y =

 (-(a*(- a*b^2+32*b +16*c))/16)^(1/2)/2-(a*b)/8+2

 2-(-(a*(- a*b^2+32*b +16*c))/16)^(1/2)/2-(a*b)/8