Acessibilidade / Reportar erro

The single machine earliness and tardiness scheduling problem: lower bounds and a branch-and-bound algorithm

Abstract

This paper addresses the single machine scheduling problem with a common due date aiming to minimize earliness and tardiness penalties. Due to its complexity, most of the previous studies in the literature deal with this problem using heuristics and metaheuristics approaches. With the intention of contributing to the study of this problem, a branch-and-bound algorithm is proposed. Lower bounds and pruning rules that exploit properties of the problem are introduced. The proposed approach is examined through a computational comparative study with 280 problems involving different due date scenarios. In addition, the values of optimal solutions for small problems from a known benchmark are provided. Mathematical subject classification: 90C11, 62P30, 90B35.

single machine; common due date; earliness and tardiness; lower bound; branch-and-bound


The single machine earliness and tardiness scheduling problem: lower bounds and a branch-and-bound algorithm* * This research has been financially supported by "Funda ção de Amparo à Pesquisa do Estado de São Paulo" - FAPESP (Grant 06/03496-3) and by "Conselho Nacional de Desenvolvimento Científico e Tecnológico" - CNPq (Grants 486124/2007-0 and 307399/2006-0).

Débora P. Ronconi; Márcio S. Kawamura

Department of Production Engineering, EPUSP, University of São Paulo Av. Prof. Almeida Prado, 128, Cidade Universitária, 05508-900, São Paulo SP, Brazil E-mails: dronconi@usp.br / marcio.kawamura@poli.usp.br

ABSTRACT

This paper addresses the single machine scheduling problem with a common due date aiming to minimize earliness and tardiness penalties. Due to its complexity, most of the previous studies in the literature deal with this problem using heuristics and metaheuristics approaches. With the intention of contributing to the study of this problem, a branch-and-bound algorithm is proposed. Lower bounds and pruning rules that exploit properties of the problem are introduced. The proposed approach is examined through a computational comparative study with 280 problems involving different due date scenarios. In addition, the values of optimal solutions for small problems from a known benchmark are provided.

Mathematical subject classification: 90C11, 62P30, 90B35.

Key words: single machine, common due date, earliness and tardiness, lower bound, branch-and-bound.

1 Introduction

Scheduling problems involving both earliness and tardiness costs have received significant attention in recent years. This type of problem became more important with advent of lean production principles, including the just-in-time (JIT) concept. According to JIT, earliness and tardiness are considered harmful to profitability and, for this reason, must be minimized: tardiness causes loss of customer goodwill and damage reputation, as well as delay of payments, while earliness causes inventory carrying costs and possible loss of product quality. Probably based on this motivation, many authors have considered the scheduling problem aiming to minimize earliness and tardiness in the delivery of goods. Comprehensive surveys on the common due date assignment and scheduling problems can be found in [8, 2].

The problem of scheduling jobs with a common due date in a single machine has been studied by several authors. From a practical point of view, accordingto [7], customer orders with a combination of goods to be delivered at a specified time, export shipping and chemical or physical mixtures containing some ingredients with a short half-life period are some examples where jobs are to be scheduled in a single machine and be delivered on a common due date.

In this problem, there are n jobs available at time zero to be processed on a single machine and to be delivered on a common due date d. Each job i requires exactly one operation and its processing time pi is known. If a job i is completed before the due date, its earliness is given by Ei = d - Ci, where Ci is the completion time of job i. Conversely, if a job i is completed after the desired date, its tardiness is given by Ti = Ci - d. Each job i has its own unit earliness penalty αi and unit tardiness penalty βi. Preemption is not allowed and the initial processing time is not necessarily at time zero, when all jobs are available. The objective of the problem is to obtain an optimal schedule that minimizes the sum of earliness and tardiness penalties.

The common due date can be restrictive or unrestrictive. A due date is called unrestrictive if its optimal value has to be calculated or if its given value does not influence the optimal schedule. In the cases where the given due date is greater than or equal to the sum of processing times of all jobs available, this due date is unrestrictive [7]. In [9] it was demonstrated that this scheduling problem is NP-hard even with the unrestrictive due date and αi = βi. In [15] the author addressed the particular unrestrictive case in which αi = βi = 1 for all jobs, that can be solved by a polynomial algorithm of O(n log n) complexity. The case in which penalties are independent of the jobs (αi = α and βi = β for all i) can also be treated by polynomial algorithms [21]. For the general case (no restrictions on the penalties), a branch-and-bound algorithm that is capable of solving instances with up to 15 jobs was proposed in [6]. By that time, a 0-1 quadratic model for this problem was presented and solved with a specialized branch-and-bound algorithm [5]. Another exact approach that combined column generation with a Lagrangean relaxation algorithm was presented in [1]. This strategy solves problems with up to 125 jobs.

