A semismooth equation approach to the solution of nonlinear complementarity problems

In this paper we present a new algorithm for the solution of nonlinear complementarity problems. The algorithm is based on a semismooth equation reformulation of the complementarity problem. We exploit the recent extension of Newton's method to semismooth systems of equations and the fact that the natural merit function associated to the equation reformulation is continuously differentiable to develop an algorithm whose global and quadratic convergence properties can be established under very mild assumptions. Other interesting features of the new algorithm are an extreme simplicity along with a low computational burden per iteration. We include numerical tests which show the viability of the approach.     

Power series solution to the pendulum equation

This note gives a power series solution to the pendulum equation that enables to investigate the system in an analytical way only, i.e. to avoid numeric methods. A method of determining the number of the terms for getting a required relative error is presented that uses bigger and lesser geometric series. The solution is suitable for modelling the motion in spreadsheets up to the amplitude of 90°.     

Initial-boundary-value problems for the generalized multi-term time-fractional diffusion equation

In this paper, the initial-boundary-value problems for the generalized multi-term time-fractional diffusion equation over an open bounded domain G×(0,T), G∈Rn are considered. Based on an appropriate maximum principle that is formulated and proved in the paper, too, some a priory estimates for the solution and then its uniqueness are established. To show the existence of the solution, first a formal solution is constructed using the Fourier method of the separation of the variables. The time-dependent components of the solution are given in terms of the multinomial Mittag-Leffler function. Under certain conditions, the formal solution is shown to be a generalized solution of the initial-boundary-value problem for the generalized time-fractional multi-term diffusion equation that turns out to be a classical solution under some additional conditions. Another important consequence from the maximum principle is a continuously dependence of the solution on the problem data (initial and boundary conditions and a source function) that - together with the uniqueness and existence results - makes the problem under consideration to a well-posed problem in the Hadamard sense.     

Francesco Carlini: Kepler's equation and the asymptotic solution to singular differential equations

Carlini's career was mainly dedicated to astronomy, but he was also a particularly skilled mathematician. In this article we collect and analyse his mathematical contributions in detail. In particular, in his important Memoir of the year 1817 devoted to Kepler's equation he introduced an innovative idea to solve ordinary differential equations with singular perturbations by means of asymptotic expansions. In the same Memoir also appeared, five years before Laplace's contributions, what is usually called the Laplace limit constant. Furthermore, Carlini published other mathematical Memoirs anticipating, 70 years in advance, the importance of complex branches of the Lambert's special function.     

双边配给问题的Shapley解及其在博物馆通票问题中的应用

双边配给问题描述了现实生活中一类带有二部图结构的稀缺资源配置问题, 例如, 在自然灾害期间救援物资的配给; 电力和天然气等自然资源按需分配; 高校引进人才调配等。本文通过求解线性规划, 并从联盟边际贡献的角度出发定义了双边配给问题的一个Shapley解。之后, 通过合作对策模型和解的公理化方法说明新解的合理性。首先, 建立双边配给问题的合作对策模型, 论证了新解与双边配给合作对策的Shapley值一致; 其次, 证明了Shapley解是唯一满足优先一致性的有效配给方案。最后, 将Shapley解应用于博物馆通票问题的研究, 探讨了博物馆合作制定通票后所得单票和通票收益的分配方式。     

The penetration function and its application to microscale problems

The penetration function measures the effect of the boundary data on the energy of the solution of a second order linear elliptic PDE taken over an interior subdomain. Here the coefficients of the PDE are functions of position and often represent the material properties of non homogeneous media with microstructure. The penetration function is used to assess the accuracy of global-local approaches for recovering local solution features from coarse grained solutions such as those delivered by homogenization theory. AMS subject classification (2000)  65N15, 78M40     

双边配给问题的Shapley解及其在博物馆通票问题中的应用

双边配给问题描述了现实生活中一类带有二部图结构的稀缺资源配置问题, 例如, 在自然灾害期间救援物资的配给; 电力和天然气等自然资源按需分配; 高校引进人才调配等。本文通过求解线性规划, 并从联盟边际贡献的角度出发定义了双边配给问题的一个Shapley解。之后, 通过合作对策模型和解的公理化方法说明新解的合理性。首先, 建立双边配给问题的合作对策模型, 论证了新解与双边配给合作对策的Shapley值一致; 其次, 证明了Shapley解是唯一满足优先一致性的有效配给方案。最后, 将Shapley解应用于博物馆通票问题的研究, 探讨了博物馆合作制定通票后所得单票和通票收益的分配方式。     

Interactive bicriterion solution method and its application to critical path method problems

An interactive solution method is developed for bicriterion mathematical programming (BCMP) problems. The new method, called the dichotomous bicriterion mathematical programming (DBCMP) method, combines Tchebycheff theory and the existing paired comparison method (PCM). The DBCMP method is then compared with the PCM method based on critical path method problems with two conflicting objectives: minimizing the total crashing cost and minimizing the total project completion time. The extension of the DBCMP method to BCMP problems with multiple decision makers is also discussed.     

Approximate solutions to Dirichlet control problems for the wave equation in Sobolev classes and dual observation problems

Dual control and observation problems for the wave equation with variable coefficients subject to Dirichlet boundary conditions are solved by a variational method. This method was earlier proposed by the author for an approximate analysis of linear equations with nonuniform perturbations of the operator. Explicit bounds on the constant that are required to implement the method are obtained using the correct solvability property of the dual observation problem. Finite-dimensional approximations of the control and observation problems are obtained by the difference method preserving the duality relation. The convergence of approximate solutions is established in the norms of the corresponding dual spaces.     

A new approximation for the generalized fractional derivative and its application to generalized fractional diffusion equation

