function question_13
%2016-12-13
clear;clc
%You can change the value of k.
k=0.4;
exitflag=0;
f=@(x) gamma(2/x)^2/(gamma(1/x)*gamma(3/x))-k;
j=1;
for i=1:20
[x,fval,exitflag,output]=fsolve(f,unifrnd(-10,10));
if exitflag>0 & fval<1e-10
    X(j)=x;
    j=j+1;
end
end
disp(X);