The restrictive version of this problem is NP-hard even with αi = βi = 1 [10, 12]. Due to its complexity, many authors addressed this problem using heuristic and metaheuristic approaches (see, for example, [13, 17, 7, 11, 18]). Most of these methods based their search strategies on the following properties.

For the restrictive common due date case with general penalties, there is an optimal solution with these properties:

1. No idle times are inserted between consecutive jobs [4];

2. The schedule is V-Shaped, that is, jobs that complete on or before the due date are sequenced in a non-increasing order of the pi / αi ratio. The jobs that start on or after the due date are sequenced in a non-decreasing order of the pi / βi ratio. Note that there may be a straddling job, i.e., a job whose processing is started before and finished after the due date (see [12, 3]);

3. There is an optimal schedule in which either the processing time of the first job starts at time zero or one job is completed on the due date. The proof is similar to the one presented in [12].

The development of exact algorithms for the restrictive common due date case with special characteristics was considered by some authors. A pseudopolynomial dynamic programming algorithm for this problem with αi = βi for all jobs was presented in [12], while a pseudopolynomial algorithm for the constraint problem (machine idle time is forbidden) with αi = α and βi = β was proposed by [14]. In [19] the constraint problem with the general penalties case and different due dates related to each job was considered. The author presented a branch-and-bound algorithm that makes use of a lower bound based on Lagrangean relaxation. More recently, algorithms based on dynamic programming and branch-and-bound schemes for the unconstrained problem with uniform penalties α and β were presented [20].

Aiming to contribute with the study of this problem with restrictive common due date and general penalties, this paper addresses the development of a specific branch-and-bound algorithm. Lower bounds and pruning rules that exploit the properties of the problem are introduced. Computational tests are presented and the performance of the proposed algorithm is analysed through a comparative study with 280 problems involving different due date scenarios. In addition, the values of optimal solutions for small problems from a known benchmark are provided. Similar approaches were successfully applied in [22] and [23] to the flowshop environment with blocking aiming to minimize the makespan and the tardiness criteria, respectively.

This paper is organized as follows: the next section presents a mathematical model. Section 3 describes the search strategies of the branch-and-bound algorithm and the proposed lower bounds. Section 4 shows the computational results, while the last section summarizes the main results.

Notation. For all v we denote v+ = max{v, 0}.

2 Mathematical model

The following mixed integer linear programming (MILP) formulation, slightly modified from the model presented in [3], can be used to obtain optimal solutions to this problem.

Parameters:

d: common due date;

αi: earliness penalty of job i per time unit;

βi: tardiness penalty of job i per time unit;

pi: processing time of job i;

R: sufficiently large number.

Variables:

xik: 1, if job i is sequenced (not necessarily directly) prior to job k. 0, otherwise;

Ci: completion time of job i;

Ei: earliness of job i;

Ti: tardiness of job i.

Model:

subject to

Equation (1) represents the objective function to be minimized, i.e., the sum of tardiness and earliness penalties. In [3] the tardiness and earliness are calculated through the following restrictions:

Ti>d - Ci, i = 1, 2,..., n,

Ei>Ci - d, i = 1, 2,..., n.

Alternatively, we calculate tardiness and earliness through constraints ofform (2). Note that the presented model has n restrictions less than the original one. Constraints of form (3) and (4) indicate the completion time of each job: if job i is sequenced prior to job k, xik = 1 and, consequently, restriction (3) gives Ci < Ck - pk and, due the addition of constant R, restriction (4) is not restrictive. On the other hand, if xik = 0, restriction (4) becomes Ck < Ci - pi and restriction (3) is not restrictive. Restriction (5) assures that initial time of each job i is not negative. The set of restrictions (6) and (7) defines the non-negativity of variables Ti and Ei, while restriction (8) defines the variable xik as binary.

3 Branch-and-bound algorithm

The proposed branch-and-bound algorithm is composed of two diverse strategies. This separation is based on Property 3 (see Section 1), which states that there is an optimal solution in which either the processing of the first job starts at time zero or one job is completed on the due date.

