{-# OPTIONS -Wall #-}

import SimpleVec ( vec )
import Electricity ( elementaryCharge )
import Lorentz ( ParticleFieldState(..), animatePFS, defaultPFS )

main :: IO ()
main = animatePFS 1e-5 30 0.05
       ( defaultPFS { mass          = 1.673e-27  -- proton in kg
                    , charge        = elementaryCharge
                    , velocity      = vec 0 2000 0
                    , electricField = \_ -> vec 0 20 0
                    , magneticField = \_ -> vec 0  0 0.01 } )