add the floquet solver test

src/QuantumToolbox.jl

@@ -1,6 +1,7 @@
 module QuantumToolbox
 
 ## Standard Julia libraries
+using Revise
 using LinearAlgebra
 using SparseArrays
 
@@ -124,6 +125,7 @@ include("time_evolution/mcsolve.jl")
 include("time_evolution/ssesolve.jl")
 include("time_evolution/smesolve.jl")
 include("time_evolution/time_evolution_dynamical.jl")
+include("time_evolution/floquet_solver.jl")
 
 ## Other functionalities
 include("correlations.jl")

src/time_evolution/floquet_implementation.jl

@@ -0,0 +1,19 @@
+# Replicate the plots in Qutip Floquet tutorial
+
+# First box
+
+δ = 0.2 * 2 * π
+ϵ0 = 1.0 * 2 * π
+A = 2.5 * 2 * π
+ω = 1.0 * 2 * π
+H0 = - δ/2.0 * sigmax() - ϵ0/2.0 * sigmaz()
+H1 = A/2.0 * sigmaz()
+H = (H0, (H1, t -> cos(ω * t)))
+
+# Second box
+
+T = 2π / ω
+floquet_basis = FloquetBasis(H,T)
+f_energies = equasi(floquet_basis)
+f_energies
+