clear
close all
clc
Eco_n1=load("Eco_n1.txt");
[ns, y, a, beta, alpha, theta, LT, KT ] = Calibrar(Eco_n1);
w=1;
r=1;
p=ones(1,ns);
ytot=y;
for j=1:ns
pp(j)=(1/beta(j))*((r/alpha(j))^alpha(j)*(w/(1-alpha(j))^(1-alpha(j))));
for i=1:ns
pp(j)=pp(j)+a(i,j)*p(j);
end
end
xl=(ytot./beta).*((1-alpha).*r./(alpha.*w)).^alpha;
pinv=0;
for i=1:2
pinv=pinv+a(i,1+ns)*p(i);
end
Ycon=w*LT+r*KT;
for i=1:2
c(i)=theta(i)*Ycon/p(i);
end
inv=theta(ns+1)*Ycon/pinv;
for j=1:ns
for i=1:ns
x(i,j)=a(i,j)*ytot(j);
end
end
xinv=a(:,ns+1).*inv;
ff=zeros(1,2*ns+1);
ff(1)=LT-sum(xl);
for i=1:2
ff(1+i)=ytot(i)-c(i)-xinv(i)-sum(x(i,:),2);
ff(1+ns+i)=pp(i)-p(i);
end