matice_vzdalenosti <- matrix(NA, 5, 5)

matice_vzdalenosti[1,] <- c(0.000, sqrt(26), sqrt(98), sqrt(25), sqrt(10)) 
matice_vzdalenosti[2,] <- c(sqrt(26), 0.000, sqrt(40), sqrt(13), sqrt(20))
matice_vzdalenosti[3,] <- c(sqrt(98), sqrt(40), 0.000, sqrt(25), sqrt(52))
matice_vzdalenosti[4,] <- c(sqrt(25), sqrt(13), sqrt(25), 0.000, sqrt(5))
matice_vzdalenosti[5,] <- c(sqrt(10), sqrt(20), sqrt(52), sqrt(5), 0.000)
matice_vzdalenosti

vektor_vzdalenosti_bodu <- c(sqrt(18), sqrt(8), sqrt(32), 1, 2)
vektor_vzdalenosti_bodu

c0 <- 2.5
c1 <- 7.5
a <- 10

sfericka_matice <- c0 + c1 * (3/2 * matice_vzdalenosti/a - 1/2 * (matice_vzdalenosti/a)^3)
sfericka_matice

matice_A <- cbind(rbind(sfericka_matice, 1), 1)
matice_A[6,6] <- 0
matice_A

sfericky_vektor <- c0 + c1 * (3/2 * vektor_vzdalenosti_bodu/a - 1/2 * (vektor_vzdalenosti_bodu/a)^3)
sfericky_vektor

vektor_b <- c(sfericky_vektor, 1)
vektor_b

vahy <-solve(matice_A, vektor_b)
vahy

sum(vahy[1:5])