In the first strategy (Search Strategy 1), schedules that have a job being completed exactly on the due date are explored, while in the second one (Search Strategy 2), schedules that start at time zero are investigated. These strategies will be applied in a sequential form and their combined execution guarantees that the algorithm covers the entire solution space. The initial incumbent solution is provided by a constructive heuristic, HRM, proposed in [11].

3.1 Search Strategy 1

This search strategy looks for the best solution that has a job being completed exactly on the due date; the starting time of this solution may be different from time zero. The sequence of jobs is treated as two subsequences: in one subsequence, the last job finishes on the due date and, in the other, the first job starts on the due date.

The first level of the search tree corresponds to the number of jobs that finish their processing before or on the due date (nfb). For the remaining levels, only two nodes will be generated representing the relative position of a job in relation to the due date. These positions are indicated by the label Before (the analyzed job finishes its processing before or on the due date) or After (the job starts its processing after or on the due date). It was assumed that level 2 corresponds to job 1; level 3 corresponds to job 2, and so on. When the node Before is created, the algorithm checks if there is enough space to allocate the new job. The order of the jobs in each subsequence follows the V-shape (Property 2).

Figure 1 illustrates an example of the proposed search tree with four jobs.Note that, in the first level, solutions with nfb = 4 and nfb = 0 are not part of the tree. If feasible schedules with all non-tardy jobs can be obtained, the considered due date is unrestrictive. On the other hand, a schedule with all late jobs always can be improved with a left shift of the sequence.


In the first level of the tree, the father node will be the one with the largest number of non-tardy jobs. This criterion aims to foster the elimination of solutions that have the minimum sum of nfb processing times greater than the period of time available before the due date. The subtree below the chosen node will be investigated before the selection of another node in this level.

For each father node in the first level, a lower bound is calculated. The proposed lower bound, LB1a, is presented in the following proposition:

Proposition 1. Consider a set J of n jobs that is composed of two subsequences, one with n fb unknown jobs starting before the due date, and another with n-n fb unknown jobs completed after the due date. Assume that there is no idle time between consecutive jobs and the job in position n fb is completed on the due date. Then, the weighted sum of earliness and tardiness penalties of the jobs is greater than or equal to:

where

is the q-th smallest processing time among jobs in J.

Proof. Consider a known sequence of a set J with n jobs where nfb jobs are completed before or on the due date and n - nfb jobs are completed afterthe due date. The weighted sum of earliness and tardiness penalties of this sequence is given by:

where π(y) is the job that is allocated in position y.

It can be observed that:

Since LB1a is smaller than or equal to the weighted sum of earliness and tardiness penalties of an arbitrary sequence defined as in the statement of the proposition, we conclude that the thesis holds.

In the next levels, the node to be branched is the one that contains the subsequence that is closest to be completed. This proximity is measured by the smallest number of jobs to be positioned in each subsequence. It is calculated for each node using the following expression:

proximity = min {nfb - xb, n - nfb - xa},

where xb and xa are the number of jobs already positioned before and after the due date, respectively. Note that, if all jobs in any side of the due date have been fixed, the sequence of the remaining jobs is also determined by Property 2. In case of ties, the algorithm selects the node with the maximum lower bound LB1b which considers the jobs already positioned and it is evaluated according to the following expression:

