From 5a43af5e67a06e601595e34dd1e933d801b11fb2 Mon Sep 17 00:00:00 2001 From: Austen Lamacraft Date: Mon, 21 May 2018 20:21:34 +0200 Subject: [PATCH] Working now --- src/TracyWidom.jl | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/TracyWidom.jl b/src/TracyWidom.jl index 51bbf7f..b328093 100644 --- a/src/TracyWidom.jl +++ b/src/TracyWidom.jl @@ -8,7 +8,7 @@ function F2(s::Real, N::Integer) nodes, weights = gausslegendre(N) sqrt_weights = sqrt.(weights) weights_matrix = kron(transpose(sqrt_weights),sqrt_weights) - K_matrix = [K2tilde(ξ,η,s) for ξ in nodes, η in nodes] + K_matrix = [_K2tilde(ξ,η,s) for ξ in nodes, η in nodes] det(eye(N) - weights_matrix .* K_matrix) end @@ -22,6 +22,6 @@ end _ϕ(ξ, s) = s + 10*tan(π*(ξ+1)/4) _ϕprime(ξ) = (5π/2)*(sec(π*(ξ+1)/4))^2 -_K2tilde(ξ,η,s) = sqrt(ϕprime(ξ) * ϕprime(η)) * airy_kernel(ϕ(ξ,s), ϕ(η,s)); +_K2tilde(ξ,η,s) = sqrt(_ϕprime(ξ) * _ϕprime(η)) * _airy_kernel(_ϕ(ξ,s), _ϕ(η,s)); end # module