A new stable and accurate numerical method is used to solve the equations with Neumann boundaries in cylindrical coordinates. The simulation analysis agrees with experimental measurements of end-plate current and agrees well with the results of the conformational state of the acetylcholine receptor as a function of time and acetylcholine concentration of earlier investigations with a smaller error compared to experiments. An asymmetric emission of acetylcholine in the synaptic cleft and the subsequent effects on open receptor population is simulated. Sensitivity of the open receptor dynamics...