• Journal of applied mathematics & informatics
• /
• v.16 no.1_2
• /
• pp.91-104
• /
• 2004

A modification of the gradient method of convex programming is introduced. Also, we describe symbolic implementation of the gradient method and its modification by means of the programming language MATHEMATICA. A few numerical examples are reported.

• The Journal of Information Systems
• /
• v.5
• /
• pp.21-41
• /
• 1996

Linear programming has become an important tool in decision-making of modern business management. This remarkable growth can be traced to the pioneering efforts of many individuals and research organizations. The popular using of personal computers make it very easy to process those complicated linear programming models. Furthermore advanced linear programming software packages assist us to solve L.P. models without any difficult process. Even though the advanced L.P. professional packages, the needs of more detailed deterministic elements for business decisions have forced us to apply dynamic approaches for more resonable solutions. For the purpose of these problems applying to the "Mathematica" packages which is composed of mathematic tools, the simplex processes show us the flexible and dynamic decision elements included to any other professional linear programming tools. Especially we need proper dynamic variables to analyze the shadow prices step by step. And applying SAS(Statistical Analysis System) packages to the L.P. problems, it is also one of the best way to get good solution. On the way trying to the other L.P. packages which are prepared for Spreadsheets i.e., MS-Excel, Lotus-123, Quatro etc. can be applied to linear programming models. But they are not so much useful for the problems. Calculating simplex tableau is an important method to interpret L.P. format for the optimal solution. In this paper we find out that the more detailed and efficient techniques to interpret useful software of mathematica and SAS for business decision making of linear programming. So it needs to apply more dynamic technique of using of Mathematica and SAS multiple software to get more efficient deterministic factors for the sophiscated L.P. solutions.

• The Journal of Korean Association of Computer Education
• /
• v.21 no.6
• /
• pp.39-47
• /
• 2018

This paper dealt with investigating the number of minimal spanning ladders originated from ladder game and their properties as well as the related computational thinking aspects. The author modified the filtering techniques to enhance Mathematica project where a new type of graph was generated based on the algorithm using a generator of firstly found minimal spanning graph by repeatedly applying independent ladder operator to a subsequence of ladder sequence. The newly produced YC graphs had recursive and hierarchical graph structures and showed the properties of edge-symmetric. As the computational complexity increased the author divided the whole search space into the each floor of the newly generated minimal spanning graphs for the (5, 10) YC graph and the higher (6, 15) YC graph. It turned out that the computational thinking capabilities such as data visualization, abstraction, and parallel computing with Mathematica contributed to enumerating the new YC graphs in order to investigate their structures and properties.

• Structural Engineering and Mechanics
• /
• v.28 no.2
• /
• pp.221-238
• /
• 2008

Numerical solution to buckling analysis of beams and columns are obtained by the method of differential quadrature (DQ) and harmonic differential quadrature (HDQ) for various support conditions considering the variation of flexural rigidity. The solution technique is applied to find the buckling load of fully or partially embedded columns such as piles. A simple semi- inverse method of DQ or HDQ is proposed for determining the flexural rigidities at various sections of non-prismatic column ( pile) partially and fully embedded given the buckling load, buckled shape and sub-grade reaction of the soil. The obtained results are compared with the existing solutions available from other numerical methods and analytical results. In addition, this paper also uses a recently developed technique, known as the differential transformation (DT) to determine the critical buckling load of fully or partially supported heavy prismatic piles as well as fully supported non-prismatic piles. In solving the problem, governing differential equation is converted to algebraic equations using differential transformation methods (DT) which must be solved together with applied boundary conditions. The symbolic programming package, Mathematica is ideally suitable to solve such recursive equations by considering fairly large number of terms.

• Journal of the Chungcheong Mathematical Society
• /
• v.19 no.2
• /
• pp.133-139
• /
• 2006

