Contents

USC 2025

clear
close all
clc

data

data=load("Eco_n1.txt");

nf=size(data,1);

ns=nf-2;

y=zeros(1,ns);
l=size(y);
k=size(y);
alpha=size(y);
beta=size(y);

c=zeros(1,ns+1);
theta=size(c);

for i=1:ns
    y(i)=sum(data(:,i));
end

LT=sum(data(ns+1,:));
KT=sum(data(ns+2,:));

Income=LT+KT;

c=data(1:ns,ns+1)';
c(ns+1)=sum(data(1:ns,ns+2));

for i=1:ns+1
    theta(i)=c(i)./Income;
end

for i=1:ns
    l(i)=data(ns+1,i);
    k(i)=data(ns+2,i);
    y(i)=sum(data(i,:));

    alpha(i)=(k(i)/l(i))/(1+(k(i)/l(i)));
    beta(i)=y(i)./(k(i)^alpha(i)*l(i)^(1-alpha(i)));
end

a=data(1:ns,1:ns)./y;
a(:,ns+1)=data(1:ns,ns+2)./c(ns+1);

if abs(sum(theta)-1) > 10^(-9)
    disp('error thetas')
    return
end

if  abs(sum(sum(beta.*k.^alpha.*l.^(1-alpha)-y))) > 10^(-4)
    disp('error alphas betas')
end

check_a=zeros(size(a));

for i=1:ns
    for j=1:ns
        check_a(i,j)=data(i,j)-a(i,j).*y(1,j);
    end
end

if  sum(sum(check_a)) > 10^(-9)
    disp('error input output coef')
end

Table

xl=(y./beta).*((1-alpha)./(alpha)).^alpha;
xk=(y./beta).*(alpha./((1-alpha))).^(1-alpha);

x=zeros(ns,ns);
for i=1:ns
    for j=1:ns
        x(i,j)=a(i,j).*y(1,j);
    end
end

c

Income2=sum(xl)+sum(xk);
c=theta.*Income2;

c sectors

c_x=c(1,1:ns)';
inv_x=a(:,ns+1).*c(1,ns+1);


data_calibra1=[x  c_x inv_x];
data_calibra2=[xl  zeros(1,ns)];
data_calibra3=[xk  zeros(1,ns)];

data_calibra=[data_calibra1;
    data_calibra2;
    data_calibra3];

error=(data_calibra-data);

if  sum(sum(error)) > 10^(-9)
    disp('eureka')
end

Published with MATLAB® R2024b