# A Computational Model forq-Bernstein Quasi-Minimal Bézier Surface

Authors: Daud Ahmad; M. Khalid Mahmood; Qin Xin; Ferdous M. O. Tawfiq; Sadia Bashir; Arsha Khalid
Journal: Journal of Mathematics
 (2022)

---

## Abstract

A computational model is presented to find theq-Bernstein quasi-minimal Bézier surfaces as the extremal of Dirichlet functional, and the Bézier surfaces are used quite frequently in the literature of computer science for computer graphics and the related disciplines. The recent work [1–5] on q-Bernstein–Bézier surfaces leads the way to the new generalizations of q-Bernstein polynomial Bézier surfaces for the related Plateau–Bézier problem. The q-Bernstein polynomial-based Plateau–Bézier problem is the minimal area surface amongst all the q-Bernstein polynomial-based Bézier surfaces, spanned by the prescribed boundary. Instead of usual area functional that depends on square root of its integrand, we choose the Dirichlet functional. Related Euler–Lagrange equation is a partial differential equation, for which solutions are known for a few special cases to obtain the corresponding minimal surface. Instead of solving the partial differential equation, we can find the optimal conditions for which the surface is the extremal of the Dirichlet functional. We workout the minimal Bézier surface based on the q-Bernstein polynomials as the extremal of Dirichlet functional by determining the vanishing condition for the gradient of the Dirichlet functional for prescribed boundary. The vanishing condition is reduced to a system of algebraic constraints, which can then be solved for unknown control points in terms of known boundary control points. The resulting Bézier surface is q-Bernstein–Bézier minimal surface.


---

## Body

## 1. Introduction

