From ce731bcbeece143df6878dfafd6b826220faa843 Mon Sep 17 00:00:00 2001
From: chrishavlin <chris.havlin@gmail.com>
Date: Fri, 12 Jul 2024 12:50:58 -0500
Subject: [PATCH] fix vectorization

---
 vbr/vbrCore/functions/Q_andrade_analytical.m | 1 +
 1 file changed, 1 insertion(+)

diff --git a/vbr/vbrCore/functions/Q_andrade_analytical.m b/vbr/vbrCore/functions/Q_andrade_analytical.m
index 4a3c011..2a0d6fd 100644
--- a/vbr/vbrCore/functions/Q_andrade_analytical.m
+++ b/vbr/vbrCore/functions/Q_andrade_analytical.m
@@ -64,6 +64,7 @@
 
   % loop over frequency
   tau_m = tau.maxwell(1:n_th);
+  Ju_in = Ju_in(1:n_th);
   for i_f = 1:n_freq
     w = 2*pi*f_vec(i_f) ;
     % get linear index of J1, J2, etc.