diff --git a/src/overlapping/trapezoid_gradients.jl b/src/overlapping/trapezoid_gradients.jl
index f5c1f91ecd222885e28aadcb1bc0ead48b57d2b5..b199e49785dac7e81047a274a9d99349c6bdd2ba 100644
--- a/src/overlapping/trapezoid_gradients.jl
+++ b/src/overlapping/trapezoid_gradients.jl
@@ -135,7 +135,7 @@ qvec(g::TrapezoidGradient, ::Nothing, ::Nothing) = δ(g) .* gradient_strength(g)
 inverse_slice_thickness(g::TrapezoidGradient) = isnothing(g.pulse) ? nothing : inverse_bandwidth(g) .* gradient_strength(g) .* 1000
 
 function variables(tg::TrapezoidGradient) 
-    list = [qvec, δ, gradient_strength, duration, rise_time, flat_time]
+    list = [slew_rate, qvec, δ, gradient_strength, duration, rise_time, flat_time]
     if !isnothing(tg.pulse)
         push!(list, inverse_slice_thickness)
     end