# lunghezza curva P(t)=(x(t),y(t)) nel piano, da P(a) a P(b) x1 <- function(t) eval(D(body(x),"t")) y1 <- function(t) eval(D(body(y),"t")) F <- function(t) sqrt(x1(t)^2+y1(t)^2) x <- function(t) ... y <- function(t) ... a <- ... b <- ... integrate(F,a,b) # # lunghezza curva P(t)=(x(t),y(t),z(t)) nello spazio x1 <- function(t) eval(D(body(x),"t")) y1 <- function(t) eval(D(body(y),"t")) z1 <- function(t) eval(D(body(z),"t")) F <- function(t) sqrt(x1(t)^2+y1(t)^2+z1(t)^2) x <- function(t) ... y <- function(t) ... z <- function(t) ... a <- ... b <- ... integrate(F,a,b)