Simulation memristor (charge-controlled) with no lineal relationship flux(q) = (1/3)q^3