LB1b = max{LB1a, LB'1b},

where

A node is fathomed if its lower bound is greater than or equal to the current incumbent solution value.

3.2 Search Strategy 2

This strategy addresses solutions whose initial time is zero. The algorithmbuilds a branch-and-bound (b&b) tree where each node represents a partial sequence. When a node is branched, one or more nodes are generated by adding one more job to the partial sequence associated with the node being branched. This sequence is constructed from end to beginning, i.e., schedule construction starts by positioning the last job to be processed and it continues fixing jobs adjacently until it reaches the first job to be processed. Note that the first fixed job will always finish its processing in time instant equivalent to the sum of all processing times, while the last fixed job will start at time zero. According to [16], this approach delivers better results than the natural one (constructingit from beginning to end). Figure 2 illustrates this construction order.


To reduce the search space, only schedules whose jobs are positioned in V-shape (Property 2) and without idle time between consecutive jobs (Property 1) are assembled.

The node that will be branched (father node) is the one with the largest sum of processing times of the already fixed jobs. This criterion aims to favor partial sequences whose starting time is closer to the due date. In case of ties, the algorithm chooses the node with the largest lower bound. When the starting time of a partial sequence is equal to or smaller than the due date, the best position of the remaining jobs is determined by Property 2 (V-shape) and the node is fathomed.

For each generated node, a lower bound on the weighted sum of earliness and tardiness penalties is computed. A node is fathomed if this estimate is greater than or equal to the current incumbent solution.

The proposed lower bound (LB2) considers the penalties caused by the known jobs in the partial sequence. It also considers the estimates of the penalty generated by the job that will finish its processing at the starting time of the partial sequence and by the subsequence of unknown jobs that start its processing at time zero. The following lemma is needed to establish the lower bound.

Lemma 1. Let A = a1, a2,..., at and B = b1, b2,..., bt sequences of positive integer numbers and define Z(A, B) = . If a1> a2> ... > at and b1< b2< ... < bt then Z(A, B) is minimized.

Proof. Assume that Z is minimized by sequences A and B that do not satisfy the hypothesis. We will show that there is another pair of sequences A' and B' satisfying the hypothesis that also minimizes Z. As A and B do not satisfy the hypothesis then there exists at least one index v such that av < av+1 or bv > bv+1. We are left to three different cases:

i) av < av+1 and bv < bv+1,

ii) av < av+1 and bv > bv+1,

iii) av > av+1 and bv > bv+1.

Interchanging av and av+1 in case i), interchanging av and av+1 and bv and bv+1 in case ii), or interchanging bv and bv+1 in case iii), it is easy to see that the associated value of Z evaluated at the modified sequences does not increase (in fact, it remains unchanged in case ii)). By continuing this pairwise interchange we can obtain ordered sequences A' and B' satisfying the hypothesis that minimize Z, as we wanted to prove.

Proposition 2. Consider n jobs and a known subsequence R of these jobs that starts its processing at

> d and finishes at
. Let D be the set of the remaining jobs. The weighted sum of earliness and tardiness penalties of the complete schedule is greater than or equal to:

where |D| is the cardinality of D and, among jobs in D,

is the q-th smallest earliness penalty and
is the b-th largest processing time.

Proof. Consider a known sequence W with n jobs that starts at time zero. The weighted sum of earliness and tardiness penalties of this sequence is given by:

Let R and D be defined as in the statement of this proposition and let π(y) be the job that is allocated in position y and u the job that finishes its processing at instant S. Then f(W) can be expressed as:

As all penalties are non-negative and βu> mink D {βk}, the following relations can be established:

Since implies that

it follows that:

We know that:

Using Lemma 1 (with a = (d - )+ and b = , for ℓ = 1,..., |D|) we have that:

Since LB2 is smaller than or equal to the weighted sum of earliness and tardiness penalties of an arbitrary sequence defined as in the statement of the proposition we conclude that the thesis holds.

4 Computational experiments

The proposed b&b algorithm was applied to two different sets of problems. In the first experiment, the instances were generated according to [3], so that they do not depend on the computer used. The authors of [3] were the pioneers in considering solutions starting at time instants different from zero. Seven different numbers of jobs n ∈ { 5, 10, 15, 20, 25, 30, 35 } and four restrictive factors h ∈ { 0.2, 0.4, 0.6, 0.8} were considered. The factor h indicates how jammed the production line is at the beginning of the schedule and it is used in the definition of the common due date, according to the expression: d = .

The processing times are integers uniformly distributed in [1, 20], the earliness penalties in [1,10] and the tardiness penalties in [1, 15]. There are 10 instances to be tested for each problem size and each restrictive factor, totaling 7 ×4 ×10 = 280 problems. The computer code was written using C language and the experiments were run on an Intel Core 2 Duo with a 2.40 GHz processor and 2.0 Gb of RAM memory. To prevent excessive computation time, the algorithm was stopped after 1 hour of CPU time for each problem. All instances are available at [25].

Table 1 presents the results when the b&b algorithm was applied to this first set of problems. The columns Min, Max and Average represent, respectively, the minimum, maximum and average CPU time obtained in each combination of number of jobs and restrictive factor. The column Number of solved problems indicates the quantity of instances solved within the time limit. As it can be seen, all the problems with n < 25 jobs were solved in an average time of less than 23 seconds.

