In the present work, we developed a numerical analysis for an electroosmotic flow circulating in a rectangular microchannel considering electrolyte viscosity as a function of the induced electric field; 
which is also reflected in the slip condition imposed on the system walls, since the slip length is a function of the fluid viscosity. It should be clarified this is an entirely hydrodynamic problem, and 
for this reason there are no induced pressure gradients, because we are in the presence of a purely electroosmotic flow, where the fluid motion is due only to electrokinetic forces. Based on these comments, 
the problem is centered on high induced potentials, enabling viscoelectric effect analysis in the electroosmotic flow, which leads to significant increases in velocity and volumetric flow profiles compared to 
the case where the viscosity is a constant and there is no slip condition. Due to analytical analysis limitations, we implemented a dimensionless equation scheme defined by the continuity equation, 
the momentum equations in the x and y direction, the Poisson-Boltzmann equation, and the charge conservation equation to obtain the velocity and volumetric flow rate profiles mentioned above. This model 
is described in its variational form in order to implement the finite element technique using free software, FreeFem++.