functions.jl (400B)
1 #/usr/bin/julia 2 3 t(n) = [(2*(i-1)/(n-1) - 1) for i in 1:n] 4 5 function f(x) 6 s = 0 7 for (k, xk) in enumerate(x) 8 s = s .+ (xk.^2 / (8^(k-1))) 9 end 10 return (1 .+ s).^(-1) 11 end 12 13 function g(x) 14 s1 = 0 15 s2 = 0 16 for (k, xk) in enumerate(x) 17 s1 = s1 .+ (xk.^2 / (8^(k-1))) 18 s2 = s2 .+ (4 * pi * xk)/(4^(k-1)) 19 end 20 return sqrt.(s1) .* (1 .+ 1/2 * cos.(s2)) 21 end