In this paper, the stability of a cubic functional equation in the setting of intuitionistic random normed spaces is proved. We first introduce the notation of intuitionistic random normed spaces. Then, by virtue of t...In this paper, the stability of a cubic functional equation in the setting of intuitionistic random normed spaces is proved. We first introduce the notation of intuitionistic random normed spaces. Then, by virtue of this notation, we study the stability of a cubic functional equation in the setting of these spaces under arbitrary triangle norms. Furthermore, we present the interdisciplinary relation among the theory of random spaces, the theory of intuitionistic spaces, and the theory of functional equations.展开更多
Let G be an Abelian group and letρ：G×G→[0,∞） be a metric on G. Let E be a normed space. We prove that under some conditions if f：G→E is an odd function and Cx：G→E defined by Cx（y）：=2 f （x＋y）＋2 f ...Let G be an Abelian group and letρ：G×G→[0,∞） be a metric on G. Let E be a normed space. We prove that under some conditions if f：G→E is an odd function and Cx：G→E defined by Cx（y）：=2 f （x＋y）＋2 f （x-y）＋12 f （x）-f （2x＋y）-f （2x-y） is a cubic function for all x∈G, then there exists a cubic function C：G→E such that f？C is Lipschitz. Moreover, we investigate the stability of cubic functional equation 2 f （x＋y）＋2 f （x-y）＋12 f （x）-f （2x＋y）-f （2x-y）=0 on Lipschitz spaces.展开更多
In this paper, we determine the general solution of the functional equation f1 （2x ＋ y） ＋ f2（2x - y） = f3（x ＋ y） ＋ f4（x - y） ＋ f5（x） without assuming any regularity condition on the unknown functions f1...In this paper, we determine the general solution of the functional equation f1 （2x ＋ y） ＋ f2（2x - y） = f3（x ＋ y） ＋ f4（x - y） ＋ f5（x） without assuming any regularity condition on the unknown functions f1,f2,f3, f4, f5 ： R→R. The general solution of this equation is obtained by finding the general solution of the functional equations f（2x ＋ y） ＋ f（2x - y） = g（x ＋ y） ＋ g（x - y） ＋ h（x） and f（2x ＋ y） - f（2x - y） = g（x ＋ y） - g（x - y）. The method used for solving these functional equations is elementary but exploits an important result due to Hosszfi. The solution of this functional equation can also be determined in certain type of groups using two important results due to Szekelyhidi.展开更多
In this paper, the direct method and the fixed point alternative method are implemented to give Hyers-Uiam-Rassias stability of the functional equation 6f（x＋y）-6f（x-y）＋4f（3y）=3f（x＋2y）-3f（x-2y）＋9f（2y） i...In this paper, the direct method and the fixed point alternative method are implemented to give Hyers-Uiam-Rassias stability of the functional equation 6f（x＋y）-6f（x-y）＋4f（3y）=3f（x＋2y）-3f（x-2y）＋9f（2y） in fuzzy Banach spaces. We can find the range of approximate solutions obtained using the direct method are less than those obtained by using the fixed point alternative method for the above and the functional equation.展开更多
In this paper, we will find out the general solution and investigate the generalized Hyers-Ulam-Rassias stability problem for the following cubic functional equation2f（x ＋ 2y） ＋ f（2x - y） = 5f（x ＋ y） ＋ 5f（x...In this paper, we will find out the general solution and investigate the generalized Hyers-Ulam-Rassias stability problem for the following cubic functional equation2f（x ＋ 2y） ＋ f（2x - y） = 5f（x ＋ y） ＋ 5f（x - y）＋ 15f（y）in the spirit of Hyers, Ulam, Rassias and Gavruta.展开更多
In this paper,we obtain the general solution and stability of the Jensen-cubic functional equation f（（x1＋x2）/2,2y1＋y2）＋f（（x1＋x2）/2,2（y1-y2）） = f（x1,y1 ＋y2）＋f（x1,y1-y2）＋6f（x1,y1＋ f（x2,y1y2）＋f...In this paper,we obtain the general solution and stability of the Jensen-cubic functional equation f（（x1＋x2）/2,2y1＋y2）＋f（（x1＋x2）/2,2（y1-y2）） = f（x1,y1 ＋y2）＋f（x1,y1-y2）＋6f（x1,y1＋ f（x2,y1y2）＋f（x2,y1-y2）＋6f（x2,y1）.展开更多
In this paper, we investigate the general solution and the stability of a cubic functional equation f（x ＋ ny） ＋ f（x - ny） ＋ f（nx） = n^2 f（x ＋ y） ＋ n^2 f（x - y）＋ （n^3 - 2n^2 ＋ 2）f（x）,where n ≥ 2 i...In this paper, we investigate the general solution and the stability of a cubic functional equation f（x ＋ ny） ＋ f（x - ny） ＋ f（nx） = n^2 f（x ＋ y） ＋ n^2 f（x - y）＋ （n^3 - 2n^2 ＋ 2）f（x）,where n ≥ 2 is an integer. Furthermore, we prove the stability by the fixed point method.展开更多
Reliability allocation of computerized numerical controlled（CNC）lathes is very important in industry.Traditional allocation methods only focus on high-failure rate components rather than moderate failure rate compon...Reliability allocation of computerized numerical controlled（CNC）lathes is very important in industry.Traditional allocation methods only focus on high-failure rate components rather than moderate failure rate components,which is not applicable in some conditions.Aiming at solving the problem of CNC lathes reliability allocating,a comprehensive reliability allocation method based on cubic transformed functions of failure modes and effects analysis（FMEA）is presented.Firstly,conventional reliability allocation methods are introduced.Then the limitations of direct combination of comprehensive allocation method with the exponential transformed FMEA method are investigated.Subsequently,a cubic transformed function is established in order to overcome these limitations.Properties of the new transformed functions are discussed by considering the failure severity and the failure occurrence.Designers can choose appropriate transform amplitudes according to their requirements.Finally,a CNC lathe and a spindle system are used as an example to verify the new allocation method.Seven criteria are considered to compare the results of the new method with traditional methods.The allocation results indicate that the new method is more flexible than traditional methods.By employing the new cubic transformed function,the method covers a wider range of problems in CNC reliability allocation without losing the advantages of traditional methods.展开更多
Modelling and simulation of projectile flight is at the core of ballistic computer software and is essential to the study of performance of rifles and projectiles in various engagement conditions.An effective and repr...Modelling and simulation of projectile flight is at the core of ballistic computer software and is essential to the study of performance of rifles and projectiles in various engagement conditions.An effective and representative numerical model of projectile flight requires a relatively good approximation of the aerodynamics.The aerodynamic coefficients of the projectile model should be described as a series of piecewise polynomial functions of the Mach number that ideally meet the following conditions:they are continuous,differentiable at least once,and have a relatively low degree.The paper provides the steps needed to generate such piecewise polynomial functions using readily available tools,and then compares Piecewise Cubic Hermite Interpolating Polynomial(PCHIP),cubic splines,and piecewise linear functions,and their variant,as potential curve fitting methods to approximate the aerodynamics of a generic small arms projectile.A key contribution of the paper is the application of PCHIP to the approximation of projectile aerodynamics,and its evaluation against a set of criteria.Finally,the paper provides a baseline assessment of the impact of the polynomial functions on flight trajectory predictions obtained with 6-degree-of-freedom simulations of a generic projectile.展开更多
Hermite interpolation is a very important tool in approximation theory and nu- merical analysis, and provides a popular method for modeling in the area of computer aided geometric design. However, the classical Hermit...Hermite interpolation is a very important tool in approximation theory and nu- merical analysis, and provides a popular method for modeling in the area of computer aided geometric design. However, the classical Hermite interpolant is unique for a prescribed data set, and hence lacks freedom for the choice of an interpolating curve, which is a crucial requirement in design environment. Even though there is a rather well developed fractal theory for Hermite interpolation that offers a large flexibility in the choice of interpolants, it also has the short- coming that the functions that can be well approximated are highly restricted to the class of self-affine functions. The primary objective of this paper is to suggest a gl-cubic Hermite in- terpolation scheme using a fractal methodology, namely, the coalescence hidden variable fractal interpolation, which works equally well for the approximation of a self-affine and non-self-affine data generating functions. The uniform error bound for the proposed fractal interpolant is established to demonstrate that the convergence properties are similar to that of the classical Hermite interpolant. For the Hermite interpolation problem, if the derivative values are not actually prescribed at the knots, then we assign these values so that the interpolant gains global G2-continuity. Consequently, the procedure culminates with the construction of cubic spline coalescence hidden variable fractal interpolants. Thus, the present article also provides an al- ternative to the construction of cubic spline coalescence hidden variable fractal interpolation functions through moments proposed by Chand and Kapoor [Fractals, 15（1） （2007）, pp. 41-53].展开更多
This paper deals with a Unit Commitment (UC) problem of a power plant aimed to find the optimal scheduling of the generating units involving cubic cost functions. The problem has non convex generator characteristics, ...This paper deals with a Unit Commitment (UC) problem of a power plant aimed to find the optimal scheduling of the generating units involving cubic cost functions. The problem has non convex generator characteristics, which makes it very hard to handle the corresponding mathematical models. However, Teaching Learning Based Optimization (TLBO) has reached a high efficiency, in terms of solution accuracy and computing time for such non convex problems. Hence, TLBO is applied for scheduling of generators with higher order cost characteristics, and turns out to be computationally solvable. In particular, we represent a model that takes into account the accurate higher order generator cost functions along with ramp limits, and turns to be more general and efficient than those available in the literature. The behavior of the model is analyzed through proposed technique on modified IEEE-24 bus system.展开更多
In this paper,the kernel of the cubic spline interpolation is given.An optimal error bound for the cu- bic spline interpolation of lower smooth functions is obtained.
In this paper,the forecasting equations of a 2nd-order space-time differential remainder are deduced from the Navier-Stokes primitive equations and Eulerian operator by Taylor-series expansion.Here we introduce a cubi...In this paper,the forecasting equations of a 2nd-order space-time differential remainder are deduced from the Navier-Stokes primitive equations and Eulerian operator by Taylor-series expansion.Here we introduce a cubic spline numerical model(Spline Model for short),which is with a quasi-Lagrangian time-split integration scheme of fitting cubic spline/bicubic surface to all physical variable fields in the atmospheric equations on spherical discrete latitude-longitude mesh.A new algorithm of"fitting cubic spline—time step integration—fitting cubic spline—……"is developed to determine their first-and2nd-order derivatives and their upstream points for time discrete integral to the governing equations in Spline Model.And the cubic spline function and its mathematical polarities are also discussed to understand the Spline Model’s mathematical foundation of numerical analysis.It is pointed out that the Spline Model has mathematical laws of"convergence"of the cubic spline functions contracting to the original functions as well as its 1st-order and 2nd-order derivatives.The"optimality"of the 2nd-order derivative of the cubic spline functions is optimal approximation to that of the original functions.In addition,a Hermite bicubic patch is equivalent to operate on a grid for a 2nd-order derivative variable field.Besides,the slopes and curvatures of a central difference are identified respectively,with a smoothing coefficient of 1/3,three-point smoothing of that of a cubic spline.Then the slopes and curvatures of a central difference are calculated from the smoothing coefficient 1/3 and three-point smoothing of that of a cubic spline,respectively.Furthermore,a global simulation case of adiabatic,non-frictional and"incompressible"model atmosphere is shown with the quasi-Lagrangian time integration by using a global Spline Model,whose initial condition comes from the NCEP reanalysis data,along with quasi-uniform latitude-longitude grids and the so-called"shallow atmosphere"Navier-Stokes primitive equations in the spherical coordinates.The Spline Model,which adopted the Navier-Stokes primitive equations and quasi-Lagrangian time-split integration scheme,provides an initial ideal case of global atmospheric circulation.In addition,considering the essentially non-linear atmospheric motions,the Spline Model could judge reasonably well simple points of any smoothed variable field according to its fitting spline curvatures that must conform to its physical interpretation.展开更多
A closed but approximate formula of Green’s function for an arbitrary aggregate of cubic crystallites is given to derive the e?ective elastic sti?ness tensor of the polycrystal. This formula, which includes thr...A closed but approximate formula of Green’s function for an arbitrary aggregate of cubic crystallites is given to derive the e?ective elastic sti?ness tensor of the polycrystal. This formula, which includes three elastic constants of single cubic crystal and ?ve texture coe?cients, accounts for the e?ects of the orientation distribution function (ODF) up to terms linear in the tex- ture coe?cients. Thus it is expected that our formula would be applicable to arbitrary aggregates with weak texture or to materials such as aluminum whose single crystal has weak anisotropy. Three examples are presented to compare predictions from our formula with those from Nishioka and Lothe’s formula and Synge’s contour integral through numerical integration. As an applica- tion of Green’s function, we brie?y describe the procedure of deriving the e?ective elastic sti?ness tensor for an orthorhombic aggregate of cubic crystallites. The comparison of the computational results given by the ?nite element method and our e?ective elastic sti?ness tensor is made by an example.展开更多
The equation of state(EOS) for hard-sphere fluid derived from compressibility routes of Percus-Yevick theory(PYC) is extended. The two parameters are determined by fitting well-known virial coefficients of pure fluid....The equation of state(EOS) for hard-sphere fluid derived from compressibility routes of Percus-Yevick theory(PYC) is extended. The two parameters are determined by fitting well-known virial coefficients of pure fluid.The extended cubic EOS can be directly extended to multi-component mixtures, merely demanding the EOS of mixtures also is cubic and combining two physical conditions for the radial distribution functions at contact(RDFC) of mixtures.The calculated virial coefficients of pure fluid and predicted compressibility factors and RDFC for both pure fluid and mixtures are excellent as compared with the simulation data. The values of RDFC for mixtures with extremely large size ratio 10 are far better than the BGHLL expressions in literature.展开更多
基金supported by the Natural Science Foundation of Yibin University (No. 2009Z003)
文摘In this paper, the stability of a cubic functional equation in the setting of intuitionistic random normed spaces is proved. We first introduce the notation of intuitionistic random normed spaces. Then, by virtue of this notation, we study the stability of a cubic functional equation in the setting of these spaces under arbitrary triangle norms. Furthermore, we present the interdisciplinary relation among the theory of random spaces, the theory of intuitionistic spaces, and the theory of functional equations.
文摘Let G be an Abelian group and letρ：G×G→[0,∞） be a metric on G. Let E be a normed space. We prove that under some conditions if f：G→E is an odd function and Cx：G→E defined by Cx（y）：=2 f （x＋y）＋2 f （x-y）＋12 f （x）-f （2x＋y）-f （2x-y） is a cubic function for all x∈G, then there exists a cubic function C：G→E such that f？C is Lipschitz. Moreover, we investigate the stability of cubic functional equation 2 f （x＋y）＋2 f （x-y）＋12 f （x）-f （2x＋y）-f （2x-y）=0 on Lipschitz spaces.
文摘In this paper, we determine the general solution of the functional equation f1 （2x ＋ y） ＋ f2（2x - y） = f3（x ＋ y） ＋ f4（x - y） ＋ f5（x） without assuming any regularity condition on the unknown functions f1,f2,f3, f4, f5 ： R→R. The general solution of this equation is obtained by finding the general solution of the functional equations f（2x ＋ y） ＋ f（2x - y） = g（x ＋ y） ＋ g（x - y） ＋ h（x） and f（2x ＋ y） - f（2x - y） = g（x ＋ y） - g（x - y）. The method used for solving these functional equations is elementary but exploits an important result due to Hosszfi. The solution of this functional equation can also be determined in certain type of groups using two important results due to Szekelyhidi.
文摘In this paper, the direct method and the fixed point alternative method are implemented to give Hyers-Uiam-Rassias stability of the functional equation 6f（x＋y）-6f（x-y）＋4f（3y）=3f（x＋2y）-3f（x-2y）＋9f（2y） in fuzzy Banach spaces. We can find the range of approximate solutions obtained using the direct method are less than those obtained by using the fixed point alternative method for the above and the functional equation.
基金Korea Research Foundation Grant KRF-2007-313-C00033
文摘In this paper, we will find out the general solution and investigate the generalized Hyers-Ulam-Rassias stability problem for the following cubic functional equation2f（x ＋ 2y） ＋ f（2x - y） = 5f（x ＋ y） ＋ 5f（x - y）＋ 15f（y）in the spirit of Hyers, Ulam, Rassias and Gavruta.
文摘In this paper,we obtain the general solution and stability of the Jensen-cubic functional equation f（（x1＋x2）/2,2y1＋y2）＋f（（x1＋x2）/2,2（y1-y2）） = f（x1,y1 ＋y2）＋f（x1,y1-y2）＋6f（x1,y1＋ f（x2,y1y2）＋f（x2,y1-y2）＋6f（x2,y1）.
文摘In this paper, we investigate the general solution and the stability of a cubic functional equation f（x ＋ ny） ＋ f（x - ny） ＋ f（nx） = n^2 f（x ＋ y） ＋ n^2 f（x - y）＋ （n^3 - 2n^2 ＋ 2）f（x）,where n ≥ 2 is an integer. Furthermore, we prove the stability by the fixed point method.
基金Supported by National Natural Science Foundation of China(Grant Nos.51135003,51205050,U1234208)Key National Science & Technology Special Project on"High-Grade CNC Machine Tools and Basic Manufacturing Equipments"(Grant No.2013ZX04011011)+1 种基金Research Fund for the Doctoral Program of Higher Education of China(Grant No.20110042120020)Fundamental Research Funds for the Central
文摘Reliability allocation of computerized numerical controlled（CNC）lathes is very important in industry.Traditional allocation methods only focus on high-failure rate components rather than moderate failure rate components,which is not applicable in some conditions.Aiming at solving the problem of CNC lathes reliability allocating,a comprehensive reliability allocation method based on cubic transformed functions of failure modes and effects analysis（FMEA）is presented.Firstly,conventional reliability allocation methods are introduced.Then the limitations of direct combination of comprehensive allocation method with the exponential transformed FMEA method are investigated.Subsequently,a cubic transformed function is established in order to overcome these limitations.Properties of the new transformed functions are discussed by considering the failure severity and the failure occurrence.Designers can choose appropriate transform amplitudes according to their requirements.Finally,a CNC lathe and a spindle system are used as an example to verify the new allocation method.Seven criteria are considered to compare the results of the new method with traditional methods.The allocation results indicate that the new method is more flexible than traditional methods.By employing the new cubic transformed function,the method covers a wider range of problems in CNC reliability allocation without losing the advantages of traditional methods.
文摘Modelling and simulation of projectile flight is at the core of ballistic computer software and is essential to the study of performance of rifles and projectiles in various engagement conditions.An effective and representative numerical model of projectile flight requires a relatively good approximation of the aerodynamics.The aerodynamic coefficients of the projectile model should be described as a series of piecewise polynomial functions of the Mach number that ideally meet the following conditions:they are continuous,differentiable at least once,and have a relatively low degree.The paper provides the steps needed to generate such piecewise polynomial functions using readily available tools,and then compares Piecewise Cubic Hermite Interpolating Polynomial(PCHIP),cubic splines,and piecewise linear functions,and their variant,as potential curve fitting methods to approximate the aerodynamics of a generic small arms projectile.A key contribution of the paper is the application of PCHIP to the approximation of projectile aerodynamics,and its evaluation against a set of criteria.Finally,the paper provides a baseline assessment of the impact of the polynomial functions on flight trajectory predictions obtained with 6-degree-of-freedom simulations of a generic projectile.
基金partially supported by the CSIR India(Grant No.09/084(0531)/2010-EMR-I)the SERC,DST India(Project No.SR/S4/MS:694/10)
文摘Hermite interpolation is a very important tool in approximation theory and nu- merical analysis, and provides a popular method for modeling in the area of computer aided geometric design. However, the classical Hermite interpolant is unique for a prescribed data set, and hence lacks freedom for the choice of an interpolating curve, which is a crucial requirement in design environment. Even though there is a rather well developed fractal theory for Hermite interpolation that offers a large flexibility in the choice of interpolants, it also has the short- coming that the functions that can be well approximated are highly restricted to the class of self-affine functions. The primary objective of this paper is to suggest a gl-cubic Hermite in- terpolation scheme using a fractal methodology, namely, the coalescence hidden variable fractal interpolation, which works equally well for the approximation of a self-affine and non-self-affine data generating functions. The uniform error bound for the proposed fractal interpolant is established to demonstrate that the convergence properties are similar to that of the classical Hermite interpolant. For the Hermite interpolation problem, if the derivative values are not actually prescribed at the knots, then we assign these values so that the interpolant gains global G2-continuity. Consequently, the procedure culminates with the construction of cubic spline coalescence hidden variable fractal interpolants. Thus, the present article also provides an al- ternative to the construction of cubic spline coalescence hidden variable fractal interpolation functions through moments proposed by Chand and Kapoor [Fractals, 15（1） （2007）, pp. 41-53].
文摘This paper deals with a Unit Commitment (UC) problem of a power plant aimed to find the optimal scheduling of the generating units involving cubic cost functions. The problem has non convex generator characteristics, which makes it very hard to handle the corresponding mathematical models. However, Teaching Learning Based Optimization (TLBO) has reached a high efficiency, in terms of solution accuracy and computing time for such non convex problems. Hence, TLBO is applied for scheduling of generators with higher order cost characteristics, and turns out to be computationally solvable. In particular, we represent a model that takes into account the accurate higher order generator cost functions along with ramp limits, and turns to be more general and efficient than those available in the literature. The behavior of the model is analyzed through proposed technique on modified IEEE-24 bus system.
文摘In this paper,the kernel of the cubic spline interpolation is given.An optimal error bound for the cu- bic spline interpolation of lower smooth functions is obtained.
文摘In this paper,the forecasting equations of a 2nd-order space-time differential remainder are deduced from the Navier-Stokes primitive equations and Eulerian operator by Taylor-series expansion.Here we introduce a cubic spline numerical model(Spline Model for short),which is with a quasi-Lagrangian time-split integration scheme of fitting cubic spline/bicubic surface to all physical variable fields in the atmospheric equations on spherical discrete latitude-longitude mesh.A new algorithm of"fitting cubic spline—time step integration—fitting cubic spline—……"is developed to determine their first-and2nd-order derivatives and their upstream points for time discrete integral to the governing equations in Spline Model.And the cubic spline function and its mathematical polarities are also discussed to understand the Spline Model’s mathematical foundation of numerical analysis.It is pointed out that the Spline Model has mathematical laws of"convergence"of the cubic spline functions contracting to the original functions as well as its 1st-order and 2nd-order derivatives.The"optimality"of the 2nd-order derivative of the cubic spline functions is optimal approximation to that of the original functions.In addition,a Hermite bicubic patch is equivalent to operate on a grid for a 2nd-order derivative variable field.Besides,the slopes and curvatures of a central difference are identified respectively,with a smoothing coefficient of 1/3,three-point smoothing of that of a cubic spline.Then the slopes and curvatures of a central difference are calculated from the smoothing coefficient 1/3 and three-point smoothing of that of a cubic spline,respectively.Furthermore,a global simulation case of adiabatic,non-frictional and"incompressible"model atmosphere is shown with the quasi-Lagrangian time integration by using a global Spline Model,whose initial condition comes from the NCEP reanalysis data,along with quasi-uniform latitude-longitude grids and the so-called"shallow atmosphere"Navier-Stokes primitive equations in the spherical coordinates.The Spline Model,which adopted the Navier-Stokes primitive equations and quasi-Lagrangian time-split integration scheme,provides an initial ideal case of global atmospheric circulation.In addition,considering the essentially non-linear atmospheric motions,the Spline Model could judge reasonably well simple points of any smoothed variable field according to its fitting spline curvatures that must conform to its physical interpretation.
基金Project supported by the Natural Science Foundation of Jiangxi Province (No. 0450035).
文摘A closed but approximate formula of Green’s function for an arbitrary aggregate of cubic crystallites is given to derive the e?ective elastic sti?ness tensor of the polycrystal. This formula, which includes three elastic constants of single cubic crystal and ?ve texture coe?cients, accounts for the e?ects of the orientation distribution function (ODF) up to terms linear in the tex- ture coe?cients. Thus it is expected that our formula would be applicable to arbitrary aggregates with weak texture or to materials such as aluminum whose single crystal has weak anisotropy. Three examples are presented to compare predictions from our formula with those from Nishioka and Lothe’s formula and Synge’s contour integral through numerical integration. As an applica- tion of Green’s function, we brie?y describe the procedure of deriving the e?ective elastic sti?ness tensor for an orthorhombic aggregate of cubic crystallites. The comparison of the computational results given by the ?nite element method and our e?ective elastic sti?ness tensor is made by an example.
基金Supported by the Science and Technology Foundation of State Key Laboratory for Shock Wave and Detonation Physics under Grant No.9140C670103120C6702the Program for Excellent Talents of Sichuan Province of China under Grant No.2011JQ0053University Electronic Science and Technology of China under Grant No.23601008
文摘The equation of state(EOS) for hard-sphere fluid derived from compressibility routes of Percus-Yevick theory(PYC) is extended. The two parameters are determined by fitting well-known virial coefficients of pure fluid.The extended cubic EOS can be directly extended to multi-component mixtures, merely demanding the EOS of mixtures also is cubic and combining two physical conditions for the radial distribution functions at contact(RDFC) of mixtures.The calculated virial coefficients of pure fluid and predicted compressibility factors and RDFC for both pure fluid and mixtures are excellent as compared with the simulation data. The values of RDFC for mixtures with extremely large size ratio 10 are far better than the BGHLL expressions in literature.
