#d0T0 =-(1A*1h01*(0T0-1T1)                    )/1K0
#d1T1 =-(1A*1h12*(1T1-1T2) - 1A*1h01*(1T1-0T0))/1K1
#d1T2 =-(1A*1h23*(1T2-1T3) + 1A*1h12*(1T2-1T1))/1K2

K = 10000
ref = c('01','02','03','04')
A = c(100, 100, 100,100)
U = c(5.0,5.0,5.0,5.0)
K = c(1000,1000,1000,1000)

df <- data.frame(ref, A, U, K)

print(df)


#y = c(rep(20, nrow(df)*5+2))
#y = rbind(y)
#vars = vector()

#vars = append(vars, "NTN")

# initial node temps
NTN=20
XTX=20
for(j in 1:5) {
    for(i in 1:nrow(df)) { 
        assign(paste0("S",i,"T",j),20)

    }
}
#print("S1T1")
#print(as.list(.GlobalEnv))
#vars = append(vars, "XTX")

#colnames(y) = vars

#print("f")
#print(y)
ls(all.names = TRUE)
#print(state[1:nrow(df)])


#XTX <- read.csv("iso.csv", header=FALSE,skip = 5)
#XTX <- approxfun(
#	x = XTX$V1,
#	y = XTX$V6, 
#	method = "linear", 
#	rule = 2)

#NTN=20
#XTX=20
Xh=5
Nh=5


#dNTN = (-(sum(df$A,Nh,y[,'NTN'])-sum(df$A,Nh,y[2:nrow(df)]))/K)

#dS1T1 = (-df$A[1]*((1/df$U[1]/6)*(y[,()]

# surface 1
dS1T1 = (-df$A[1]*((1/df$U[1])/6)*((get("S1T1")-get("S1T2")) + (Nh*(get("S1T1"))-get("NTN"))))
dS1T2 = (-df$A[1]*((2/df$U[1])/6)*((get("S1T2")-get("S1T3")) + (Nh*(get("S1T2"))-get("S1T1"))/(df$K[1]/3)))

print(dS1T1)
#for(i in 1:nrow(df)) {
#        assign(paste0("dS",i,"T2"), (-df$A[i]*((2/df$U[i])/6)*(get(paste0("S",i,"T2"))-get(paste0("S",i,"T3"))) +
#        (df$A[i]*((1/df$U[i])/6)*(get(paste0("S",i,"T2"))-get(paste0("S",i,"T1")))/(df$K[i]/3))))
#}

#for(i in 1:nrow(df)) {
#        assign(paste0("dS",i,"T3"), (-df$A[i]*((2/df$U[i])/6)*(get(paste0("S",i,"T3"))-get(paste0("S",i,"T4"))) +
#        (df$A[i]*((2/df$U[i])/6)*(get(paste0("S",i,"T3"))-get(paste0("S",i,"T2")))/(df$K[i]/3))))
#}

#for(i in 1:nrow(df)) {
#        assign(paste0("dS",i,"T4"), (-df$A[i]*((2/df$U[i])/6)*(get(paste0("S",i,"T4"))-get(paste0("S",i,"T5"))) +
#        (df$A[i]*((1/df$U[i])/6)*(get(paste0("S",i,"T4"))-get(paste0("S",i,"T3")))/(df$K[i]/3))))
#}

#for(i in 1:nrow(df)) {
#        assign(paste0("dS",i,"T5"), (-df$A[i]*(Xh*(get(paste0("S",i,"T5"))-get(paste0("XTX"))) +
#        (df$A[i]*((1/df$U[i])/6)*(get(paste0("S",i,"T5"))-get(paste0("S",i,"T4")))))))
#}




#ls(all.names = TRUE)
#print(dS1T2)