这就是一个非线性方程组:
k = 2; % 这里输入k f=@(x,q)q-x*(2/(k+1)*(1+(k-1)/2*x^2))^((k+1)/(2*(k-1))); q = 0:0.1:1;%这里输入q向量 x = zeros(size(q)); for i = 1:1:numel(q) x(i) = fsolve(@(x)f(x,q(i)),0.4); end x %x即是你想要的
如有帮助,还望采纳哟