man/doc/f1.R

In gripp: General Inverse Problem Platform

# Test function for the Direct Problem
# F(x) = At^2+Bt+C
#
A<-1
B<-3
C<-50
t_min <- 0
t_max<-100
t_passo<-5
result <- seq(from=t_min, to=t_max, by=t_passo)
t <- t_min
i <- 1
while(t<=t_max)
 {
 f <- A*t*t+B*t+C
 t <- t+t_passo
 result[i] <- f
 i <- i+1
 }
folder <- getwd()
setwd(.GlobalEnv$result_folder)
write.table(result, file = "result.dat", append = FALSE, quote = TRUE, sep = " ",eol = "\n", na = "NA", dec = ".", row.names = FALSE,col.names = FALSE, qmethod = c("escape", "double"),fileEncoding = "")
setwd(folder)