Very high-magnification microlensing events provide chances to measure limb darkening of distant stars. We use the Finite Element Method (FEM) as an inversion tool for discretization and inversion of the magnification-limb darkening integral equation. This method makes no explicit assumption about the shape of brightness profile more than the flatness of the profile near the centre of the stellar disk. From the simulation, we investigate the accuracy and stability of this method and we use regularization techniques to stabilize it. Finally, we apply this method to the single lens, high magnification transit events of OGLE-2004-BLG-254 (SAAO I ),