我试图通过矢量元素之间的循环逐元素求解方程来求解向量方程。我目前正在使用vpasolve来找到解决方案,但是当我循环遍历向量时,我正在努力保存每行的输出。如何将最后一个循环的输出保存到向量100x1中?谢谢您的帮助!我是matlab的新手。
output = [1:100]'
govt = ones(100,1)*.3
labor=zeros(100,1)
consumption=zeros(100,1)
gamma=.6
eta=.7
beta=.96
delta=.5
for k=1:length(output)
if mod(output(k),2)==0
output(k)=.8
else output(k)=1.2
end
end
%%%%%%my problem is down here%%%%%%%%
save=zeros(100,1)
for k=1:length(output)
vpasolve(((x^eta)*(1+delta)+(x^(eta+1))*delta)/(((output(k)*x-govt(k))^(-gamma))*(1+delta)+((output(k)*x-govt(k))^(-gamma+1))*delta)==output(k), x)
save(k)=int64(x)
end