(* Piecewise linear DISCONTINOUS function.*)

W2[ϕ_]:=(Piecewise[{{43.25*(10^(-4)), ϕ≥ 1.5*(10^(-4))}, {9.33*ϕ-9.67*(10^(-4)), 0.5*(10^(-4))≤ ϕ<1.5*(10^(-4))}, {-5.005*(10^(-4)), -0.5*(10^(-4))< ϕ<0.5*(10^(-4))}, {-9.33*ϕ-9.67*(10^(-4)), -1.5*(10^(-4))< ϕ≤ -0.5*(10^(-4))}, {43.25*(10^(-4)), ϕ≤ -1.5*(10^(-4))}}])

Plot[W2[ϕ],{ϕ,-2*(10^(-4)),2*(10^(-4))},AxesLabel->{"ϕ {weber}","W(ϕ)"},PlotLabel->"DISCONTINUOUS Memristance Function"]

(* Plot charge as a function of flux *)

Q2[ϕ_]=∫W2[ϕ]\[DifferentialD]ϕ

Plot[Q2[ϕ],{ϕ,-2*(10^(-4)),2*(10^(-4))},AxesLabel->{"ϕ {weber}","Q(ϕ)"},PlotLabel->"Charge function corresonding to discontinuous memristance function"]