In order to evaluate the proposed b&b we solve the same problem set using the MILP model described in Section 2 with the solver CPLEX 11.0 with its default parameter values. The computation time for each test instance was also limited to 1 hour. Analyzing Table 1 it can be noted that the CPLEX solver was not able to prove optimality within the allowed execution time in all instances with more than 10 jobs while the b&b algorithm found the optimal solution for all instances with up to 30 jobs. Moreover, it can be observed that in all instances solved by both methods, the CPU time of the b&b was smaller than the one presented by the CPLEX solver. These results indicate the efficiency of the proposed algorithm in reducing the search space. This good performance was expected since the lower bounds and pruning rules proposed in this work are specific for the single machine earliness and tardiness scheduling problem. These essential components of the presented b&b are able to exploit properties of the problem, as, for example, the fact that if a sequence of jobs before or after the common due date is determined the complete schedule is already defined.

Aiming to contribute to the study of this problem, in a second experiment, the b&b algorithm was applied to instances with 10 and 20 jobs from the benchmark problems presented in [3]. These benchmark values are frequently used in literature (see, for example [11] and [18]). For the 10-jobs instances, the benchmark provides optimal solutions as well as the b&b algorithm. Considering problems with 20 jobs, the proposed algorithm obtained optimal solutions within the time limit for all problems. Table 2 shows the solution values for the 40 instances generated by the b&b algorithm.

We also computed the percentage difference of the optimal value (Fo) in relation to the benchmark value of Biskup and Feldman (FBF), as follows:

Table 2 shows the %Diff for each considered problem. It can be observed that, in the worst case scenario, the benchmark value is 11.37% greater than the optimal value. These optimal solutions are interesting for future comparison purposes.

5 Final remarks

This paper considered the single machine scheduling problem with restrictive common due date involving tardiness and earliness penalties. This type of problem became more important with the advent of the lean production principles, including the just-in-time (JIT) concept. Due to its complexity, most of the authors addressed this problem using heuristic and metaheuristic approaches.

In this study, a branch-and-bound algorithm was proposed to find optimal solutions to this problem. In the development of the algorithm, the use of problem properties was important for the development of new lower bounds and pruning rules that have enhanced the efficiency of the proposed method.

An implementation of the method was tested in 280 problems generated as presented in [3]. The proposed b&b outperformed the CPLEX optimization software. This software was unable to prove optimality in all instances with n > 15 (200 instances) within the time limit, while the b&b algorithm found the optimal solution in 92% of the considered instances. These results indicate the efficiency of the algorithm, mainly due to the elimination of inferior qualitysolutions through the use of the proposed lower bounds that exploit characteristics of the considered problem.

In addition, the optimal solution values obtained for the benchmark problems suggested in [3] revealed that these reference values can be improved by up to 11.37%. These results can be used to evaluate the performance of heuristics and meta-heuristics developed for this problem.

As an extension of this study, we suggest the use of properties (see [24]) in the development of a lower bound for more general cases, such as in the flowshop problem with multiple machines.

Acknowledgements. The authors would like to thank the anonymous referees whose comments helped a lot to improve this paper.

Received: 18/VI/09.

Accepted: 26/XI/09.

