# Skalarprodukt
sprod <- function(y1, y2) {
     n1  <- length(y1)
     n2  <- length(y2)
     if (n1 != n2) { return("Fehler") }
     s <- 0
     for (i in 1:n1) {
       s <- s + y1[i]*y2[i]
     }
     s
}

# eine Rechnung
y1 <- 1:3
y2 <- 4:6
sprod(y1,y2)