x <- c(10.4, 5.6, 3.1, 6.4, 21.7)

y <- (2*x-5)^2

z <- c(x, 0, x)

sum(y)
mean(y)
sd(y)

seq(-5, 5, by=.2) -> s1

seq(length=51,from=-5, by=.2) -> s2

s3 <- 1:10

h <- 1:16

z2 <- array(h, dim=c(2,4,2))
z2

19 -> z2[1,2,1]

z3 <- array(1:4, dim=c(2,2))