# CAM-107/09.

  • [1] M. van den Akker, H. Hoogeveen and S. van de Velde, Combining column generation and Lagrangean relaxation to solve a single-machine common due date. INFORMS Journal on Computing, 14 (2002), 37-51.
  • [2] K.R. Baker and G.D. Scudder, Sequencing with earliness and tardiness penalties: A review. Operations Research, 38 (1990), 22-36.
  • [3] D. Biskup and M. Feldmann, Benchmarks for scheduling on a single-machine against restrictive and unrestrictive common due dates. Computers and Operations Research, 28 (2001), 787-801.
  • [4] T.C.E. Cheng and H.G. Kahlbacher, A proof for the longest-job-first policy in one-machine scheduling. Naval Research Logistics, 38 (1991), 715-720.
  • [5] P. De, J.B. Ghosh and C.E. Wells, Solving a generalized model for CON due date assignment and sequencing. International Journal of Production Economics, 34 (1994), 179-185.
  • [6] P. Dileepan, Common due date scheduling problem with separate earliness and tardiness penalties. Computers and Operations Research, 20 (1993), 179-184.
  • [7] M. Feldmann and D. Biskup, Single-machine scheduling for minimizing earliness and tardiness penalties by meta-heuristic approaches. Computers and Industrial Engineering, 44 (2003), 307-323.
  • [8] V. Gordon, J.M. Proth and C. Chu, A survey of the state-of-art of common due date assignment and scheduling research. European Journal of Operational Research, 139 (2002), 1-25.
  • [9] N.G. Hall and M.E. Posner, Earliness-tardiness scheduling problems, I: weighted deviation of completion times about a common due date. Operations Research, 39 (1991), 836-846.
  • [10] N. Hall, G.W. Kubiak and S.P. Sethi, Earliness-tardiness scheduling problems, II: deviations of completion times about a restrictive common due date. Operations Research, 39 (1991), 847-856.
  • [11] C.M. Hino, D.P. Ronconi and A.B. Mendes, Minimizing earliness and tardiness penalties in a single-machine problem with a common due date. European Journal of Operational Research, 160 (2005), 190-201.
  • [12] J.A. Hoogeven and S.L. van de Velde, Scheduling around a small common due date. European Journal of Operational Research, 55 (1991), 237-242.
  • [13] R.J.W. James, Using tabu search to solve the common due date early/tardy machine scheduling problem. Computers and Operations Research, 24 (1997), 199-208.
  • [14] H.G. Kahlbacher, Scheduling with monotonous earliness and tardiness penalties. European Journal of Operational Research, 64 (1993), 258-277.
  • [15] J.J. Kanet, Minimizing the average deviation of job completion times about a common due date. Naval Research Logistics Quartely, 28 (1981), 643-651.
  • [16] M.S. Kawamura and D.P. Ronconi, Aplicação do método branch-and-bound na programação de tarefas em uma única máquina com data de entrega comum sob penalidades de adiantamento e atraso. Technical Report 0606, Escola Politécnica, Universidade de São Paulo (2006).
  • [17] C.Y. Lee and S.J. Kim, Parallel genetic algorithm for the earliness-tardiness job scheduling problem with general penalty weights. Computers and Industrial Engineering, 28 (1995), 231-243.
  • [18] C-J. Liao and C-C. Cheng, A variable neighborhood search for minimizing single machine weighted earliness and tardiness with common due date. Computers and Industrial Engineering, 52 (2007), 404-413.
  • [19] C.-F. Liaw, A branch-and-bound algorithm for the single machine earliness and tardiness scheduling problem. Computers and Operations Research, 26 (1999), 679-693.
  • [20] S.A. Mondal and A.K. Sen, Single machine weighted earliness-tardiness penalty problem with a common due date. Computers and Operations Research, 28 (2001), 649-669.
  • [21] S.S. Panwalkar, M.L. Smith and A. Seidmann, Common due date assignment to minimize total penalty for one machine problem. Operations Research, 30 (1982), 391-399.
  • [22] D.P. Ronconi, A Branch-and-Bound Algorithm to Minimize the Makespan in a Flowshop with Blocking. Annals of Operations Research, 138 (2005), 53-65.
  • [23] D.P. Ronconi and V.A. Armentano, Lower Bounding Schemes for Flowshops with Blocking In-Process. Journal of the Operational Research Society, 52 (2001), 1289-1297.
  • [24] C.S. Sakuraba, D.P. Ronconi and F. Sourd, Scheduling in a two-machine flowshop for the minimization of the mean absolute deviation from a common due date. Computers and Operations Research, 36 (2009), 60-72.
  • [25] http://www.poli.usp.br/pro/docentes/dronconi/
    » link
  • *
    This research has been financially supported by "Funda ção de Amparo à Pesquisa do Estado de São Paulo" - FAPESP (Grant 06/03496-3) and by "Conselho Nacional de Desenvolvimento Científico e Tecnológico" - CNPq (Grants 486124/2007-0 and 307399/2006-0).
  • Publication Dates

    • Publication in this collection
      22 July 2010
    • Date of issue
      June 2010

    History

    • Received
      18 June 2009
    • Accepted
      26 June 2009
    Sociedade Brasileira de Matemática Aplicada e Computacional Sociedade Brasileira de Matemática Aplicada e Computacional - SBMAC, Rua Maestro João Seppe, nº. 900 , 16º. andar - Sala 163, 13561-120 São Carlos - SP Brasil, Tel./Fax: 55 16 3412-9752 - São Carlos - SP - Brazil
    E-mail: sbmac@sbmac.org.br