 approach are that, first it reduces the given problems to those of solving a system of algebraic equations, and in general converges exponentially and almost always supplies the most terse representation of a smooth solution([@a13], [@a14], [@aa26]).

In this article, we use shifted GJPs on $\Omega$, which are mutually orthogonal with respect to the shifted weight function $\delta^{\alpha,\beta}(v)=(2-2v)^\alpha(2v)^\beta$ on $\Omega$ where $\alpha, \beta$ belong to one of the following index sets $$\begin{aligned}
{\mathcal N_1}&=&\{(\alpha,\beta); \alpha, \beta \le -1, ~\alpha,
\beta \in \mathbb Z\},\quad~~~~~~~~~~~~~~ {\mathcal
N_2}=\{(\alpha,\beta); \alpha \le -1, \beta > -1,~~\alpha \in
\mathbb Z, \beta \in \mathbb R\},
 \\
{\mathcal N_3}&=&\{(\alpha,\beta); \alpha>-1, \beta \le -1,~\alpha
\in \mathbb R, \beta \in \mathbb Z\},\quad {\mathcal
N_4}=\{(\alpha,\beta); \alpha, \beta
> -1,~\alpha, \beta \in \mathbb R\},\end{aligned}$$ where the symbol $\mathbb Z$ is the collection of all integer numbers. The main advantage of GJPs is that these polynomials, with indexes corresponding to the number of homogeneous initial conditions in a given FIDE, are the natural basis functions to the Galerkin approximation of this problem([@a15], [@a16]).

The organization of this paper is as follows: we begin by reviewing some preliminaries which are required for establishing our results in Section 2. In Section 3, we introduce the discrete Galerkin method based on the GJPs and its application to (\[6x\]). Numerical solvability of the algebraic system obtained from discrete Galerkin discretization of a special case of (\[6x\]) with $0<q<\frac{1}{2}$ and $\bar p(v)=1$ based on GJPs is given in Section 4. Convergence analysis of the proposed scheme is provided in Section 5. Numerical experiments are carried out in Section 6.

Preliminaries and Notations
===========================

In this section, we
approach are that, first it reduces the given trouble to those of clear a system of algebraic equations, and in general converges exponentially and about always supplies the most crisp theatrical performance of a smooth solution([@a13 ], [ @a14 ], [ @aa26 ]). 

 In this article, we use shifted GJPs on $ \Omega$, which are mutually orthogonal with respect to the shift weight function $ \delta^{\alpha,\beta}(v)=(2 - 2v)^\alpha(2v)^\beta$ on $ \Omega$ where $ \alpha, \beta$ belong to to one of the following index sets $ $ \begin{aligned } 
 { \mathcal N_1}&=&\{(\alpha,\beta); \alpha, \beta \le -1, ~\alpha, 
 \beta \in \mathbb Z\},\quad~~~~~~~~~~~~~~ { \mathcal 
 N_2}=\{(\alpha,\beta); \alpha \le -1, \beta > -1,~~\alpha \in 
 \mathbb Z, \beta \in \mathbb R\ }, 
  \\ 
 { \mathcal N_3}&=&\{(\alpha,\beta); \alpha>-1, \beta \le -1,~\alpha 
 \in \mathbb R, \beta \in \mathbb Z\},\quad { \mathcal 
 N_4}=\{(\alpha,\beta); \alpha, \beta 
 > -1,~\alpha, \beta \in \mathbb R\},\end{aligned}$$ where the symbol $ \mathbb Z$ is the collection of all integer numbers. The main advantage of GJPs is that these polynomials, with index corresponding to the number of homogeneous initial conditions in a give FIDE, are the natural basis functions to the Galerkin approximation of this problem([@a15 ], [ @a16 ]). 

 The arrangement of this paper is as follows: we begin by reviewing some preliminaries which are command for establishing our results in Section 2. In Section 3, we introduce the discrete Galerkin method based on the GJPs and its application to (\[6x\ ]). Numerical solvability of the algebraic system obtained from discrete Galerkin discretization of a special case of (\[6x\ ]) with $ 0 < q<\frac{1}{2}$ and $ \bar p(v)=1 $ based on GJPs is given in Section 4. Convergence psychoanalysis of the propose scheme is provide in Section 5. Numerical experiments are carried out in incision 6. 

 Preliminaries and Notations 
 = = = = = = = = = = = = = = = = = = = = = = = = = = = 

 In this section, we
 aporoach are that, first it reduces the girwn proulems tk those uf solving a system of algebcaic equaupons, and in general cunverges vxponentiqlly qnd almost always supplisd thz nost terse reptesentation mf a smooth sonugiln([@a13], [@a14], [@aa26]).

