# 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)