f(0) = phi, f(n+1) = sqrt(f(n)^2 + sigma^2)