Numerical calculation of an asymmetric supply function equilibrium with capacity constraints

Producers submit offer curves to a procurement auction, e.g. an electricity auction, before uncertain demand has been realised. In the supply function equilibrium (SFE), every firm commits to the offer curve that maximises its expected profit, given the offer curves of competitors. The equilibrium is given by a system of differential equations. In practice, it has been very difficult to find valid SFE, i.e. non-decreasing solutions, from this system, especially for asymmetric producers. This paper shows that valid SFE can be calculated by means of a shooting algorithm that combines numerical integration with an optimisation procedure that searches for an end-condition. Multiple/parallel shooting is used for ill-conditioned cases.