By combining the classical Newton's method with the pseudo-secant method, pseudo-secant-Newton's method is constructed and its order and rate of convergence are investigated. Given a function $f:\mathbb{R}{\rightarrow}\mathbb{R}$ that has a simple real zero ${\alpha}$ and is sufficiently smooth in a small neighborhood of ${\alpha}$, the convergence behavior is analyzed near ${\alpha}$ for pseudo-secant-Newton's method. The order of convergence is shown to be cubic and the rate of convergence is proven to be $\(\frac{f^{{\prime}{\prime}}(\alpha)}{2f^{\prime}(\alpha)}\)^2$. Numerical experiments show the validity of the theory presented here and are confirmed via high-precision programming in Mathematica.

• Journal of the Chungcheong Mathematical Society
• /
• v.22 no.2
• /
• pp.217-227
• /
• 2009

Assuming that a given nonlinear function f : $\mathbf{R}{\rightarrow}\mathbf{R}$ has a zero $\alpha$with integer multiplicity $m{\geq}1$ and is sufficiently smooth in a small neighborhood of $\alpha$, we define extended leap-frogging Newton's method. We investigate the order of convergence and the asymptotic error constant of the proposed method as a function of multiplicity m. Numerical experiments for various test functions show a satisfactory agreement with the theory presented in this paper and are throughly verified via Mathematica programming with its high-precision computability.

• Journal of applied mathematics & informatics
• /
• v.16 no.1_2
• /
• pp.221-229
• /
• 2004

It is known that the parametric boundary equation for the main component in the Mandelbrot set represents a cardioid. We derive an epicy-cloidal boundary equation of the main component in the degree-n bifurcation set by extending the parameter which describes the cardioid in the Mandelbrot set. Computational results as well as some useful properties are presented together with the programming source codes written in Mathematica. Various boundaries are displayed for $2\leqn\leq7$7 and show a good agreement with the theory presented here. The known boundary equation enables us to significantly reduce the construction time for the degree-n bifurcation set.

• Journal of the Chungcheong Mathematical Society
• /
• v.23 no.4
• /
• pp.845-852
• /
• 2010

A nonlinear algebraic equation f(x) = 0 is considered to find a root with integer multiplicity $m{\geq}1$. A variant of Newton-secant method for a multiple root is proposed below: for n = 0, 1, $2{\cdots}$ $$x_{n+1}=x_n-\frac{f(x_n)^2}{f^{\prime}(x_n)\{f(x_n)-{\lambda}f(x_n-\frac{f(x_n)}{f^{\prime}(x_n)})\}$$, $$\lambda=\{_{1,\;if\;m=1.}^{(\frac{m}{m-1})^{m-1},\;if\;m{\geq}2$$ It is shown that the method has third-order convergence and its asymptotic error constant is expressed in terms of m. Numerical examples successfully verified the proposed scheme with high-precision Mathematica programming.

• Journal of applied mathematics & informatics
• /
• v.40 no.1_2
• /
• pp.267-281
• /
• 2022

Herein, an algorithm for efficient evaluation of oscillatory Fourier-integrals with Jacobi-Cauchy type singularities is suggested. This method is based on the use of the traditional Clenshaw-Curtis (CC) algorithms in which the given function is approximated by the truncated Chebyshev series, term by term, and the oscillatory factor is approximated by using Bessel function of the first kind. Subsequently, the modified moments are computed efficiently using the numerical steepest descent method or special functions. Furthermore, Algorithm and programming code in MATHEMATICA^® 9.0 are provided for the implementation of the method for automatic computation on a computer. Finally, selected numerical examples are given in support of our theoretical analysis.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.1
• /
• pp.35-43
• /
• 1996

A computer program for automatic deriving the symbolic equations of motion for robots using the programming language MATHEMATICA has been developed. The program, developed based on the Lagrange formalism, is applicable to the closed chain robots as well as the open chain robots. The closed chains are virtually cut open, and the kinematics and dynamics of the virtual open chain robot are analyzed. The constraints are applied to the virtually cut joints. As a result, the spatial closed chain robot can be considered as a tree structured open chain robot with kinematic constraints. The topology of tree structured open chain robot is described by a FATHER array. The FATHER array of a link indicates the link that is connected in the direction of base link. The constraints are represented by Lagrange multipliers. The parallel robot, DELTA, having three-dimensional closed chains is modeled and simulated to illustrate the approach.