Pricing dericative securities in incomplete markets

We propose the algorithms for pricing American and European options in incomplete markets. We consider a non-self-financing replicating portfolio and minimize the hedging error consisting of the self-financing error of the portfolio dynamics and the error of the option's payoff replication. We treat the pricing problem as regression with constraints and reduce it to a quadratic minimization problem. The algorithms of pricing American and European options differ in imposing one additional type of constraints. Prices of options for different initial and strike prices can be found in one optimization run. The algorithms create a table representing the option price as a function of time and the underlying stock price for the whole lifetime of the option. We illustrate the numerical performance of the algorithms with options on futures contracts in natural gas market.