In this paper, we derive a fourth order approximation for the generalized fractional derivative that is characterized by a scale function z(t) and a weight function w(t) . Combining the new approximation with compact finite difference method, we develop a numerical scheme for a generalized fractional diffusion problem. The stability and convergence of the numerical scheme are proved by the energy method, and it is shown that the temporal and spatial convergence orders are both 4. Several numerical experiments are provided to illustrate the efficiency of our scheme.     

On the solution set of convex problems and its numerical application

The solution set of a convex problem is enough interested. In this paper we have a discussion to determine it. The conditions and proofs in this note essentially differs from the previous well-known works. We also have a new result determining the solution of minimum norm for quadratic problems.     

The planar radiosity equation and its numerical solution

This article gives properties of the planar radiosity equationand methods for its numerical solution. Regularity propertiesof the radiosity solution are examined, including both the effectsof corners and the effects of the visibility function. Theseare taken into account in the design of collocation methodswith piecewise polynomial approximating functions. Numericalexamples conclude the paper.     

Comparing proposals for the solution of data analysis problems in a knowledge-based system

In a knowledge-based system, which aims at supporting persons who are interested in the analysis of special data, the problem can arise that a whole set of proposals is generated in answer to a question of a user. Such proposals are based on appropriate interconnections between user wishes, available original data as well as derived data obtained by application of adequate methods, the methods mentioned, and data analysis objectives. We use graphical visualizations of proposals to outline how the system would cope with the underlying situation.In this paper, special attention is paid to the concept of knowledge-based comparisons of proposals when propagation of certainty factors is used for a-priori judgments of proposals generated (before suggested proposals are performed). After-wards, a-posteriori judgments of proposals considered (after solutions have been computed by application of selected proposals) can be based on goodness of fit criteria derived from chosen outputs.     

Applications to give an analytical solution to the Black Scholes equation

A simple ?-function and its dynamic equation is presented. An application to give an analytical solution to the Black Scholes equation is presented.     

The uniqueness of the solution to the diffusion equation with a damping term

Consider a non-linear diffusion equation with a damping term. If the diffusion coefficient is positive, then the solutions are not unique generally. However, if the diffusion coefficient degenerates, the situation may change. In this paper, not only the existence of the weak solution is established, but also the uniqueness of the weak solutions is proved, even the boundary value condition is not imposed. The conclusions imply that, on the boundary, the degeneracy of diffusion coefficient can eliminate the action from the damping term.     

A matroid algorithm and its application to the efficient solution of two optimization problems on graphs

We address the problem of finding a minimum weight baseB of a matroid when, in addition, each element of the matroid is colored with one ofm colors and there are upper and lower bound restrictions on the number of elements ofB with colori, fori = 1, 2,,m. This problem is a special case of matroid intersection. We present an algorithm that exploits the special structure, and we apply it to two optimization problems on graphs. When applied to the weighted bipartite matching problem, our algorithm has complexity O(|EV|+|V| ²log|V|). HereV denotes the node set of the underlying bipartite graph, andE denotes its edge set. The second application is defined on a general connected graphG = (V,E) whose edges have a weight and a color. One seeks a minimum weight spanning tree with upper and lower bound restrictions on the number of edges with colori in the tree, for eachi. Our algorithm for this problem has complexity O(|EV|+m ² |V|+ m|V| ²). A special case of this constrained spanning tree problem occurs whenV ^* is a set of pairwise nonadjacent nodes ofG. One must find a minimum weight spanning tree with upper and lower bound restrictions on the degree of each node ofV ^*. Then the complexity of our algorithm is O(|VE|+|V ^* V| ²). Finally, we discuss a new relaxation of the traveling salesman problem.This report was supported in part by NSF grant ECS 8601660.     

A improved F‐expansion method and its application to the Zhiber–Shabat equation

Based on the F‐expansion method and Exp‐function method, an improved F‐expansion method is introduced. As illustrative examples, the exact solutions expressed by exponential function, hyperbolic functions, logarithmic function, and other type of functions for the Zhiber–Shabat equation are derived. Some previous results are extended. Copyright © 2012 John Wiley & Sons, Ltd.     

A new approach to numerical solution of fixed-point problems and its application to delay differential equations

In this paper we consider a certain approximation of fixed-points of a continuous operator A mapping the metric space into itself by means of finite dimensional ε(h)-fixed-points of A. These finite dimensional functions are obtained from functions defined on discrete space grid points (related to a parameter h→0) by applying suitably chosen extension operators p_h. A theorem specifying necessary and sufficient conditions for existence of fixed-points of A in terms of ε(h)-fixed-points of A is given. A corollary which follows the theorem yields an approximate method for a fixed-point problem and determines conditions for its convergence. An example of application of the obtained general results to numerical solving of boundary value problems for delay differential equations is provided.Numerical experiments carried out on three examples of boundary value problems for second order delay differential equations show that the proposed approach produces much more accurate results than many other numerical methods when applied to the same examples.     

Ergodic BSDE with unbounded and multiplicative underlying diffusion and application to large time behaviour of viscosity solution of HJB equation

We study ergodic backward stochastic differential equations (EBSDEs), for which the underlying diffusion is assumed to be multiplicative and of linear growth. The fact that the forward process has an unbounded diffusion is balanced with an assumption of weak dissipativity for its drift. Moreover, the forward equation is assumed to be non-degenerate. We study the existence and uniqueness of EBSDEs and we apply our results to an ergodic optimal control problem. In particular, we show the large time behaviour of viscosity solution of Hamilton–Jacobi–Bellman equation with an exponential rate of convergence when the underlying diffusion is multiplicative and unbounded.