for S=1:nX
solutions{S}=solve(Qo(S)*(1+R(S))*(No(S)-N)-V*(a*J+((X(S)*k*N)/(Kn+N))),(sqrt((Lstar^1.76)+5.2*(Nstar-Nstar_min))-(Lstar^0.88))/2.6==jstar^0.88,N,J);
temp = struct2cell(solutions{S});
vpa(horzcat(temp{:}))
end
This will use one column for each variable, but in the general case you will have difficulty figuring out exactly what order the variables will appear in -- and this makes no aattempt to label them on output.