creates a single section called soma with total length of just L=1 micron, divided into 10,000 pieces each 0.1 nm long. Perhaps you meant something more like L=1001 nseg=1001?

I'm guessing your mechanisms enae or eke are accumulation mechanisms for sodium and potassium that make their concentrations change in response to the currents? If so (and I realize this may be the point of your investigation), keep in mind that without a pump or other stabilization mechanism, after a sufficiently long period of time of continuous firing, the intracellular potassium concentration will increase and the sodium concentration will decrease to the point where the cell is no longer able to fire.

Instead of printing out the membrane potentials at each time step, it's probably better to use a Vector to record the state variables (potentials, concentrations, etc) that you're interested in so that you can analyze and plot them later.

A few other hints:

Start with simple models (e.g. can you record an action potential for a single compartment cell with Hodgkin-Huxley channels?) before proceeding to more complicated ones, especially if your just learning a tool.

You generally want to use an odd value for nseg because that allows you to easily refine your discretization (by increasing nseg by a factor of 3) and preserve all the previous node locations so you can check for convergence.

Thank you for the comments. I meant can anyone check whether the result make sense or not.
Yes, the mechanism enae and eke to calculate reversal sodium and potassium, and mechanism nadifl and kdifl to calculate sodium and potassium longitudinal diffusion ( I do not know the exactly the basic equation for the mechanism nadifl and kdifl and that is what I would like to ask about it also )
If anyone would like me to send this code for him to check I will do it, and this is the purpose of my post because I am new user for NEURON.
In regard to the dimension I need to create a huge change of concentration so I have chosen the small dimension

The equations for nadifl and kdifl are in the mod files that you're using; they're not part of NEURON itself. Without seeing them, they are likely to just be LONGITUDINAL_DIFFUSION based mechanisms, possibly with radial diffusion as well.

If by "small dimension" you're talking about the 0.1 microns for the diameter of a thin axon, okay. If you're talking about the 1/10000 microns = 0.1 nm length for the segments, that's way too small. You're discretizing into pieces that are thinner than a hydrogen atom.

The COMPARTMENT line indicates the volume per unit length (this will be needed later). The second line indicates that nai exhibits one dimensional diffusion at a rate D (in that file 0.6 um2/ms by default via Fick's law). The third line scales the current per unit area by the surface area (per unit length) and turns it into a change of mass for nai. Together with the COMPARTMENT specification of the volume per unit length, this finally turns into a change in concentration.

Note that with a total length of only 1 micron, diffusion will ensure a near-uniformity of concentration. Note also that if you have a line like the ~nai line that converts currents to changes in concentration, you cannot include another mechanism that does the same thing. NEURON by default updates reversal potentials (here, ena) for you whenever the corresponding concentrations (here, nai) are updated in a mod file.

Anyone can help me to do Na/K restoration process, which pushes the concentrations back to equilibrium, because right now the sodium and potassium go up at a certain time and stay constant. is there mechanism to push them back to the resting value.

Actually the basic equation is conservation of mass. Rewriting in english
rate of change of total mass of solute S in compartment C = (sum of all fluxes of S into C)/(volume of C)
Diffusion is only one of many possible mechanisms that can generate a flux of S.

right now the sodium and potassium go up at a certain time and stay constant. is there mechanism to push them back to the resting value