In this article, we use shifeed GJPx ln $\Omega$, which are iutuzlly orthogonal with respect to ths shifttd weight function $\delta^{\alpha,\beta}(v)=(2-2v)^\alpha(2v)^\betw$ on $\Omega$ where $\alpha, \beta$ belont to ine of the fullowing index sets $$\befin{aligned}
{\mathcal N_1}&=&\{(\alpha,\beta); \aloha, \bzta \le -1, ~\alpya,
\vetw \in \mathbb V\},\quad~~~~~~~~~~~~~~ {\iathcal
N_2}=\{(\alpha,\beta); \al[ha \le -1, \beta > -1,~~\alpha \ik
\mathub Z, \beta \in \mathbb R\},
 \\
{\matical N_3}&=&\{(\alpha,\beta); \alphw>-1, \beta \le -1,~\ampha
\in \mathbb R, \bwtq \in \kathtb Z\},\dyad {\mauhcel
N_4}=\{(\zlpha,\bfta); \alpha, \betz
> -1,~\alpha, \betq \in \mathbb R\},\end{alibnqe}$$ where the sgmbol $\iaehbb Z$ is the collection of all integer nujbers. The main advantagw of GJPs is that thede polynoiials, with indexes corresponding to the number of homojeveobw iniguap conditions in a given FIDE, are the natural fzsos functions to bhe Galerkin approcilayyon of this ptoblem([@a15], [@a16]).

Tge organization of this pwper us as foljows: we begin by reviewing some preliminarivs wyich are required yor establisking oor resilts in Section 2. In Seccion 3, se introducf the disdfete Galerkin meghoc tased on the GJPs and its wpplicatiin tp (\[6x\]). Numdricsl soldability ov the algebraic system ohtaingd frok discrete Galerkin discretization of a s'xcial case of (\[6x\]) winh $0<q<\frac{1}{2}$ cnd $\bav p(v)=1$ based on GTPs is given iu Sectiou 4. Conxergence ahalysis of the pro[osed scheme hd provided ii Section 5. Nunerixal expdfiments are catried out in Sectiob 6.

Preliminaries anv Nojafions
===========================

In this seeuiob, we
approach are that, first it reduces the to of solving system of algebraic exponentially almost always supplies most terse representation a smooth solution([@a13], [@a14], [@aa26]). In article, we use shifted GJPs on $\Omega$, which are mutually orthogonal with respect the shifted weight function $\delta^{\alpha,\beta}(v)=(2-2v)^\alpha(2v)^\beta$ on $\Omega$ where $\alpha, \beta$ belong to one the index $$\begin{aligned} N_1}&=&\{(\alpha,\beta); \alpha, \beta \le -1, ~\alpha, \beta \in \mathbb Z\},\quad~~~~~~~~~~~~~~ {\mathcal N_2}=\{(\alpha,\beta); \alpha \le -1, \beta -1,~~\alpha \in \mathbb Z, \beta \in \mathbb R\}, {\mathcal N_3}&=&\{(\alpha,\beta); \alpha>-1, \beta -1,~\alpha \in \mathbb R, \beta \mathbb {\mathcal N_4}=\{(\alpha,\beta); \beta -1,~\alpha, \in \mathbb R\},\end{aligned}$$ the symbol $\mathbb Z$ is the collection of all integer numbers. The main advantage of GJPs is these polynomials, corresponding to number homogeneous conditions in a are the natural basis functions to of this problem([@a15], [@a16]). The organization of this is as we begin by reviewing some preliminaries are required for establishing our results in Section In Section 3, we introduce the discrete Galerkin method based on the GJPs and its (\[6x\]). Numerical solvability of algebraic system obtained discrete discretization a case of with $0<q<\frac{1}{2}$ and $\bar p(v)=1$ based on GJPs is given in 4. Convergence analysis of the proposed scheme is provided in Numerical are carried out Section 6. Preliminaries and =========================== this section, we
 approach are that, first it redUces the givEn proBleMs tO tHose Of soLving a system of ALgebRaic equations, and in geneRal coNvERges EXpOnentIally anD AlMOSt aLwAyS suPpLIeS the mOst Terse rePresentatiOn oF a Smooth solutiON([@a13], [@A14], [@aa26]).

In this aRtiCle, we use shifTed gJPs on $\omEga$, WHich aRe mUtualLy orthOGonal wIth respecT tO The shiFTed weigHT FuNctiOn $\delta^{\alpha,\beta}(v)=(2-2V)^\AlPHa(2v)^\beta$ on $\Omega$ Where $\aLpHA, \bETA$ beLonG to one of thE fOllowINg index SEtS $$\BEGin{ALigned}
{\mathcal n_1}&=&\{(\alpha,\beta); \aLPha, \Beta \le -1, ~\AlPha,
\BEta \in \mAthbb z\},\qUAd~~~~~~~~~~~~~~ {\mAthcal
N_2}=\{(\alphA,\betA); \alpha \le -1, \bEta > -1,~~\alpHA \in
\mathBB Z, \beta \iN \mathbB R\},
 \\
{\mAthCal N_3}&=&\{(\ALpHa,\BetA); \aLPha>-1, \BEtA \le -1,~\ALphA
\in \mathbB R, \BeTa \in \mAthbB z\},\QUAd {\maThcAl
N_4}=\{(\aLpha,\bEta); \alpha, \beta
> -1,~\aLphA, \betA \In \mAthbb r\},\end{aLignEd}$$ Where The symBol $\maThBb Z$ is the collectIon oF all integEr nUmBerS. THe maiN AdvantAge Of GjPs is thAt these POlyNoMIALs, With indexes correspOnDINg To the numBer of hOMoGeNEous initIaL coNditIONs in a GiveN fIdE, are the NaturaL BaSiS functiOnS to the gaLerKin ApproXImatIon of tHis problEm([@a15], [@a16]).

THE organization oF This paper is as FOlLOWs: WE begIn bY reviewing sOme pRElimInarIEs WhiCH are rEquirEd FOr EStablishing our resulTs In SectIon 2. In section 3, we intrOduce the diSCREte GalerKin mEThOD based on the GJPS and iTs applicatIOn to (\[6x\]). NumEricaL solvabiLity of the ALGebraic sYstEm oBtaIneD FRoM discrete GaleRKIn diScRetizatIon Of a specIal CasE of (\[6X\]) wiTh $0<Q<\frac{1}{2}$ and $\bAr p(v)=1$ baseD oN GjPS iS giVen in sEction 4. CoNvErgEnCe aNalysIS of the PropoSed sChEmE Is pRovided IN SECTion 5. nuMeRicaL exPeRimenTs arE CarRied out In Section 6.

preLIminArIeS and NotAtions
===========================

In this sEcTion, we
 approach are that, firstit reduces  thegiv enpr oble ms t o those of sol v inga system of algebraicequat io n s, a n din ge neral c o nv e r ges  e xp one nt i al ly an d a lmost a lways supp lie sthe most ter s erepresenta tio n of a smoot h s olutio n( [@a 1 3], [ @a1 4], [ @aa26] ) .

Inthis arti cl e , we u s e shift e d  G JPson $\Omega$, whic h  a r e mutually ort hogona lw it h res pec t to the s hi ftedw eight f u nc t i o n $ \ delta^{\alpha ,\beta}(v)= ( 2-2 v)^\al ph a(2 v )^\bet a$ on  $ \ Ome ga$ where $ \alp ha, \beta $ belo n g to on e  of the  follo win g i ndex se ts  $$ \b e gin { al ign e d}{\mathca lN_ 1}&=& \{(\ a l p h a,\b eta ); \ alpha , \beta \le - 1,~\al p ha, 
\bet a \in  \ma th bb Z\ },\qua d~~~~ ~~ ~~~~~~~~ {\math calN_2}=\{(\ alp ha ,\b et a); \ a lpha \ le-1,  \beta> -1,~~ \ alp ha \ i n\mathbb Z, \beta \ in \ ma thbb R\} ,
 \\{ \m at h cal N_3} &= &\{ (\al p h a,\be ta); \a lpha>-1,  \beta \l e-1,~\al ph a
\in\m ath bbR, \b e ta \ in \ma thbb Z\} ,\qua d  {\mathcal
N_4 } =\{(\alpha,\b e ta ) ;  \ a lpha , \ beta
> -1,~ \alp h a, \ beta \i n \ m athbb  R\}, \e n d{ a ligned}$$ where the  s ymbol$\mat hbb Z$ is the  collectio n o f all in tege r  n u mbers. The mai n adv antage ofG JPs is t hat t hese pol ynomials, w ith inde xes  co rre spo n d in g to the numb e r  ofho mogeneo usinitial  co ndi tio nsin  a givenFIDE, ar eth ena tur al ba s is funct io nsto  th e Gal e rkin a pprox imat io no f t his pro b le m ( [@a1 5] ,[@a1 6]) .
Theorga n iza tion of  this pap eri s as  f ol lows: w e begin by re vi ewing some  p rel iminar i e s whichare required for establ i shing o urresul ts i n Section  2.  In Se cti o n 3, w e intr oduce  t hed i scret e Ga ler ki n method b a s edon th eGJPs  and it s application to ( \ [6x \]). Numerica l s olva b i li tyo ft heal g ebr a i c system obtain ed from di sc r et e Galerkin dis cr etizati on of a  spec i al case  of (\[6x \]) with$0 <q<\ f r ac{ 1}{2}$ and  $\bar p (v)=1$ ba s ed on GJ Ps is  gi ven in  S ect ion 4 . Conv e rge nce a nalysi sof the  prop os ed schem e is provided in Sectio n 5. N umeri cal  experime nts are  carriedoutin Section  6. 

P relim ina r ies a nd N o ta tio n s
=== ==== = ========= = == === = = == 

In this s e c t ion , we
 approach_are that,_first it reduces the_given problems_to_those of_solving_a system of_algebraic equations, and_in general converges exponentially_and almost always_supplies_the most terse representation of a smooth solution([@a13], [@a14], [@aa26]).

In this article, we use_shifted_GJPs on_$\Omega$,_which_are mutually orthogonal with respect_to the shifted weight function_$\delta^{\alpha,\beta}(v)=(2-2v)^\alpha(2v)^\beta$ on_$\Omega$ where $\alpha, \beta$ belong to one of_the_following index sets_$$\begin{aligned}
{\mathcal N_1}&=&\{(\alpha,\beta); \alpha, \beta \le -1, ~\alpha,
\beta \in \mathbb_Z\},\quad~~~~~~~~~~~~~~ {\mathcal
N_2}=\{(\alpha,\beta); \alpha \le -1, \beta_> -1,~~\alpha \in
\mathbb_Z,_\beta_\in \mathbb R\},
 \\
{\mathcal_N_3}&=&\{(\alpha,\beta); \alpha>-1, \beta \le -1,~\alpha
\in \mathbb_R, \beta \in \mathbb Z\},\quad {\mathcal
N_4}=\{(\alpha,\beta);_\alpha, \beta
> -1,~\alpha, \beta \in \mathbb R\},\end{aligned}$$_where the symbol $\mathbb Z$ is_the collection of all integer_numbers. The_main advantage of GJPs is_that these polynomials,_with indexes_corresponding to the_number of homogeneous initial conditions in_a given FIDE,_are the natural basis functions to_the_Galerkin approximation of_this_problem([@a15],_[@a16]).

The organization_of this paper_is_as follows:_we_begin by reviewing some preliminaries which_are_required for establishing our results in Section_2. In Section 3,_we_introduce the discrete Galerkin_method based on the GJPs_and its application to (\[6x\]). Numerical_solvability of_the algebraic_system obtained from discrete Galerkin discretization of a special case of_(\[6x\]) with $0<q<\frac{1}{2}$ and $\bar p(v)=1$_based on GJPs is_given in_Section_4. Convergence analysis_of_the proposed_scheme is provided in Section 5. Numerical_experiments are_carried out in Section 6.

Preliminaries and_Notations
===========================

In this section, we