We observe the tendency of nature to be minimum in constructing the path predicting models of different objects around us. For instance, the principles of flexibility involves the optimization of processes based on the observation that a business model should change and grow with time. Other examples are the optimal designs in solid and fluid mechanics, electromagnet of an electromagnetic particle, and the mathematical models for gravity theory are the phenomena in nature that access to physical systems for optimization. These phenomena are analysed in the optimization theory, which include the calculus of variations, game theory, decision theory, linear programming, control theory, network analysis, and Markov chains.One of the active research areas in the optimization theory [1] is calculus of variations in which we try to find a best suitable function, which could be a curve or a surface subject to certain constraint usually expressed in the integral form called the functional, which is in fact a critical point of that functional, a maxima or minima of the functional under consideration based on the mathematical analysis of minimization principals on multidimensional function spaces. The constraints can be in the form of integrals (functionals based on unknown function of one or more parameters), differential constraints (they arise for instance in case of rolling wheels, flying aircrafts, in situations arising from dynamics of any mechanical system, and so on) or the algebraic structure in constraints. This work is related to finding the surfaces subject to constraints in the form of certain integral called the functional, and one of such functionals is Dirichlet functional with appropriate boundary conditions. The Dirichlet functional being the integral of non-negative quantity is itself non-negative. The variational problem of solving the Dirichlet functional for the minimal vector function xs,t subject to certain boundary conditions is the Euler–Lagrange partial differential equation ∇2xs,t=0, which is Laplace equation in this case. Such a vector function xs,t has minimal Dirichlet energy and is said to be harmonic. This can be done by finding the extremal of area functional of a surface xs,t, but it involves square root of the integrand, and instead, we can use variational approach for finding such a surface as the extremal of Dirichlet functional and we obtain a quasi-minimal surface. A minimal surface is a surface for which mean curvature of the surface vanishes for all possible parameterizations of the same surface.In curve theory, the familiar examples are that of shortest distance between two points in a plane, which is a straight line in case when no constrains are involved; however, in case of given constraints, the shortest path or the distance between two points is usually called a geodesic on a surface, and the related ancient well-known problem is the brachistochrone problem. In optics, one of the related problems is the Fermat’s principle (1657), to find the path in least possible time. Fermat explained that the laws of geometric optics, lens design such as reflection, refraction, focusing, and aberrations could be explained from the geometric and analytical properties of this principle. The natural generalization of the minimal curve or geodesic problem is the minimal surface problem [2, 3]. It consists of searching the surface x⊂ℝ3 with possible bounded area spanned by a closed contour Γ. This essentially means setting the surface area functional A=∫ℛdx as an objective of extremization functional over all possible x⊂ℝ3 with prescribed contour Γ to achieve a minimal area surface. Minimal surface problem is referred to as Plateau–Bézier problem [4, 5] in the honour of physicist of Belgian Joeseph Plateau, who established that minimal surfaces can be associated with soap films spanned by closed wire frames in 1849. He explained further through his experiments that a minimal surface can be achieved in the form of a thin soap film spanned by a wire frame by immersing it into soapy water and displacing it back cautiously. The soap film itself is a surface of minimal area, and the wire frame in shape of closed contour Γ serves as a spanning curve.The soap films and bubbles have been a source of fascination from aesthetical, physical, and mathematical point of view that is why it is an active field of research for several hundred years, finding its applications in different disciplines of science. The earliest significant work is that of Euler, who in 1744, while searching for rotational surfaces of minimal area, proved that a minimal surface is planer if and only if its Gaussian curvature is identically equal to zero and a minimal surface is always locally saddled-shaped. Later in 1762, Lagrange considered the minimal surface problem by deriving associated Euler–Lagrange equations, which is a quadratic partial differential equationpde. The solutions of such partial differential equations are the functions that minimize a given functional. However, he was not able to find out a general solution of the equation other than a plane. In 1776, Meusnier illustrated the condition for minimality for a specific form of a regular surface, known as Monge’s patch. The outcome of the minimal condition for the Monge’s patch is a quasilinear, quadratic partial differential equation derived earlier by Lagrange. He also showed that catenoid and helicoid satisfy Euler–Lagrange equations too and proved that mean curvature vanishing surfaces are in fact minimal surfaces. In 1830, Heinrich Scherk used the Euler–Lagrange equation to find the nontrivial examples of complete minimal surfaces. The first golden era of minimal surface theory began in mid-19th century after the discovery made by J. Plateau. In 1867, Schwarz [6] found the solution of Plateau–Bézier problem for a general quadrilateral by an appeal to methods in complex analysis. Specifically, he discovered CLP (crossed layers of parallels) surface Schwarz T for tetragonal, Schwarz D for diamond, Schwarz H for hexagonal, and Schwarz P for primitive surfaces. Weierstrass and Enneper [4] developed representation formulas for minimal surface in conformal parametrization as a pair of holomorphic and meromorphic functions. Weierstrass and Enneper parametric representation of minimal surfaces serves as bridge between the geometry and the complex analysis [7]. Later on, achievements of Jesse Douglas [8], Rado [9], Tonelli [10], Courant [7, 11], Wang [12] (minimal surfaces based modelling of atoms and molecules), and others contributed many revolutionary results in the subsequent years in minimal surfaces spanned by given boundary curves. General solution of Plateau–Bézier problem remained challenging for it involves extremization of area functional that is highly nonlinear in its integrand. In 1931, American mathematician J. Douglas [8] and, in 1933, Hungarian Tibor Radó [9] showed the existence of a minimal surface spanned by a closed curve by exchanging the area integral by rather a simpler integral the so-called Douglas-Dirichlet integral. Revolutionary achievements in subsequent years by Courant [11], Max. Shiffman [13], Morse [14], Tompkins [14], Robert-Osserman [15], and others are worth mentioning that contributed many significant results to theory of minimal surfaces. A detailed account of related results of measure theory and differential geometry for Plateau–Bézier problem can be seen in reference [16]. Struwe [17] points out the parametric minimal surfaces and the developments in Plateau–Bézier problem with a concise treatment of constant mean curvature surfaces. Harrison [18] considered Plateau–Bézier problem as test problem for differential chains Theory.The properties of minimal surfaces provide a support for shape modelling and shape fairing techniques. Minimal surface theory is an area of contemporary research and is of great importance in engineering design, computer-aided geometric designs (CAGD), architectural design, and biology that includes foams, domes, and cell membranes, and so on. In 2004, Monterde [19] used tensor product Bézier surfaces to find an approximate minimal surface giving a continuous surface spanned by a given boundary curve, called the Plateau–Bézier problem. They extremized the Dirichlet functional by finding the gradient of the functional with respect to the control points. Monterde and Ugail [20], in 2006, introduced a different functional, namely, the general bi-quadratic functional, which could be reduced to a functional introduced by Farin and Hansford [21], standard biharmonic functional introduced by Schneider and Kobbelt [22], or Bloor and Wilson’s modified biharmonic functional [23]. In 2011, Monterde [24] constructed triangular Bézier patches by introducing a variational approach, and in 2015, Monterde and Ugail [25] found the quasi-harmonic surface as the extremal of quasi-harmonic functional in an attempt to solve the related Plateau–Bézier problem. Farin and Handford [26] obtained combinatorial relations between q-Genocchi numbers and polynomials with weight α and β and integral representation of weighted q-Bernstein polynomials. Bernstein polynomial is an important class of polynomials related to special polynomials like Bernoulli polynomials, Euler polynomials, and so on, in the theory of analytical numbers, and they find their applications in designing smooth curves and surfaces called Bézier curves and Bézier surfaces. They play a significant role for constructing the surfaces of different shapes and desired characteristics that depend on the choice of Bernstein polynomials, namely, the classical Bernstein polynomials, shifted knots Bernstein polynomials, and q-Bernstein polynomials, generally called the modified Bernstein polynomials. An extension of q-Bernstein polynomials in surface theory is p,q-Bernstein polynomials, basic properties, and generating functions for Bernstein polynomial-related results with the help of p,q-calculus, which can be seen in references [27, 28]. Araci [29] and Jang et al. [30] study the q-analogue of Euler numbers and polynomials naturally arising from the p-adic fermionic integrals on Zp and investigate some properties for these numbers and polynomials. We use one of these Bernstein polynomials, recently introduced by Kim [31] to formulate the related Bézier surfaces for extreme values of Dirichlet functional resulting in a quasi-minimal surface.The minimal surfaces occur naturally in the developing fields, namely, CAGD, (computer-aided geometric design), CG (computer graphics), and CAD (computer-aided design) in mathematical. The minimal surfaces are known through the variation of the respective area functional. The area functional involves square-root in its integrand, and it is not always possible to find the solution of the PDE as the outcome of the vanishing condition of its gradient for a surface to be a minimal surface. However, the variation of other functionals is more useful for some specific desired features of a surface based on the restriction or the constraint. These constraints are usually an integral to be minimized, for instance, an energy integral, for example, the Dirichlet functional,rms of mean curvature, quasi-harmonic functional, and biharmonic functional. These variational techniques can be applied to a famous class of surfaces known as Bézier surfaces, which have their own importance for their useful properties and applications in CAGD, CG, and CAD. There are two main categories of Bézier patches, namely, the rectangular and triangular Bézier patches. For a given Bézier control net of points Pjk, a rectangular Bézier surface is defined as(1)xs,t=∑j=0n∑k=0mℬjnsℬkmtPjk,where ℬjns are the Bernstein polynomials. The rectangular Bézier surfaces (1) are based on the univariate Bernstein polynomials:(2)ℬjns=njsj1−sn−j,which named after the Sergei Natanovich Bernstein [32], where(3)nj=n!j!n−j!are the usual binomial coefficients. The triangular Bézier patches Ss,t,r=∑j+k+l=nPjklBj,k,lns,t,r are expressed in terms of bivariate Bernstein polynomials Bj,k,lns,t,r=nj,k,lsjtkrl. One of the problems in CAGD, for example, is to construct a surface of least area amongst all the surfaces spanned by the prescribed boundary, called the Plateau–Bézier problem. A surface of least area can be determined as one of the critical points of the area functional by finding constraints on the interior control points as the outcome of vanishing condition for the gradient of area functional; however, the area functional involves square root in its integrand and makes it difficult to find the surface as one of its critical points. The area functional can be replaced by some other suitable energy functional having similar properties as that of area functional or may be some other desired feature of the surface. It is then possible to find the gradient of this chosen energy functional (constraint functional in the form of an integral). The vanishing condition of the gradient of the functional gives linear algebraic constraints on the unknown control points in terms of known boundary control points. One of the widely used class of surfaces is known as the Bézier surfaces based on Bernstein polynomials as has been done by Monterde [19,20], and he found the quasi-minimal Bézier surface as the extremal of Dirichlet functional, other related works for ansatz method, and for the vanishing condition of certain energy functionals for the Bézier surfaces based on modified Bernstein polynomials and Coons patch one can see Ahmad et al. [33–38].As mentioned above, one of the widely used restrictions is to find out the minimal Bézier surface as the extremal of various energy functionals by the vanishing condition of gradient of such a functional. We choose one of such restrictions, namely, the Bézier surface withq-Bernstein polynomials, and find the corresponding minimal surface as the extremal of Dirichlet functional by finding the vanishing condition of gradient of this Dirichlet functional for Bézier surface with q-Bernstein polynomials, which gives us constraints on the interior control points in terms of boundary control points for which the new surface is minimal. It is to be remarked that the minimizers of a certain chosen functional are possibly the candidate functions for a quasi-minimal surface obtained from the zero-functional gradient of the surface xs,t, which is equivalent to finding the extremal points of a function of several variables in calculus. The functional gradient equated to zero gives us linear constraints on the unknown interior control points depending on the known boundary control points. The surface spanned by these new interior control points with prescribed border is minimal, which we refer to as the q-Bernstein quasi-minimal Bézier surface. The applications of these surfaces are in building and material sciences, and formation of a surface in CAGD, for example, in modelling, a problem related to protest investigation and mechanics of cell material. Bézier models can assist in designing the computational-based prediction models, for instance, for human-engineered representations for their possible applications for machine learning capabilities [39–42].Our aim is to solve Plateauq-Bernstein–Bézier problem, which is to find the q-Bernstein–Bézier surface as one of the critical points of Dirichlet functional instead of area functional, and the resulting surface has minimal area from amongst all the possible q-Bernstein–Bézier surfaces constructed from the given prescribed boundary control points. The q-Bernstein–Bézier surface can be written as follows:(4)xs,t=∑j=0m∑k=0nQj,qmsQk,qntPjk,where Qjnsq are the q-Bernstein polynomials introduced by Kim [31] in 2011. The Kim’s q-Bernstein polynomials are different from the Philips version of q-Bernstein polynomials (1997) [26] that depend on q-integers, and later in 2003, Oruc and Phillips obtained related Bézier curves [43] and some interesting properties. We intend to find the Kim’s q-Bernstein quasi-minimal Bézier surface as the extremal of the Dirichlet functional.(5)Dxs,t=12∫ℛxs2+xt2dudv,for Kim’s q-Bernstein–Bézier surface (equation (4)). The vanishing condition of a functional gradient generates a system of algebraic conditions on the unknown inner control points as boundary control points.The rest of the article is organized as follows: in the forthcoming Section2, we give few definitions of differential geometry-related quantities for curves and surface, the minimal surfaces, q-Bernstein polynomials, and Kim type [31] q-Bernstein–Bézier surfaces, specific classes of surfaces, their basic construction schemes, properties of the general Bézier curves with kim operator [31], the generalized rational Bézier surfaces, and the q-Bernstein–Bézier surfaces. In Section 3, we develop a technique to find the quasi-minimal surfaces corresponding to q-Bernstein–Bézier surfaces. In the same section, we have included the illustrative examples for the bi-quadratic and bi-cubic q-Bernstein–Bézier surfaces as the application of the technique developed. Finally, Section 4 includes the final remarks and the future prospects of the work.


