We revisit the American put and call option valuation problems. We derive analytical formulas for the option prices and approximate ordinary differential equations for the optimal exercise boundaries. Numerical simulations yield accurate option prices and comparable computational speeds when benchmarked against the binomial method for calculating option prices. Our approach is based on the Mellin transform and an adaptation of the Kármán–Pohlausen technique for boundary layers in fluid mechanics.
