diff --git a/Project.toml b/Project.toml
index b7765fa..0849756 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,7 +1,7 @@
 name = "MicroMagnetic"
 uuid = "cef16ca0-16a8-11ef-389e-9fbcf1974e83"
 authors = ["Weiwei Wang <wangweiwei@ahu.edu.cn>", "Boyao Lv <lvboyao@mail.ustc.edu.cn>"]
-version = "0.3.6"
+version = "0.3.7"
 
 [deps]
 KernelAbstractions = "63c18a36-062a-441e-b654-da1e3ab1ce7c"
diff --git a/ext/CairoMakieExt.jl b/ext/CairoMakieExt.jl
index 45cc8d8..44e3c0f 100644
--- a/ext/CairoMakieExt.jl
+++ b/ext/CairoMakieExt.jl
@@ -248,6 +248,7 @@ function MicroMagnetic.jld2movie(jld_file; framerate=12, output=nothing, figsize
     function update_function(i)
         index = @sprintf("m/%d", i)
         m = reshape(data[index], 3, nx, ny, nz)
+        empty!(ax)
         return plot_m(m; dx=dx, dy=dy, fig=fig, ax=ax, kwargs...)
     end
 
diff --git a/src/llg/llg.jl b/src/llg/llg.jl
index d0a7d92..6e75084 100644
--- a/src/llg/llg.jl
+++ b/src/llg/llg.jl
@@ -242,6 +242,7 @@ where p_perp = p - (m.p)m
         #the above part is the standard LLG equation
 
         @inbounds aj = ut * a_J[I]
+        bj = ut * bj
         c1::T = (aj + alpha * bj) / (1 + alpha * alpha)
         c2::T = (bj - alpha * aj) / (1 + alpha * alpha)