- About this Journal
- Abstracting and Indexing
- Aims and Scope
- Annual Issues
- Article Processing Charges
- Articles in Press
- Author Guidelines
- Bibliographic Information
- Citations to this Journal
- Contact Information
- Editorial Board
- Editorial Workflow
- Free eTOC Alerts
- Publication Ethics
- Reviewers Acknowledgment
- Submit a Manuscript
- Subscription Information
- Table of Contents
Abstract and Applied Analysis
Volume 2013 (2013), Article ID 120849, 10 pages
Combined Heat and Power Dynamic Economic Dispatch with Emission Limitations Using Hybrid DE-SQP Method
1Department of Mathematics, Faculty of Science, King Abdulaziz University, P.O. Box 80203, Jeddah 21589, Saudi Arabia
2Department of Mathematics, Faculty of Science, Al-Azhar University, Assiut 71511, Egypt
3Centre of New Energy Systems, Department of Electrical, Electronic and Computer Engineering, University of Pretoria, Pretoria 0002, South Africa
Received 28 August 2013; Accepted 1 October 2013
Academic Editor: Jinde Cao
Copyright © 2013 A. M. Elaiw et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Combined heat and power dynamic economic emission dispatch (CHPDEED) problem is a complicated nonlinear constrained multiobjective optimization problem with nonconvex characteristics. CHPDEED determines the optimal heat and power schedule of committed generating units by minimizing both fuel cost and emission simultaneously under ramp rate constraints and other constraints. This paper proposes hybrid differential evolution (DE) and sequential quadratic programming (SQP) to solve the CHPDEED problem with nonsmooth and nonconvex cost function due to valve point effects. DE is used as a global optimizer, and SQP is used as a fine tuning to determine the optimal solution at the final. The proposed hybrid DE-SQP method has been tested and compared to demonstrate its effectiveness.
Recently, combined heat and power (CHP) units, known as cogeneration or distributed generation, have played an increasingly important role in the utility industry. CHP units can provide not only electrical power but also heat to the customers. While the efficiency of the normal power generation is between 50% and 60%, the power and heat cogeneration increases the efficiency to around 90% . Besides thier high efficiency, CHP units reduce the emission of gaseous pollutants (SO2, , CO, and) by about 13–18% .
In order to utilize the integrated CHP system more CO2 economically, combined heat and power economic dispatch (CHPED) problem is applied. The objective of the CHPED problem is to determine both power generation and heat production from units by minimizing the fuel cost such that both heat and power demands are met, while the combined heat and power units are operated in a bounded heat versus power plane. For most CHP units the heat production capacities depend on the power generation. This mutual dependency of the CHP units introduces a complication to the problem . In addition, considering valve point effects in the CHPED problem makes the problem nonsmooth with multiple local optimal point which makes finding the global optimal challenging.
In the literature, several optimization techniques have been used to solve the CHPED problem with complex objective functions or constraints such as Lagrangian relaxation (LR) [4, 5], semidefinite programming (SDP) , augmented Lagrange combined with Hopfield neural network , harmony search (HS) algorithm [1, 8], genetic algorithm (GA) , ant colony search algorithm (ACSA) , mesh adaptive direct search (MADS) algorithm , self adaptive real-coded genetic algorithm (SARGA) , particle swarm optimization (PSO) [2, 12], artificial immune system (AIS) , bee colony optimization (BCO) , differential evolution , and evolutionary programming (EP) . In [2, 13–15], the valve point effects and the transmission line losses are incorporated into the CHPED problem.
In the CHPED formulation the ramp rate limits of the units are neglected. Plant operators, to avoid life-shortening of the turbines and boilers, try to keep thermal stress on the equipments within the safe limits. This mechanical constraint is usually transformed into a limit on the rate of change of the electrical output of generators. Such ramp rate constraints link the generator operation in two consecutive time intervals. Combined heat and power dynamic economic dispatch (CHPDED) problem is an extension of CHPED problem where the ramp rate constraint is considered. The primary objective of the CHPDED problem is to determine the heat and power schedule of the committed units so as to meet the predicted heat and electricity load demands over a time horizon at minimum operating cost under ramp rate constraints and other constraints . Since the ramp rate constraints couple the time intervals, the CHPDED problem is a difficult optimization problem. If the ramp rate constraints are not included in the optimization problem, the CHPDED problem is reduced to a set of uncoupled CHPED problems that can easily be solved. In the literature an overwhelming number of reported works deal with CHPED problem; however, the CHPDED problem has only been considered in .
The traditional dynamic economic dispatch (DED) problem which considers only thermal units that provide only electric power has been studied by several authors (see the review paper ). The emission has been taken into the traditional (DED) formulation in three main approaches. The first approach is to minimize the fuel cost and treat the emission as a constraint with a permissible limit (see, e.g., [19–21]). This formulation, however, has a severe difficulty in getting the trade-off relations between cost and emission . The second approach handles both fuel cost and emission simultaneously as competing objectives [23–25]. The third approach treats the emission as another objective in addition to fuel cost objective. However, the multiobjective optimization problem is converted to a single-objective optimization problem by linear combination of both objectives [19, 26–30]. In the second and third approaches, the dynamic dispatch problem is referred to as dynamic economic emission dispatch (DEED) which is a multiobjective optimization problem, which minimizes both fuel cost and emission simultaneously under ramp rate constraint and other constraints [19, 24]. In this paper, we incoroporate the CHP units into the DEED problem. Combined heat and power dynamic economic emission dispatch (CHPDEED) is formulated with the objective to determine the unit power and heat production so that the system’s production cost and emission are simultaneously minimized, while the power and heat demands and other constraints are met . The emission has been taken into consideration in the CHPED and CHPDED in [17, 31], respectively. In , both fuel cost and emission are simultaneously handled as competing objectives and the multiobjective problem is solved using an enhanced firefly algorithm (FA). In the present paper, the multiobjective optimization problem is converted into a single-objective optimization using the weighting method. This approach yields meaningful result to the decision maker when solved many times for different values of the weighting factor. In , the simulation results for test system are shown, but the data of the heat demand is not explicitly tabulated; instead it is expressed graphically (see Figure 12 in ). In this case a comparison of our proposed method and FA cannot be performed. In our paper, all the data and the solutions of the test system are available for comparison.
Differential evolution algorithm (DE), which was proposed by Storn and Price  is a population based stochastic parallel search technique. DE uses a rather greedy and less stochastic approach to problem solving compared to other evolutionary algorithms. DE has the ability to handle optimization problems with nonsmooth/nonconvex objective functions . Moreover, it has a simple structure and a good convergence property, and it requires a few robust control parameters . DE has been applied to the CHPED and CHPDED problems with non-smooth and non-convex cost functions in [15, 33], respectively.
The DE shares many similarities with evolutionary computation techniques such as genetic algorithms (GA) techniques. The system is initialized with a population of random solutions and searches for optima by updating generations. DE has evolution operators such as crossover and mutation. Although DE seem to be good methods to solve the CHPDEED problem with non-smooth and non-convex cost functions, solutions obtained are just near global optimum with long computation time. Therefore, hybrid methods such as DE-SQP can be effective in solving the CHPDEED problems with valve point effects.
The main contributions of the paper are as follows. (1) A multi-objective optimization problem is formulated using CHPDEED approach. The multi-objective optimization problem is converted into a single-objective optimization using the weighting method. (2) Hybrid DE-SQP method is proposed and validated for solving the CHPDEED problem with nonsmooth and nonconvex objective function. DE is used as a base level search for global exploration and SQP is used as a local search to fine-tune the solution obtained from DE. (3) The effectiveness of the proposed method is shown for test systems.
2. Problem Formulation
In this section we formulate the CHPDEED problem. The system under consideration has three types of generating units, conventional thermal units (TU), CHP units, and heat-only units (H). The power is generated by conventional thermal units and CHP units, while the heat is generated by CHP units and heat-only units. The objective of the CHPDEED problem is to simultaneously minimize the system’s production cost and emission so as to meet the predicted heat and power load demands over a time horizon under ramp rate and other constraints. The following objectives and constraints are taken into account in the formulation of the CHPDEED problem.
2.1. Objective Functions
In this section, we introduce the cost and emission functions of three types of generating units, conventional thermal units which produce power only, CHP units which produce both heat and power, and heat-only units which produce heat only.
2.1.1. Conventional Thermal Units
Cost. The cost function curve of a conventional thermal unit can be approximated by a quadratic function . Power plants commonly have multiple valves which are used to control the power output of the unit. When steam admission valves in conventional thermal units are first open, a sudden increase in losses is registered which results in ripples in the cost function [18, 36]. This phenomenon is called as valve-point effects. The generator with valve-point effects has very different input-output curve compared with smooth cost function. Taking the valve-point effects into consideration, the fuel cost is expressed as the sum of a quadratic and sinusoidal functions [17, 24, 25, 37]. Therefore, the fuel cost function of the conventional thermal units is given by where , and are positive constants, and are the coefficients of conventional thermal unit reflecting valve-point effects, is the power generation of conventional thermal unit during the th time interval , is the minimum capacity of conventional thermal unit , and is the fuel cost of conventional thermal unit to produce .
Emission. The amount of emission of gaseous pollutants from conventional thermal units can be expressed as a combination of quadratic function and exponential function of the unit’s active power output . The emission function is given by where is the amount of emission from unit from producing power . Constants , , , , and are the coefficients of the th unit emission characteristics .
2.1.2. CHP Units
Cost. A CHP unit has a convex cost function in both power and heat. The form of the fuel cost function of CHP units can be given by [6, 17] the following: where is the generation fuel cost of CHP unit to produce power and heat . Constants , and are the fuel cost coefficients of CHP unit .
2.1.3. Heat-Only Units
Let be the number of dispatch intervals and the number of committed units, where is the number of conventional thermal units, is the number of the CHP units, and is the number of the heat-only units. Then the total fuel cost and amount of emission over the dispatch period are given, respectively, by where , , , , , and .
There are three kinds of constraints considered in the CHPDEED problem, that is, the equilibrium constraints of power and heat production, the capacity limits of each unit, and the ramp rate limits.
(i) Power Production and Demand Balance where and are the system power demand and transmission line losses at time (i.e., the th time interval), respectively. The B-coefficient method is one of the most commonly used by power utility industry to calculate the network losses. In this method the network losses are expressed as a quadratic function of the unit’s power outputs that can be approximated in the following: where and is the th element of the loss coefficient square matrix of size .
(ii) Heat Production and Demand Balance where is the system heat demand at time .
(iii) Capacity Limits of Conventional Thermal Units where and are the minimum and maximum power capacity of conventional thermal unit , respectively.
(iv) Capacity Limits of CHP Units where and are the minimum and maximum power limit of CHP unit , respectively, and they are functions of generated heat (. and are the heat generation limits of CHP unit which are functions of generated power .
(v) Capacity Limits of Heat-Only Units where and are the minimum and maximum heat capacity of heat-only unit , respectively.
(vi) Upper/Down Ramp Rate Limits of Conventional Thermal Units where and are the maximum ramp up/down rates for conventional thermal unit .
(vii) Upper/Down Ramp Rate Limits of CHP Units where and are the maximum ramp up/down rates for CHP unit .
2.3. The Optimization Problem
Aggregating the objectives and constraints, the CHPDEED problem can be mathematically formulated as a nonlinear constrained multi-objective optimization problem which can be converted into a single-objective optimization using the weighting method as where is a weighting factor. It will be noted that, when , problem (17) determines the optimal amount of the generated heat and power by minimizing the fuel cost regardless of emission and the problem will be referred to as combined heat and power dynamic economic dispatch (CHPDED) problem. If , then problem (17) determines the optimal amount of the generated power by minimizing the emission regardless of cost and the problem will be referred to as combined heat and power pure dynamic emission dispatch (CHPPDED).
3. Differential Evolution Method
DE is a simple yet powerful heuristic method for solving nonlinear, nonconvex, and nonsmooth optimization problems. DE algorithm is a population based algorithm using three operators; mutation, crossover, and selection to evolve from randomly generated initial population to final individual solution . In the initialization a population of NP target vectors (parents) , , is randomly generated within user-defined bounds, where is the dimension of the optimization problem. Let be the individual at the current generation . A mutant vector is generated according to with randomly chosen integer indexes . Here is the mutation factor.
According to the target vector and the mutant vector , a new trial vector (offspring) is created with where , and is the evaluation of a uniform random number between . is the crossover constant which has to be determined by the user. is a randomly chosen index from which ensures that gets at least one parameter from .
The selection process determines which of the vectors will be chosen for the next generation by implementing one-to-one competition between the offsprings and their corresponding parents. If denotes the function to be minimized, then where . The value of of each trial vector is compared with that of its parent target vector . The above iteration process of reproduction and selection will continue until a user-specified stopping criteria is met.
In this paper, we define the evaluation function for evaluating the fitness of each individual in the population in DE algorithm as follows: where and are penalty values. Then the objective is to find , the minimum evaluation value of all the individuals in all iterations. The penalty term reflects the violation of the equality constraints. Once the minimum of is reached, the equality constraints are satisfied.
4. Sequential Quadratic Programming Method
SQP method can be considered as one of the best nonlinear programming methods for constrained optimization problems . It outperforms every other nonlinear programming method in terms of efficiency, accuracy, and percentage of successful solutions over a large number of test problems. The method closely resembles Newton’s method for constrained optimization, just as is done for unconstrained optimization. At each iteration, an approximation is made of the Hessian of the Lagrangian function using Broyden-Fletcher-Goldfarb-Shanno (BFGS) quasi-Newton updating method. The result of the approximation is then used to generate a quadratic programming (QP) subproblem whose solution is used to form a search direction for a line search procedure. Since the objective function of the CHPDEED problem is non-convex and non-smooth, SQP ensures a local minimum for an initial solution. In this paper, DE is used as a global search and finally the best solution obtained from DE is given as initial condition for SQP method as a local search to fine-tune the solution. SQP simulations can be computed by the fmincon code of the MATLAB Optimization Toolbox.
5. Simulation Results
In this section we present two examples. The first example shows the efficiency of the proposed DE-SQP method for the DED problem. In the second example, the hybrid DE-SQP method is applied to the CHPDEED problem. In DE-SQP method, the control parameters are chosen as , and . The maximum number of iterations are selected as . The results represent the average of runs of the proposed method. All computations are carried out by MATLAB program.
Example 1. This example consists of ten conventional thermal units to investigate the effectiveness of the proposed DE-SQP technique in solving the DED problem with valve point effects and transmission line losses. The technical data of the units as well as the demand for the 10-unit system are taken from . The best solution of the DED problem is given in Table 1. Comparison between our proposed method (DE-SQP) and other methods is given in Table 2. It is observed that the proposed method reduces the total generation cost better than the other methods reported in the literature.
Example 2. This example is 11-unit system (eight conventional thermal units, two CHP units, and one heat-only unit) for solving the CHPDED, CHPDEED, and CHPPDED problems using DE-SQP method. We shall solve the CHPDEED problem when , in addition to the CHPDED and CHPPDED problems which correspond to and , respectively. The technical data of conventional thermal units, the matrix , and the demand are taken from the 10-unit system presented in . The 5th and 8th conventional units in  were replaced by two CHP units. The technical data of the two CHP units and the heat-only unit are taken from  and are given in Table 3. The heat demand for 24 hours is given in Table 4. The feasible operating regions of the two CHP units are given in Figures 1 and 2 (see [4, 14]).
The best solutions of the CHPDED, CHPDEED, and CHPPDED problems for DE-SQP algorithm are given in Tables 5, 6, and 7, respectively. The best cost, the amount of emission, and the transmission line losses are also given in Tables 5–7. It is seen that the cost is under CHPDED, but it increases to under CHPPDED. The emission obtained from CHPDED is , but it decreases to under CHPPDED. Under the CHPDEED problem, the cost is which is more than and less than . Moreover, the emission is which is less than and more than .
This paper presents a hybrid method combining differential evolution (DE) and sequential quadratic programming (SQP) for solving dynamic dispatch (CHPDED, CHPDEED, and CHPPDED) problems with valve-point effects including generator ramp rate limits. In this paper, DE is first applied to find the best solution. This best solution is given to SQP as an initial condition that fine tunes the optimal solution at the final. The feasibility and efficiency of the DE-SQP were illustrated by conducting case studies with system consisting of eight conventional thermal units, two CHP units, and one heat-only unit.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
This work was funded by the Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah, under Grant no. (130-107-D1434). The authors, therefore, acknowledge with thanks DSR technical and financial support.
- A. Vasebi, M. Fesanghary, and S. M. T. Bathaee, “Combined heat and power economic dispatch by harmony search algorithm,” International Journal of Electrical Power and Energy Systems, vol. 29, no. 10, pp. 713–719, 2007.
- M. Behnam, M. Mohammad, and R. Abbas, “Combined heat and power economic dispatch problem solution using particle swarm optimization with time varying acceleration coefficients,” Electric Power Systems Research, vol. 95, pp. 9–18, 2013.
- P. Subbaraj, R. Rengaraj, and S. Salivahanan, “Enhancement of combined heat and power economic dispatch using self adaptive real-coded genetic algorithm,” Applied Energy, vol. 86, no. 6, pp. 915–921, 2009.
- T. Guo, M. I. Henwood, and M. van Ooijen, “An algorithm for combined heat and power economic dispatch,” IEEE Transactions on Power Systems, vol. 11, no. 4, pp. 1778–1784, 1996.
- A. Sashirekha, J. Pasupuleti, N. H. Moin, and C. S. Tan, “Combined heat and power (CHP) economic dispatch solved using Lagrangian relaxation with surrogate subgradient multiplier updates,” Electrical Power and Energy Systems, vol. 44, pp. 421–430, 2013.
- A. M. Jubril, A. O. Adediji, and O. A. Olaniyan, “Solving the combined heat and power dispatch problem: a semi-definite programming approach,” Electric Power Components and Systems, vol. 40, pp. 1362–1376, 2012.
- V. N. Dieu and W. Ongsakul, “Augmented lagrangehopfield network for economic load dispatch with combined heat and power,” Electric Power Components and Systems, vol. 37, no. 12, pp. 1289–1304, 2009.
- E. Khorram and M. Jaberipour, “Harmony search algorithm for solving combined heat and power economic dispatch problems,” Energy Conversion and Management, vol. 52, no. 2, pp. 1550–1554, 2011.
- C. Su and C. Chiang, “An incorporated algorithm for combined heat and power economic dispatch,” Electric Power Systems Research, vol. 69, no. 2-3, pp. 187–195, 2004.
- Y. H. Song, C. S. Chou, and T. J. Stonham, “Combined heat and power economic dispatch by improved ant colony search algorithm,” Electric Power Systems Research, vol. 52, no. 2, pp. 115–121, 1999.
- S. S. Sadat Hosseini, A. Jafarnejad, A. H. Behrooz, and A. H. Gandomi, “Combined heat and power economic dispatch by mesh adaptive direct search algorithm,” Expert Systems with Applications, vol. 38, no. 6, pp. 6556–6564, 2011.
- V. Ramesh, T. Jayabarathi, N. Shrivastava, and A. Baska, “A novel selective particle swarm optimization approach for combined heat and power economic dispatch,” Electric Power Components and Systems, vol. 37, no. 11, pp. 1231–1240, 2009.
- M. Basu, “Artificial immune system for combined heat and power economic dispatch,” Electrical Power and Energy Systems, vol. 43, pp. 1–5, 2012.
- M. Basu, “Bee colony optimization for combined heat and power economic dispatch,” Expert Systems with Applications, vol. 38, no. 11, pp. 13527–13531, 2011.
- M. Basu, “Combined heat and power economic dispatch by using differential evolution,” Electric Power Components and Systems, vol. 38, no. 8, pp. 996–1004, 2010.
- K. P. Wong and C. Algie, “Evolutionary programming approach for combined heat and power dispatch,” Electric Power Systems Research, vol. 61, no. 3, pp. 227–232, 2002.
- N. Taher, A. A. Rasoul, R. Alireza, and A. Babak, “A new multi-objective reserve constrained combined heat and power dynamic economic emission dispatch,” Energy, vol. 42, pp. 530–545, 2012.
- X. Xia and A. M. Elaiw, “Optimal dynamic economic dispatch of generation: a review,” Electric Power Systems Research, vol. 80, no. 8, pp. 975–986, 2010.
- A. M. Elaiw, X. Xia, and A. M. Shehata, “Application of model predictive control to optimal dynamic dispatch of generation with emission limitations,” Electric Power Systems Research, vol. 84, no. 1, pp. 31–44, 2012.
- G. P. Granelli, M. Montagna, G. L. Pasini, and P. Marannino, “Emission constrained dynamic dispatch,” Electric Power Systems Research, vol. 24, no. 1, pp. 55–64, 1992.
- Y. H. Song and I. Yu, “Dynamic load dispatch with voltage security and environmental constraints,” Electric Power Systems Research, vol. 43, no. 1, pp. 53–60, 1997.
- M. A. Abido, “Environmental/economic power dispatch using multiobjective evolutionary algorithms,” IEEE Transactions on Power Systems, vol. 18, no. 4, pp. 1529–1537, 2003.
- M. Basu, “Dynamic economic emission dispatch using evolutionary programming and fuzzy satisfied method,” International Journal of Emerging Electric Power Systems, vol. 8, pp. 1–15, 2007.
- M. Basu, “Dynamic economic emission dispatch using nondominated sorting genetic algorithm-II,” International Journal of Electrical Power and Energy Systems, vol. 30, no. 2, pp. 140–149, 2008.
- C. X. Guo, J. P. Zhan, and Q. H. Wu, “Dynamic economic emission dispatch based on group search optimizer with multiple producers,” Electric Power Systems Research, vol. 86, pp. 8–16, 2012.
- A. M. Elaiw, X. Xia, and A. M. Shehata, “Hybrid DE-SQP and hybrid PSO-SQP methods for solving dynamic economic emission dispatch problem with valve-point effects,” Electric Power Systems Research, vol. 84, pp. 192–200, 2013.
- A. M. Elaiw, X. Xia, and A. M. Shehata, “Minimization of fuel costs and gaseous emissions of electric power generation by model predictive control,” Mathematical Problems in Engineering, vol. 2013, Article ID 906958, 15 pages, 2013.
- J. S. Alsumait, M. Qasem, J. K. Sykulski, and A. K. Al-Othman, “An improved Pattern Search based algorithm to solve the Dynamic Economic Dispatch problem with valve-point effect,” Energy Conversion and Management, vol. 51, no. 10, pp. 2062–2067, 2010.
- M. Basu, “Particle swarm optimization based goal-attainment method for dynamic economic emission dispatch,” Electric Power Components and Systems, vol. 34, no. 9, pp. 1015–1025, 2006.
- P. Nicole, T. Anshul, T. Shashikala, and P. Manjaree, “An improved bacterial foraging algorithm for combined static/dynamic environmental economic dispatch,” Applied Soft Computing, vol. 12, pp. 3500–3513, 2012.
- G. S. Piperagkas, A. G. Anastasiadis, and N. D. Hatziargyriou, “Stochastic PSO-based heat and power dispatch under environmental constraints incorporating CHP and wind power units,” Electric Power Systems Research, vol. 81, no. 1, pp. 209–218, 2011.
- R. Storn and K. Price, “Differential Evolution—a simple and efficient adaptive scheme for global optimization over continuous spaces,” Journal of Global Optimization, vol. 11, no. 4, pp. 341–359, 1997.
- A. M. Elaiw, X. Xia, and A. M. Shehata, “Hybrid DE-SQP method for solving combined heat and power dynamic economic dispatch problem,” Mathematical Problems in Engineering, vol. 2013, Article ID 982305, 7 pages, 2013.
- M. Basu, “Artificial immune system for dynamic economic dispatch,” International Journal of Electrical Power and Energy Systems, vol. 33, no. 1, pp. 131–136, 2011.
- X. Xia, J. Zhang, and A. Elaiw, “An application of model predictive control to the dynamic economic dispatch of power generation,” Control Engineering Practice, vol. 19, no. 6, pp. 638–648, 2011.
- C. K. Panigrahi, P. K. Chattopadhyay, R. N. Chakrabarti, and M. Basu, “Simulated annealing technique for dynamic economic dispatch,” Electric Power Components and Systems, vol. 34, no. 5, pp. 577–586, 2006.
- P. Attaviriyanupap, H. Kita, E. Tanaka, and J. Hasegawa, “A hybrid EP and SQP for dynamic economic dispatch with nonsmooth fuel cost function,” IEEE Transactions on Power Systems, vol. 17, no. 2, pp. 411–416, 2002.
- P. T. Boggs and J. W. Tolle, “Sequential quadratic programming,” Acta Numerica, vol. 3, no. 4, pp. 1–52, 1995.