## 2. Preliminaries

In this section, we give few geometric quantities and basics of curve and surface theory,q-Bernstein polynomials, classical Bézier and q-Bézier curves and surfaces, derivatives, integrals, mean, and Gaussian curvature of q-Bézier surfaces.Definition 1.
The area functional of a surface.
A surfacexs,t of minimal area is a surface that locally minimizes its area. This is equivalent to having zero mean curvature:(6)H=12Eg−2Ff+GeEG−F2.
For every parametrization, which is direct consequence of minimizing the area functional,(7)AP=∫ℛxs∧xtdsdt=∫ℛEG−F21/2dsdt,where E=xs,xs,F=xs,xt, and G=xt,xt are coefficients of the quadratic form Is,t=dxs,t,dxs,t, for dxs,t=xsds+xtdt, a 1−1 linear mapping of vectors ds,dt onto dxs,t=xsds+xtdt, which lies in the tangent plane. The quadratic form Is,t, called first fundamental form, is usually written as ds2=Es,tds2+2Fs,tdsdt+Gs,tdt2 for the surface xs,t.Definition 2.
Bernstein polynomials, Bézier curves, and surfaces.
Blending functions are used to define the curves and surfaces in parameterized form, and they effect the type and the shape of curves and surfaces. A set of points used to produce a curve or a surface is called a set of control points for that curve or the surface. If we denote blending functions byfjs and the control net of points by Pj, then the parameterized form of the curve for 0≤s≤1 is xs=∑j=0nfjsPj. Similarly for a net of control points Pjk for j,k=0,…,n and the blending functions fjks,t, the parameterized form of the surface is xs,t=∑j,k=0nfjks,tPjk, where 0≤s,t≤1. For a given continuous function f on the interval 0,1, the Bernstein polynomial is expressed as(8)Bnf;u=∑j=0nfjnBjnu=∑j=0nfjnnjuj1−un−k,where n is a positive integer. For each function f, the equation (8) results in a sequence of Bernstein polynomials. The continuity of the function f assures the uniform convergence of Bernstein polynomials to the function f, and this means that on the interval 0,1, limn⟶∞Bnf;u=f. It can be seen from (8) that for u=0, Bnf;0=f0, and for u=1, Bnf;1=f1. For f=1, (8) reduces to Bnf;u=Bn1;u=1, whereas for f=u, (8) reduces to Bnf;u=Bnu;u=u, and for f=u2, it is Bnf;u=Bnu2;u=n−1u2/n+u/n and so on. Bernstein [32] introduced the so-called the Bernstein polynomials, a linear combination of functions ℬjnu, ℬnu=∑j=0nβjℬjnu, also called polynomials in Bernstein form, with βk as the Bernstein coefficients or Bézier coefficients ℬknu (equations (2) and (3)), and the Bernstein operator of order n for k,n∈z0≤k≤n. For instance, the Bernstein polynomials of degree n=5, namely, ℬ05u,ℬ15u,ℬ25u,ℬ35u,ℬ45u, and ℬ55u, are shown in Figure 1.
A Bézier curve is a parametric curve, which is used in computer graphics and related fields [2,19]. The Bézier curve depends on Bernstein polynomials, which are called the blending functions or the basis of Bézier curve with a set of n+1 control points (also called Bézier points) denoted by P0,P1,P2,…,Pn. A Bézier curve of degree n is given in the form xs=∑j=0nℬjnsPj for the Bernstein polynomials ℬjns (equation (2)) of degree n for u∈0,1. Bézier surfaces xs,t (equation (1)) are the higher dimension generalization of Bézier curves for a given set of n+1,m+1 control points Pjkj,k=0m,n for the blending functions ℬjnsℬkmt=ℬj,kn,ms,t:ℛ2⟶ℛ, where ℬjns and ℬkmt are Bernstein basis functions given by the equation (2) for 0≤s,t≤1 (Figure 2).Figure 1
Polynomials in Bernstein form,ℬ05u,ℬ15u,ℬ25u,ℬ35u,ℬ45u,ℬ55u, are shown.Figure 2
q-Bernstein polynomials of various degrees are shown.Definition 3.
q-Bernstein–Bézier curves and surfaces.
Aq-Bernstein–Bézier surface (4) is a Bézier surface based on q-Bernstein polynomials, where q-Bernstein polynomials Qj,qms, the function of curve parameter s, serve as the blending functions for the Bézier surface, taken in the form,(9)Qj,qms=mjsqj1−sq−1−j+m,sqj=1−qs1−qj,where mj are the usual binomial coefficients, sqj is the Kim’s q-Bernstein operator (for detail, see Equation (1.7) of reference [31]) for 0≤q≤1. Note that limq⟶1sq=limq⟶11−qs/1−q=s and limq⟶11−sq−1=limq⟶11−q−11−s/1−q−1=1−s, so that the q-Bernstein polynomials Qj,qms reduce to the classical Bernstein polynomials ℬkns. In particular, the q-Bernstein polynomials of degrees n = 1, 2, 3 are shown in Figure 2. Plugging the q-Bernstein polynomials given by equation (9) for the surface parameters s,t in the equation (4), the Bézier surface with q-Bernstein polynomials may be rewritten as(10)xs,t=∑j,k=0m,nmjnksqjtqk1−sq−1m−j1−tq−1n−kPjk.
The bi-quadraticq-Bernstein–Bézier surface and the bi-cubic q-Bernstein–Bézier surface obtained from the equation (10) for q=0.2 for the prescribed are shown in Figure 3.Figure 3
Bi-quadratic and bi-cubicq-Bernstein–Bézier surfaces for the prescribed border.Definition 4.
Partial derivatives ofq-Bernstein–Bézier surface with respect to the control points.
The partial derivatives∂xs/∂xpqa=∂/∂s∑j,k=0m,nQj,qmsQk,qnt∂Pjk/∂xpqa of tangent vectors xs to the coordinate curves on the Bézier surface with q-Bernstein functions based on the Kim operator [31] with respect to the components xpqap=0,1,2,…,m,q=0,1,2,…n,a=1,2,3 of control points Ppq, where Pjk=xjk,yjk,zjk=xjk1,xjk2,xjk3 are the control points as mentioned above, ∂Pjk/∂xpqa is one of the standard basis vectors ea for p=j,q=k, and otherwise, for p≠j or q≠k, ∂Pjk/∂xpqa is zero vector, which can be written as(11)∂xs∂xjka=∂∂sQj,qmsQk,qntea.
Theq-Bernstein polynomial [31], for m≥j, Qj,qms, is defined as(12)Qj,qms=mjsqj1−sq−1m−j,and Qj,qms=0; otherwise, where mj is the usual binomial coefficient and(13)sqj=1−qs1−qj,1−sq−1m−j=1−q−11−s1−q−1m−j,the partial derivative of q-Bernstein polynomial Qjmsq with respect to s is given by(14)∂∂sQj,qms=mj−jsqj−11−sq−1m−j+m−jsqj1−sq−1m−j−1qslogq1−q.Definition 5.
Integral ofq-Bernstein polynomial
Theq-Bernstein polynomials (9) for the Kim operator sqj can be written as(15)Qj,qms=mj1−qs1−qj1−q−11−s1−q−1m−j.
The product of twoq-Bernstein polynomials with same parameter s is given by(16)Qj,qmsQk,qns=mjnk1−qs1−qj1−qs1−qk1−q−11−s1−q−1m−j1−q−11−s1−q−1n−k,and it ((16) can be written as(17)Qj,qmsQk,qns=mjnk1−qs1−qj+k1−q−11−s1−q−1m+n−j+k.
However, thej+k-th term of q-Bernstein polynomial of degree m+n is given by(18)Qj+k,qm+ns=m+nj+k1−qs1−qj+k1−q−11−s1−q−1m+n−j+k.
Thus, the product of twoq-Bernstein polynomials (17) can be written in the form(19)Qj,qmsQk,qns=mjnkm+nj+k−1Qj+k,qm+ns.
For the integral ofq-Bernstein polynomial based on Kim operator and the integral of product of two q-Bernstein polynomials having the same parameter s, we state the result below.Theorem 1.
The integral ofq2sQj,qms, where Qj,qms is the q-Bernstein polynomial for 0≤i≤m−j and 0≤p≤i+j, with respect to the parameter s, is given by(20)∫01q2sQj,qmsds=∑i=0m−j∑p=0i+j−1i+pmjm−jii+jp1−q−i−jqp+2−1logqp+2,where(21)Qj,qms=mjsqj1−sq−1m−j,andlimp⟶0qp−1logqp=1.Proof.
In order to find out the above integral, we shall use the following result that helps us to write the product ofsqj and 1−sq−1m−j in the following form [26]:(22)q2ssqj1−sq−1m−j=q2s∑i=0m−j−1im−jisqi+j,which can then be used to establish the given result. Let us write the integral of sqj as follows:(23)∫01q2ssqjds=∫01q2s1−qs1−qjds=1−q−jIj,where(24)Ij=∫01q2s1−qsjds.
The binomial expansion of1−qsj is given by(25)q2s1−qsj=q2s∑k=0j−1pjpqsp=∑p=0j−1pjpqp+2s,which enable us to write the integral (24) in the form(26)Ij=∑p=0j−1pjp∫01qp+2sds,and thus, the integral (26) reduces to(27)Ij=∑p=0j−1pjpqp+2slogqp+201=∑p=0j−1pjpqp+2−1logqp+2,where(28)I0=1.
Plugging the (27) in (23), we get(29)∫01sqjds=1−q−j∑p=0j−1pjpqp+2−1logqp+2,and thus,(30)∫01sqi+jds=1−q−i−j∑p=0i+j−1pi+jpqp+2−1logqp+2.
The integral of (22)(31)∫01sqj1−sq−1m−jds=∑i=0m−j−1im−ji∫01sqi+jdsalong with (30) is given by(32)∫01sqj1−sq−1m−jds=∑i=0m−j−1im−ji1−q−i−j∑p=0i+j−1pi+jpqp+2−1logqp+2,which reduces to(33)∫01sqj1−sq−1m−jds=∑l=0m−j∑p=0i+j−1i+p1−q−i−jm−jii+jpqp+2−1logqp+2,where(34)limp⟶0qp−1logqp=1.
Hence, the integral ofq-Bernstein polynomial is given by (20). The results in the corollaries below are given for later use in Section 3.Corollary 1.
The integral ofq-Bernstein polynomial with respect to the parameter s when m is replaced by 2m−2 and j by j+l−1 is given by(35)∫01q2sQj+l−1,q2m−2sds=∑i=02m−j−2∑p=0i+j+l−1−1i+p2m−2j2m−j−l−1ii+j+l−1p1−q−i−j−l+1qp+2−1logqp+2.Corollary 2.
The integral ofq-Bernstein polynomial with respect to the parameter s when m replaced by 2m−2 and j by j+l is given by(36)∫01q2sQj+l,q2m−2sds=∑i=02m−j+l−2∑p=0i+j+l−1i+p2m−2j+l2m−j+l−2ii+j+lp1−q−i−j−lqp+2−1logqp+2.Corollary 3.
The integral ofq-Bernstein polynomial with respect to the parameter s when m replaced by 2m and j by j+l is given by(37)∫01q2sQj+l,q2msds=∑i=02m−j+l∑p=0i+j+l−1i+p2mj+l2m−j+lii+j+lp1−q−i−j+lqp+2−1logqp+2.Corollary 4.
The integral ofq-Bernstein polynomial with respect to the parameter t when m=2n−2,j=k+r−1 is given by(38)∫01q2tQk+r−1,q2n−2tdt=∑i=02n−k−r−1∑p=0i+k+r−1−1i+p2n−2k+r−12n−k−r−1i×i+k+r−1p1−q−i−k−r+1qp+2−1logqp+2.Corollary 5.
The integral ofq-Bernstein polynomial with respect to the parameter t when m is replaced by 2n−2 and j by r+k is given by(39)∫01q2tQk+r,q2n−2tdt=∑i=02n−2−k+r∑p=0i+k+r−1i+p2n−2k+r2n−2−k+rii+k+rp1−q−i−k+rqp+2−1logqp+2.Corollary 6.
The integral ofq-Bernstein polynomial with respect to the parameter t when m=2n,j=r+k is given by(40)∫01q2tQk+r,q2ntdt=∑i=02n−k+r∑p=0i+k+r−1i+p2nk+r2n−k+rii+k+rp1−q−i−k+rqp+2−1logqp+2.Definition 6.
Mean and Gaussian curvature ofq-Bernstein–Bézier surface.
For the mean and Gaussian curvature of theq-Bernstein–Bézier surface, (10), we find the fundamental coefficients E,F,G and e,f,g of the q-Bernstein surface with Kim operator xqs,t, which requires computation of the partial derivatives of the surface xqs,t with respect to the surface parameters s,t. Note that the partial derivatives of the surface with respect to its surface parameters appear as terms involving the derivatives of q-Bernstein polynomials. The fundamental coefficients of the q-Bernstein–Bézier surface (10) help us to compute the mean curvature and the Gaussian curvature of the surface. The mean curvature (6) of the surface is(41)H=∑j,k=0m,njkm−jq1−t1−q1−t−j2k4qt1−qtm2m−jq1−t1−q1−t−jm2qt1−qtjkm−jq1−s1−q1−s−j2k4qs1−qsm2m−jq1−s1−q1−s−jm2qs1−qs1−q1−s1−qm−j1−q1−t1−qm−j1−qs1−qj1−qt1−qjlogq1−qPjk.
The Gaussian curvatureK=LN−ℳ2/EG−F2 of the surface is given by the following expression:(42)K=∑j,k=0m,njkm−jq1−t1−q1−t−j2k4qt1−qt2m2m−jq1−t1−q1−t−jm2qt1−qt2jkm−jq1−s1−q1−s−j2k4qs1−qs2m2m−jq1−s1−q1−s−jm2qs1−qs21−q1−s1−qm−j1−q1−t1−qm−j1−qs1−qj1−qt1−qjlogq1−qPjk.
We have discussed restricted class of surfaces, namely, the Bézier surfaces withq-Bernstein polynomials, the basic construction scheme of these surfaces, their properties based on kim operator [31], mean curvature, and Gaussian curvature of the surfaces. In the section below, we come up with a scheme for finding the quasi-minimal q-Bernstein–Bézier surface as the extremal of Dirichlet functional.


## 3. Quasi-Minimalq-Bernstein–Bézier Surfaces as the Extremal of Dirichlet Functional

This section is devoted to the problem of finding the quasi-minimalq-Bernstein–Bézier surface as the extremal of Dirichlet functional by determining the constraints on the interior control points of a given q-Bernstein–Bézier surface for the prescribed border. For this, we shall find the gradient of the Dirichlet functional (5) for the q-Bernstein–Bézier surface (10) with the prescribed border in terms of its boundary points and equate it to zero that results in a system of linear algebraic constraints. These algebraic constraints can be solved for interior control points in terms of boundary control points for a specific mesh structure. When these new interior control points are plugged in the q-Bernstein–Bézier surface, the emerging surface is a quasi-minimal surface. For illustration of the given scheme, we have included bi-quadratic and bi-cubic quasi-minimal q-Bernstein–Bézier surfaces in the section immediately after exploring the scheme.Following proposition enables us to write the linear algebraic constraints satisfied by the interior control points ofq-Bernstein–Bézier surface meshed with the boundary control points as the outcome of vanishing condition of gradient of the Dirichlet functional for the q-Bernstein–Bézier surface. The interior control points worked out from these linear algebraic constraints together with known boundary control points help us to construct the desired quasi-minimal q-Bernstein–Bézier surface. The proposition is given below.Proposition 1.
A control netP=Pjkj,k=0n,m of q-Bernstein–Bézier surface of prescribed border is the extremal of Dirichlet functional if the control points Pjkj,k=0n,m (for 0<q<1) satisfy the following constraint equation:(43)mlogq1−q2∑i,j=0m−1,nξj−1,lm−1,m−1δj+l−12m−2s−ξj,lm−1,m−1δj+l2m−2sξk,rn,nδk+r2ntΔ10Pij+nlogq1−q2∑i,j=0m,n−1ξk−1,rn−1,n−1δk+r−12n−2t−ξr,kn−1,n−1δr+k2n−2tξl,jm,mδl+j2msΔ01Plr=0,where(44)δj,qms=∑i=0m−j∑p=0j+l−1i+pmjm−jij+ip1−q−j−iqp+2−1logqp+2.Proof.
LetP=Pjkj,k=0n,m be the control net of a q-Bernstein–Bézier surface with the corresponding patch defined in equation (4). The Dirichlet functional [44], equation (5), can be written in the following convenient form:(45)DP=12∫ℛxs,xs+xt,xtdsdt.
Let us calculate the gradient of Dirichlet functional with respect to the coordinatesxjka of the control point Pjk=xjk1,xjk2,xjk3. For a∈1,2,3, j=1,2,…,m−1 and k=1,2,…,n−1, the gradient of the Dirichlet functional (45) can be written as(46)∂DP∂xjka=∫ℛ∂xs∂xjka,xs+∂xt∂xjka,xtdsdt.
We need to find the partial derivatives∂xs/∂xjka and ∂xt/∂xjka of the (46) for the q-Bernstein–Bézier surface (4). For the q-Bernstein–Bézier surface given by the equation (4),(47)∂xs∂xpqa=∂∂s∑j,k=on,mQj,qmsQk,qnt∂Pjk∂xpqa,where Pjk=xjk,yjk,zjk=xjk1,xjk2,xjk3≡xjka, a=1,2,3, and ∂Pjk/∂xpqa reduces to one of the standard basis vectors ea for p=j,q=k; otherwise, for p≠j or q≠k, ∂Pjk/∂xpqa is zero vector, and the equation (47) reduces to(48)∂xs∂xjka=∂∂sQj,qmsQk,qntea.
Theq-Bernstein polynomials Qj,qms (Kim [31]) are given by the following expression:(49)Qj,qms=mjsqj1−sq−1m−j,where mj is the usual binomial coefficients,(50)sqj=1−qs1−qj.
By virtue of the equation (50) for sqj, the q-Bernstein polynomial Qj,qms (49) can be written in the form(51)Qj,qms=mj1−qs1−qj1−q−11−s1−q−1m−j,the partial derivative of Qj,qms (51) with respect to one of (the surface) parameter s,(52)∂∂sQj,qms=mj∂∂s1−qs1−qj1−q−11−s1−q−1m−j+1−qs1−qj∂∂s1−q−11−s1−q−1m−j,for(53)∂∂s1−qs1−qj=−jqs1−q1−qs1−qj−1logq,and(54)∂∂s1−q−11−s1−q−1m−j=m−jq−11−s1−q−11−q−11−s1−q−1m−j−1logq−1,can be written as(55)∂∂sQj,qms=mj−jqs1−q1−qs1−qj−1logq1−q−11−s1−q−1m−j+1−qs1−qj×n−jq−11−s1−q−11−q−11−s1−q−1n−j−1logq−1m−j,which is then reduced to the form (using the property that logq−1=−logq).(56)∂∂sQj,qms=mj−jqs1−q1−qs1−qj−11−q−11−s1−q−1m−jlogq−m−jq−11−s1−q−11−qs1−qj1−q−11−s1−q−1m−j−1logq.
Taking out the common factors from the above expression, we get(57)∂∂sQj,qms=mj1−qs1−qj1−q−11−s1−q−1m−j−jqs1−q1−qs1−q−1−m−j×q−11−s1−q−11−q−11−s1−q−1−1logq,which can be further simplified in the form(58)∂∂sQj,qms=mj1−qs1−qj1−q−11−s1−q−1m−j−jqs1−qs−m−jq−11−s1−q−11−slogq.
By virtue of the equation (50), last equation (21) can be written as(59)∂∂sQj,qms=mjsqj1−sq−1m−j−jqs1−qs−m−jq−11−s1−q−11−slogq.
Note that for−jqs/1−qs−m−jq−11−s/1−q−11−s=−jsq−1×qs/1−q−m−j1−sq−1−1×q−11−s/1−q−1 of the above (59) can be reduced to the following form:(60)∂∂sQj,qms=mjsqj1−sq−1m−j−jsq−1×qs1−q−m−j1−sq−1−1×q−11−s1−q−1logq,in which q−11−s/1−q−1=qs−1/1−q−1=qsq−1/1−q−1=qsq−1/q−1q−1=qs/q−1, and thus, we can be written in the equation (60) as(61)∂∂sQj,qms=mjsqj1−sq−1m−j−jsq−1×qs1−q−m−j1−sq−1−1×qsq−1logq,and taking out the common factor from (61) yields(62)∂∂sQj,qms=mjsqj1−sq−1m−j−jsq−1+m−j1−sq−1−1qs1−qlogq.
Combining the termssqj1−sq−1m−j with the terms inside the parentheses of above (62), which is(63)∂∂sQj,qms=mj−jsqjsq−11−sq−1m−j+m−jsqj1−sq−1m−j1−sq−1−1qs1−qlogq,
equation (63) can be written in the following form:(64)∂∂sQj,qms=mj−jsqj−11−sq−1m−j+m−jsqj1−sq−1m−j−1qslogq1−q.
Using the (49), we note that(65)Qj−1,qm−1s=m−1j−1sqj−11−sq−1m−1−j−1,or(66)m−1j−1−1Qj−1,qm−1s=sqj−11−sq−1m−j.and in the similar way, we write(67)m−1j−1Qj,qm−1s=sqj1−sq−1m−1−j.
Plugging the (66 and 67) in the (64) allows us to write(68)∂∂sQj,qms=mj−jm−1j−1−1Qj−1,qm−1s+m−jm−1j−1Qj,qm−1sqslogq1−q.
It is to be noted that in the above (68), the binomial expressions jmjm−1j−1−1=m and m−jmjm−1j−1=m, and thus, the (68) for partial derivative of q-Bernstein polynomial can be written in lower degree polynomials given by(69)∂∂sQj,qms=mQj,qm−1s−Qj−1,qm−1sqslogq1−q.
Substituting (69) in equation (48), we obtain the relation for the partial derivative of the q-Bernstein–Bézier surface with respect to the surface parameter s as follows:(70)∂xs∂xjka=mqslogq1−qQj,qm−1s−Qj−1,qm−1sQk,qntea.
Adopting the same procedure as done for the derivation of (70), we can find the partial derivative ∂xt/∂xjkaq that, with respect to other surface parameter t,(71)∂xt∂xjka=∂∂xjka∂x∂t=∂∂t∂x∂xjka,reduces to(72)∂xt∂xjka=∂∂tQj,qmsQk,qntea=Qj,qms∂∂tQk,qntea,and we can write the partial derivative ∂/∂tQknt in the similar way as we worked out (69). Thus, the equation (70) can be written by appropriately replacing s, m, and j by t, n, and k, respectively, in the above (70) to obtain(73)∂xt∂xjka=nqtlogq1−qQj,qmsQk,qn−1t−Qk−1,qn−1tea.
Substituting equations (70) and (73) in equation (46) to get(74)∂DP∂xjka=∫ℛmlogq1−qqsQj,qm−1s−Qj−1,qm−1sQk,qntea,xs+nlogq1−qqtQj,qmsQkn−1tq−Qk−1n−1tqea,xtdsdt.
In above equation (74), we need now the partial derivatives xss,tq and xts,tq of q-Bernstein–Bézier surface xs,tq. For this purpose, we use the following expression for the q-Bernstein–Bézier surface xs,tq,(75)xs,tq=∑l=0m∑r=0nQl,qmsQr,qntPlr,and then find its partial derivative with respect to the surface parameter s and t. (This can be established by finding directly the partial derivative of Qkns.) The partial derivative of q-Bernstein–Bézier surface xs,tq with respect to its parameter s is(76)xss,tq=∑l=0m∑r=0n∂∂sQl,qmsQr,qntPlr.
Plugging the value of∂/∂sQl,qms from (69) in above (76) to obtain(77)xss,tq=mqslogq1−q∑l=0m−1∑r=0nQl,qm−1s−Ql−1,qm−1sQrntPlr,which can be written as(78)xss,tq=−mqslogq1−q∑l=0m−1∑r=0nQl,qm−1sQr,qntΔ10Plr.
By the symmetry of the resulting (78), we can write the partial derivative of q-Bernstein–Bézier surface xs,tq with respect to surface parameter t as(79)xts,tq=−nqtlogq1−q∑l=0m∑r=0n−1Ql,qmsQr,qn−1tΔ01Plr.
Substituting (78) and (79) in (74), we get(80)∂DP∂xjka=∫ℛmlogq1−q2q2sQj−1,qm−1s−Qj,qm−1sQk,qntea,∑l=0m−1∑r=0nQl,qm−1sQr,qntΔ10Plr+nlogq1−q2q2tQj,qmsQk−1,qn−1t−Qk,qn−1tea,∑l=0m∑r=0n−1Ql,qmsQr,qn−1tΔ01Plrdsdt,which can be written as(81)∂DP∂xjka=∫ℛlogq1−q2m2q2s∑l=0m−1∑r=0nQj−1,qm−1sQl,qm−1s−Qj,qm−1sQl,qm−1sQk,qntQr,qntea,Δ10Plr+n2q2t∑l=0m∑r=0n−1Qj,qmsQl,qn−1sQk−1,qn−1tQr,qn−1t−Qk,qn−1tQr,qn−1tea,Δ01Plrdsdt.
The equation (81) can be written as(82)∂DP∂xjka=Rjk+Sjk,where(83)Rjk=mlogq1−q2∫ℛMjks,tdsdt,and(84)Sjk=nlogq1−q2∫ℛNjks,tdsdt,in which Mjks,t represents(85)Mjks,t=q2s∑l=0m−1∑r=0nQj−1,qm−1sQl,qm−1s−Qj,qm−1sQl,qm−1sQk,qntQr,qntea,Δ10Plr,and Njks,t stands for(86)Njks,t=q2t∑l=0m∑r=0n−1Qk−1,qn−1tQr,qn−1t−Qk,qn−1tqQr,qn−1tQj,qmsQl,qmsea,Δ01Plr.
Equations (85) and (86) can be written as(87)Mjks,t=q2s∑l=0m−1∑r=0nm−1j−1m−1l2m−2j+l−1Qj+l,q2m−2s−m−1jm−1l2m−2l+jQl+j,q2m−2tnknr2nk+rQk+r,q2ntea,Δ10Plr,and(88)Njks,t=q2t∑l=0m∑r=0n−1n−1k−1n−1r2n−2k+r−1Qk+r−1,q2n−2t−n−1kn−1r2n−2r+kQr+k,q2n−2tmlmj2ml+jQl+j,q2msea,Δ01Plr.
The equations (87) and (88) can be rewritten as(89)Mjks,t=q2s∑l=0m−1∑r=0nm−1j−1m−1l2m−2j+l−1Qj+l,q2m−2s−m−1jm−1l2m−2l+jQl+j,q2m−2tnknr2nk+rQk+r,q2ntΔ10Plr,and(90)Njks,t=q2t∑l=0m∑r=0n−1n−1rn−1k2n−2k+rQk+r,q2n−2t−n−1k−1n−1r2n−2r+k−1Qr+k−1,q2n−2tmlmj2ml+jQl+j,q2msΔ01Plr.
Now, we find the integrals ofMjks,t (89) and that of Njks,t (90),(91)∫RMjks,tdsdt=∑l=0m−1∑r=0nξj−1,lm−1,m−1δj+l−12m−2s−ξj,lm−1,m−1δj+l2m−2sξk,rn,nδk+r2ntΔ10Plr,and(92)∫RNjks,tdsdt=∑l=0m∑r=0n−1ξk−1,rn−1,n−1δk+r−12n−2t−ξr,kn−1,n−1δr+k2n−2sξl,jm,mδl+j2msΔ01Plr,where(93)ξj,lm,m=mjml2mj+l,δj+l2ms=∫Rq2sQj+l2msds.
Plugging the (83) and (84) in (91) and (92), respectively, it gives(94)Rjk=mlogq1−q2∑i,j=0m−1,nξj−1,lm−1,m−1δj+l−12m−2s−ξj,lm−1,m−1δj+l2m−2sξk,rn,nδk+r2ntΔ10Plr,and(95)Sjk=nlogq1−q2∑i,j=0m,n−1ξk−1,rn−1,n−1δk+r−12n−2t−ξr,kn−1,n−1δr+k2n−2tξl,jm,mδl+j2msΔ01Plr.
Substituting the value ofRjk (94) and the value of Sjk (95) in (81), we find that(96)∂DP∂xjka=mlogq1−q2∑i,j=0m−1,nξj−1,lm−1,m−1δj+l−12m−2s−ξj,lm−1,m−1δj+l2m−2sξk,rn,nδk+r2ntΔ10Pij+nlogq1−q2∑i,j=0m,n−1ξk−1,rn−1,n−1δk+r−12n−2t−ξr,kn−1,n−1δr+k2n−2tξl,jm,mδl+j2msΔ01Plr.
The vanishing condition of gradient of Dirichlet function is∂DP/∂xjka=0 for the q- Bernstein–Bézier surface. Thus, the above (96) yields the equation (43).Corollary 7.
A control net,P=Pjkj,k=0n,m, is an extremal of the Dirichlet functional with the prescribed border for m=n if the following constraint on the interior control points is satisfied.(97)mlogq1−q2∑i,j=0m−1,mξj−1,lm−1,m−1δj+l−12m−2s−ξj,lm−1,m−1δj+l2m−2sξk,rm,mδk+r2mt×Δ10Plr+∑i,j=0m,m−1ξk−1,rm−1,m−1δk+r−12m−2t−ξr,km−1,m−1δr+k2m−2tξl,jm,mδl+j2msΔ01Plr=0.Corollary 8.
A bi-quadratic quasi-minimalq-Bernstein–Bézier surface as the extremal of Dirichlet functional satisfies the following constraint, which explicitly gives P11, the only interior point in this case as the linear combination of known boundary control points (by plugging n=m=2 and j=k=1 in the equation (97)).(98)q2tΔ01P00+Δ01P01+Δ01P02+Δ01P10+Δ01P11+Δ01P12+Δ01P21+Δ01P22q3−15q2+27q−37−2q2sΔ10P00+Δ10P01+Δ10P02+Δ10P10+Δ10P11+Δ10P12+Δ10P21+Δ10P221−q=0.Corollary 9.
The shift operator properties thatΔ01Pij=Pij+1−Pij and Δ10Pij=Pi+1j−Pij, for q=0.2, enable us to reduce the equation (98) for the bi-quadratic quasi-minimal q-Bernstein–Bézier surface as the extremal of Dirichlet functional as follows:(99)P11=1113−29P00+29P01−29P02+42P10+86P12−29P20+0.1P21−2P22.Figure4 represents the prescribed control net (first figure of Figure 4), related bi-quadratic quasi-minimal q-Bernstein–Bézier surface (for q=0.2, m=2, and n=2 given by second figure of Figure 4) and the mean curvature function (third figures of Figure 4) for the schematic illustration of the quasi-minimal q-Bernstein–Bézier surfaces described in Proposition 1, together with the unknown interior point P11 obtained from the vanishing condition for the gradient of Dirichlet functional (99). Corollary 10. A bi-quadratic quasi-minimal q-Bernstein–Bézier surface as the extremal of Dirichlet functional is obtained by taking n=m=2; in this case, we have only the choice that j=1 and k=1 for Pjk, which gives us P11 in terms of known boundary control points. In particular, for q=1, the above expression (98) reduces to equation (100) which is standard result obtained for P11 for bi-quadratic quasi-minimal Bézier surface by Monterde [19], and its mean curvature function of surface parameters is shown in Figure 5.(100)P11=2P00−3P01−5P12+9P02+P10+2P02+3P21+P22.Figure 4
Prescribed border, the quasi-minimal bi-quadraticq-Bernstein–Bézier surface, and the mean curvature function.Figure 5
The mean curvature function of the bi-quadratic quasi-minimalq-Bernstein–Bézier surface.Corollary 11.
A bicubicq-Bernstein–Bézier surface is an extremal of the Dirichlet functional with prescribed border for n=m=3 and j=k=1,2. For example, for q=0.2, the constraint equation (43) gives us four constraints on the interior control points P11,P12,P21,P22 that depend on the boundary control points, which are(101)P11=−11.48P00−11.48P01−11.48P02−P10−P13−0.50P20−0.50P23+7.96P31+7.97P32+7.97P33,P12=4.86P01+4.86P02+4.86P03+0.81P13+P20+P21+P22+P23−2.16P30−2.16P31−2.16P32,P21=4.86P01+4.86P02+4.86P03+0.81P13+P20+P21+P22+P23−2.16P30−2.16P31−2.16P32,P22=4.10P00+4.32P01+4.32P02−0.33P10−0.33P13+P20+P23+1.99P30+1.99P31+1.99P32+1.99P33.As an instance of a bi-cubic quasi-minimalq-Bernstein–Bézier surface for q=0.2, m=3, and n=3, Figure 6 serves as an illustrative example of the Proposition 1 for the quasi-minimal q-Bernstein–Bézier surface. In Figure 6, the boundary control points of the bicubic quasi-minimal q-Bernstein–Bézier surface are shown in its first figure, whereas the bicubic quasi-minimal q-Bernstein–Bézier surface and its mean curvature are shown in second and third figures along with the unknown interior control points P11, P12, P21, and P22, worked out from the vanishing condition of gradient of the Dirichlet functional given by equation (43) as control net of the surface. We conclude the section with the remark that the problem of finding the quasi-minimal q-Bernstein–Bézier surface as the extremal of Dirichlet functional is thus reduced to solving the algebraic constraints for the interior control points in terms of boundary control points as the outcome of vanishing condition of gradient of Dirichlet functional. For schematic illustration, bi-quadratic and bicubic quasi-minimal q-Bernstein–Bézier surfaces are given as representative examples.Figure 6
Prescribed border, the bi-cubic quasi-minimalq-Bernstein–Bézier surface, and the mean curvature function.


## 4. Conclusion

The Bézier surfaces and minimal surfaces arising as the extremal of certain energy functional appear quite frequently in the mathematical models of surface formation in computer science for computer-aided geometric design (CAGD), computer graphics, and other disciplines of mathematics. A minimal surface is defined as a surface of minimal area that has vanishing mean curvature everywhere on the surface. One of the widely used restrictions is to find out the Bézier surface of minimal area as the extremal of various energy integrals by the vanishing condition of gradient of such a functional. A class of surfaces, namely, theq-Bernstein–Bézier surfaces, have been discussed and the corresponding quasi-minimal surfaces as the extremal of Dirichlet functional are achieved by finding the linear algebraic constraints on the interior control points of the prescribed border (boundary control points) as the outcome of the vanishing condition for this functional gradient of Dirichlet functional. The unknown interior control points computed through this scheme along with given boundary control points can be used to plot the minimal surfaces. We call the quasi-minimal surface determined as the extremal of Dirichlet functional as the quasi-minimal q-Bernstein–Bézier surface. We have included the representative examples for the illustration of the scheme for bi-quadratic and bi-cubic q-Bernstein–Bézier surfaces for q=0.2. The work can be extended not only for the surfaces as the extremal of other functionals but also for the surfaces spanned by various polynomials as well as for Stancu polynomials, Hermite polynomials, Bernoulli polynomials, and others.



---

*Source: 8994112-2022-09-20.xml*