diff --git a/Leontief.r b/Leontief.r
index 23c5014..ef18799 100644
--- a/Leontief.r
+++ b/Leontief.r
@@ -18,9 +18,37 @@ TYPE1_MULTIPLIERS_TOTAL <- colSums(MULTIPLIER_OUTPUT_TYPE1 )
 DEMAND_RANGE <- ACCOUNT_RANGE[-1:-6]
 TYPE2_A <- A_MAT[-DEMAND_RANGE,-DEMAND_RANGE]
 MULTIPLIER_OUTPUT_TYPE2  <-  solve(diag(LEN-8)-TYPE2_A)
+MULTIPLIER_OUTPUT_TYPE2
 HH_VALUE_ACCOUNTS <- (nrow(MULTIPLIER_OUTPUT_TYPE2  )-5):nrow(MULTIPLIER_OUTPUT_TYPE2  )
 TYPE2_MULTIPLIERS_TOTAL <- colSums(MULTIPLIER_OUTPUT_TYPE2[-HH_VALUE_ACCOUNTS ,])
 TYPE2_MULTIPLIERS_TOTAL 
 
+###########Value Added
+
+VALUE_ADD_MAT <- colSums(I_MAT[c('EmpComp','PropInc','OthPropInc','TOPI','SUB'),])/colSums(I_MAT)
+VALUE_ADD_MAT 
+colSums(VALUE_ADD_MAT[colnames(TYPE2_A)]*MULTIPLIER_OUTPUT_TYPE2 )
+colSums(VALUE_ADD_MAT[colnames(TYPE1_A)]*MULTIPLIER_OUTPUT_TYPE1 )
+colSums(MULTIPLIER_OUTPUT_TYPE1*VALUE_ADD_MAT[colnames(TYPE1_A)] )
+
+
+
+VA_TYPE2 <- colSums(solve(diag(nrow(TYPE2_A))-TYPE2_A*VALUE_ADD_MAT[colnames(TYPE2_A)]))
+VA_TYPE1 <- colSums(solve(diag(nrow(TYPE1_A))-TYPE1_A*VALUE_ADD_MAT[colnames(TYPE1_A)]))
+
+EmpComp_MAT <- I_MAT[c('EmpComp'),]/colSums(I_MAT)
+	EmpComp_TYPE2 <- colSums(solve(diag(nrow(TYPE2_A))-TYPE2_A*EmpComp_MAT[colnames(TYPE2_A)]))
+EmpComp_TYPE2 
+	EmpComp_TYPE1 <- colSums(solve(diag(nrow(TYPE1_A))-TYPE1_A*EmpComp_MAT[colnames(TYPE1_A)]))
+min(EmpComp_TYPE1 )
+
+
+	
+PropInc_MAT <- I_MAT[c('PropInc'),]/colSums(I_MAT)
+OthPropInc_MAT <- I_MAT[c('OthPropInc'),]/colSums(I_MAT)
+TOPI_MAT <- I_MAT[c('TOPI'),]/colSums(I_MAT)
+SUB_MAT <- I_MAT[c('SUB'),]/colSums(I_MAT)
+
+