Physics based click to move

I made a method for this. Just be sure that all variables use the same time step (rigidbody.velocity for instance is in seconds while you apply your forces every fixed update which is 50 times per second).

I know this is a old question but iam sure you are not the only one wanting to do this.

         /// <summary>
         /// Tries to minimize the error. Useful when working with physics and trying to arrive at a position and similar problems.
         /// A example of such a problem would be a car that has to move to a goal and has to arrive exactly at that goal without overshooting.
         /// </summary>
         /// <param name="error"></param>The error. This is a value that indicates how much we are off our goal. Think of this as the distance to your goal
         /// <param name="deltaStep"></param>How much the error will change in the next step. Think of this as the current speed of your car.
         /// <param name="maxDeltaDelta"></param>The max amount of change that we can apply to deltaStep. Think of this as the acceleration of your car.
         /// <returns></returns>
         public static float CalculateArrival(float error, float deltaStep, float maxDeltaDelta)
         {
             double goalDeltaStep;
             var minInfluenceRange = deltaStep - maxDeltaDelta;
             var maxInfluenceRange = deltaStep + maxDeltaDelta;
             if(error > minInfluenceRange && error < maxInfluenceRange)
             {
                 goalDeltaStep = -error; //This prevents overstepping
             }
             else
             {
                 if (error >= 0)
                 {
                     goalDeltaStep = -Math.Sqrt(Math.Abs(2f * maxDeltaDelta * error));
 
                 }
                 else
                 {
                     goalDeltaStep = Math.Sqrt(Math.Abs(2f * maxDeltaDelta * error));
                 }
             }
             var difference = goalDeltaStep - deltaStep;
             return Mathf.Clamp((float)difference, -maxDeltaDelta, maxDeltaDelta);
         }