\documentclass[11pt,twoside]{article}\makeatletter \IfFileExists{xcolor.sty}% {\RequirePackage{xcolor}}% {\RequirePackage{color}} \usepackage{colortbl} \usepackage{wrapfig} \usepackage{ifxetex} \ifxetex \usepackage{fontspec} \usepackage{xunicode} \catcode`⃥=\active \def⃥{\textbackslash} \catcode`❴=\active \def❴{\{} \catcode`❵=\active \def❵{\}} \def\textJapanese{\fontspec{Noto Sans CJK JP}} \def\textChinese{\fontspec{Noto Sans CJK SC}} \def\textKorean{\fontspec{Noto Sans CJK KR}} \setmonofont{DejaVu Sans Mono} \else \IfFileExists{utf8x.def}% {\usepackage[utf8x]{inputenc} \PrerenderUnicode{–} }% {\usepackage[utf8]{inputenc}} \usepackage[english]{babel} \usepackage[T1]{fontenc} \usepackage{float} \usepackage[]{ucs} \uc@dclc{8421}{default}{\textbackslash } \uc@dclc{10100}{default}{\{} \uc@dclc{10101}{default}{\}} \uc@dclc{8491}{default}{\AA{}} \uc@dclc{8239}{default}{\,} \uc@dclc{20154}{default}{ } \uc@dclc{10148}{default}{>} \def\textschwa{\rotatebox{-90}{e}} \def\textJapanese{} \def\textChinese{} \IfFileExists{tipa.sty}{\usepackage{tipa}}{} \fi \def\exampleFont{\ttfamily\small} \DeclareTextSymbol{\textpi}{OML}{25} \usepackage{relsize} \RequirePackage{array} \def\@testpach{\@chclass \ifnum \@lastchclass=6 \@ne \@chnum \@ne \else \ifnum \@lastchclass=7 5 \else \ifnum \@lastchclass=8 \tw@ \else \ifnum \@lastchclass=9 \thr@@ \else \z@ \ifnum \@lastchclass = 10 \else \edef\@nextchar{\expandafter\string\@nextchar}% \@chnum \if \@nextchar c\z@ \else \if \@nextchar l\@ne \else \if \@nextchar r\tw@ \else \z@ \@chclass \if\@nextchar |\@ne \else \if \@nextchar !6 \else \if \@nextchar @7 \else \if \@nextchar (8 \else \if \@nextchar )9 \else 10 \@chnum \if \@nextchar m\thr@@\else \if \@nextchar p4 \else \if \@nextchar b5 \else \z@ \@chclass \z@ \@preamerr \z@ \fi \fi \fi \fi \fi \fi \fi \fi \fi \fi \fi \fi \fi \fi \fi \fi} \gdef\arraybackslash{\let\\=\@arraycr} \def\@textsubscript#1{{\m@th\ensuremath{_{\mbox{\fontsize\sf@size\z@#1}}}}} \def\Panel#1#2#3#4{\multicolumn{#3}{){\columncolor{#2}}#4}{#1}} \def\abbr{} \def\corr{} \def\expan{} \def\gap{} \def\orig{} \def\reg{} \def\ref{} \def\sic{} \def\persName{}\def\name{} \def\placeName{} \def\orgName{} \def\textcal#1{{\fontspec{Lucida Calligraphy}#1}} \def\textgothic#1{{\fontspec{Lucida Blackletter}#1}} \def\textlarge#1{{\large #1}} \def\textoverbar#1{\ensuremath{\overline{#1}}} \def\textquoted#1{‘#1’} \def\textsmall#1{{\small #1}} \def\textsubscript#1{\@textsubscript{\selectfont#1}} \def\textxi{\ensuremath{\xi}} \def\titlem{\itshape} \newenvironment{biblfree}{}{\ifvmode\par\fi } \newenvironment{bibl}{}{} \newenvironment{byline}{\vskip6pt\itshape\fontsize{16pt}{18pt}\selectfont}{\par } \newenvironment{citbibl}{}{\ifvmode\par\fi } \newenvironment{docAuthor}{\ifvmode\vskip4pt\fontsize{16pt}{18pt}\selectfont\fi\itshape}{\ifvmode\par\fi } \newenvironment{docDate}{}{\ifvmode\par\fi } \newenvironment{docImprint}{\vskip 6pt}{\ifvmode\par\fi } \newenvironment{docTitle}{\vskip6pt\bfseries\fontsize{22pt}{25pt}\selectfont}{\par } \newenvironment{msHead}{\vskip 6pt}{\par} \newenvironment{msItem}{\vskip 6pt}{\par} \newenvironment{rubric}{}{} \newenvironment{titlePart}{}{\par } \newcolumntype{L}[1]{){\raggedright\arraybackslash}p{#1}} \newcolumntype{C}[1]{){\centering\arraybackslash}p{#1}} \newcolumntype{R}[1]{){\raggedleft\arraybackslash}p{#1}} \newcolumntype{P}[1]{){\arraybackslash}p{#1}} \newcolumntype{B}[1]{){\arraybackslash}b{#1}} \newcolumntype{M}[1]{){\arraybackslash}m{#1}} \definecolor{label}{gray}{0.75} \def\unusedattribute#1{\sout{\textcolor{label}{#1}}} \DeclareRobustCommand*{\xref}{\hyper@normalise\xref@} \def\xref@#1#2{\hyper@linkurl{#2}{#1}} \begingroup \catcode`\_=\active \gdef_#1{\ensuremath{\sb{\mathrm{#1}}}} \endgroup \mathcode`\_=\string"8000 \catcode`\_=12\relax \usepackage[a4paper,twoside,lmargin=1in,rmargin=1in,tmargin=1in,bmargin=1in,marginparwidth=0.75in]{geometry} \usepackage{framed} \definecolor{shadecolor}{gray}{0.95} \usepackage{longtable} \usepackage[normalem]{ulem} \usepackage{fancyvrb} \usepackage{fancyhdr} \usepackage{graphicx} \usepackage{marginnote} \renewcommand{\@cite}[1]{#1} \renewcommand*{\marginfont}{\itshape\footnotesize} \def\Gin@extensions{.pdf,.png,.jpg,.mps,.tif} \pagestyle{fancy} \usepackage[pdftitle={Fourier Spectral Methods for Numerical Solving of the Kuramoto-Sivashinsky Equation}, pdfauthor={}]{hyperref} \hyperbaseurl{} \paperwidth210mm \paperheight297mm \def\@pnumwidth{1.55em} \def\@tocrmarg {2.55em} \def\@dotsep{4.5} \setcounter{tocdepth}{3} \clubpenalty=8000 \emergencystretch 3em \hbadness=4000 \hyphenpenalty=400 \pretolerance=750 \tolerance=2000 \vbadness=4000 \widowpenalty=10000 \renewcommand\section{\@startsection {section}{1}{\z@}% {-1.75ex \@plus -0.5ex \@minus -.2ex}% {0.5ex \@plus .2ex}% {\reset@font\Large\bfseries}} \renewcommand\subsection{\@startsection{subsection}{2}{\z@}% {-1.75ex\@plus -0.5ex \@minus- .2ex}% {0.5ex \@plus .2ex}% {\reset@font\Large}} \renewcommand\subsubsection{\@startsection{subsubsection}{3}{\z@}% {-1.5ex\@plus -0.35ex \@minus -.2ex}% {0.5ex \@plus .2ex}% {\reset@font\large}} \renewcommand\paragraph{\@startsection{paragraph}{4}{\z@}% {-1ex \@plus-0.35ex \@minus -0.2ex}% {0.5ex \@plus .2ex}% {\reset@font\normalsize}} \renewcommand\subparagraph{\@startsection{subparagraph}{5}{\parindent}% {1.5ex \@plus1ex \@minus .2ex}% {-1em}% {\reset@font\normalsize\bfseries}} \def\l@section#1#2{\addpenalty{\@secpenalty} \addvspace{1.0em plus 1pt} \@tempdima 1.5em \begingroup \parindent \z@ \rightskip \@pnumwidth \parfillskip -\@pnumwidth \bfseries \leavevmode #1\hfil \hbox to\@pnumwidth{\hss #2}\par \endgroup} \def\l@subsection{\@dottedtocline{2}{1.5em}{2.3em}} \def\l@subsubsection{\@dottedtocline{3}{3.8em}{3.2em}} \def\l@paragraph{\@dottedtocline{4}{7.0em}{4.1em}} \def\l@subparagraph{\@dottedtocline{5}{10em}{5em}} \@ifundefined{c@section}{\newcounter{section}}{} \@ifundefined{c@chapter}{\newcounter{chapter}}{} \newif\if@mainmatter \@mainmattertrue \def\chaptername{Chapter} \def\frontmatter{% \pagenumbering{roman} \def\thechapter{\@roman\c@chapter} \def\theHchapter{\roman{chapter}} \def\thesection{\@roman\c@section} \def\theHsection{\roman{section}} \def\@chapapp{}% } \def\mainmatter{% \cleardoublepage \def\thechapter{\@arabic\c@chapter} \setcounter{chapter}{0} \setcounter{section}{0} \pagenumbering{arabic} \setcounter{secnumdepth}{6} \def\@chapapp{\chaptername}% \def\theHchapter{\arabic{chapter}} \def\thesection{\@arabic\c@section} \def\theHsection{\arabic{section}} } \def\backmatter{% \cleardoublepage \setcounter{chapter}{0} \setcounter{section}{0} \setcounter{secnumdepth}{2} \def\@chapapp{\appendixname}% \def\thechapter{\@Alph\c@chapter} \def\theHchapter{\Alph{chapter}} \appendix } \newenvironment{bibitemlist}[1]{% \list{\@biblabel{\@arabic\c@enumiv}}% {\settowidth\labelwidth{\@biblabel{#1}}% \leftmargin\labelwidth \advance\leftmargin\labelsep \@openbib@code \usecounter{enumiv}% \let\p@enumiv\@empty \renewcommand\theenumiv{\@arabic\c@enumiv}% }% \sloppy \clubpenalty4000 \@clubpenalty \clubpenalty \widowpenalty4000% \sfcode`\.\@m}% {\def\@noitemerr {\@latex@warning{Empty `bibitemlist' environment}}% \endlist} \def\tableofcontents{\section*{\contentsname}\@starttoc{toc}} \parskip0pt \parindent1em \def\Panel#1#2#3#4{\multicolumn{#3}{){\columncolor{#2}}#4}{#1}} \newenvironment{reflist}{% \begin{raggedright}\begin{list}{} {% \setlength{\topsep}{0pt}% \setlength{\rightmargin}{0.25in}% \setlength{\itemsep}{0pt}% \setlength{\itemindent}{0pt}% \setlength{\parskip}{0pt}% \setlength{\parsep}{2pt}% \def\makelabel##1{\itshape ##1}}% } {\end{list}\end{raggedright}} \newenvironment{sansreflist}{% \begin{raggedright}\begin{list}{} {% \setlength{\topsep}{0pt}% \setlength{\rightmargin}{0.25in}% \setlength{\itemindent}{0pt}% \setlength{\parskip}{0pt}% \setlength{\itemsep}{0pt}% \setlength{\parsep}{2pt}% \def\makelabel##1{\upshape ##1}}% } {\end{list}\end{raggedright}} \newenvironment{specHead}[2]% {\vspace{20pt}\hrule\vspace{10pt}% \phantomsection\label{#1}\markright{#2}% \pdfbookmark[2]{#2}{#1}% \hspace{-0.75in}{\bfseries\fontsize{16pt}{18pt}\selectfont#2}% }{} \def\TheFullDate{2014-01-15 (revised: 010106, 230116 15 January 2014)} \def\TheID{\makeatother } \def\TheDate{2014-01-15} \title{Fourier Spectral Methods for Numerical Solving of the Kuramoto-Sivashinsky Equation} \author{}\makeatletter \makeatletter \newcommand*{\cleartoleftpage}{% \clearpage \if@twoside \ifodd\c@page \hbox{}\newpage \if@twocolumn \hbox{}\newpage \fi \fi \fi } \makeatother \makeatletter \thispagestyle{empty} \markright{\@title}\markboth{\@title}{\@author} \renewcommand\small{\@setfontsize\small{9pt}{11pt}\abovedisplayskip 8.5\p@ plus3\p@ minus4\p@ \belowdisplayskip \abovedisplayskip \abovedisplayshortskip \z@ plus2\p@ \belowdisplayshortskip 4\p@ plus2\p@ minus2\p@ \def\@listi{\leftmargin\leftmargini \topsep 2\p@ plus1\p@ minus1\p@ \parsep 2\p@ plus\p@ minus\p@ \itemsep 1pt} } \makeatother \fvset{frame=single,numberblanklines=false,xleftmargin=5mm,xrightmargin=5mm} \fancyhf{} \setlength{\headheight}{14pt} \fancyhead[LE]{\bfseries\leftmark} \fancyhead[RO]{\bfseries\rightmark} \fancyfoot[RO]{} \fancyfoot[CO]{\thepage} \fancyfoot[LO]{\TheID} \fancyfoot[LE]{} \fancyfoot[CE]{\thepage} \fancyfoot[RE]{\TheID} \hypersetup{citebordercolor=0.75 0.75 0.75,linkbordercolor=0.75 0.75 0.75,urlbordercolor=0.75 0.75 0.75,bookmarksnumbered=true} \fancypagestyle{plain}{\fancyhead{}\renewcommand{\headrulewidth}{0pt}} \date{} \usepackage{authblk} \providecommand{\keywords}[1] { \footnotesize \textbf{\textit{Index terms---}} #1 } \usepackage{graphicx,xcolor} \definecolor{GJBlue}{HTML}{273B81} \definecolor{GJLightBlue}{HTML}{0A9DD9} \definecolor{GJMediumGrey}{HTML}{6D6E70} \definecolor{GJLightGrey}{HTML}{929497} \renewenvironment{abstract}{% \setlength{\parindent}{0pt}\raggedright \textcolor{GJMediumGrey}{\rule{\textwidth}{2pt}} \vskip16pt \textcolor{GJBlue}{\large\bfseries\abstractname\space} }{% \vskip8pt \textcolor{GJMediumGrey}{\rule{\textwidth}{2pt}} \vskip16pt } \usepackage[absolute,overlay]{textpos} \makeatother \usepackage{lineno} \linenumbers \begin{document} \author[1]{Gentian Zavalani} \affil[1]{ Polytechnic University of Tirana/ Albania} \renewcommand\Authands{ and } \date{\small \em Received: 14 December 2013 Accepted: 1 January 2014 Published: 15 January 2014} \maketitle \begin{abstract} In this paper I presented a numerical technique for solving Kuramoto-Sivashinsky equation, based on spectral Fourier methods. This equation describes reaction diffusion problems, and the dynamics of viscous-fuid films flowing along walls. After we wrote the equation in Furie space, we get a system. In this case, the exponential time differencing methods integrate the system much more accurately than other methods since the exponential time differencing methods assume in their derivation that the solution varies slowly in time. When evaluating the coefficients of the exponential time differencing and the exponential time differencing Runge Kutta methods via the?Cauchy integral?. All computational work is done with Matlab package. \end{abstract} \keywords{discrete fourier transform, exponential time differencing, exponential time differencing runge kutta methods, cauchy integral, kuramoto-sivashinsky eq} \begin{textblock*}{18cm}(1cm,1cm) % {block width} (coords) \textcolor{GJBlue}{\LARGE Global Journals \LaTeX\ JournalKaleidoscope\texttrademark} \end{textblock*} \begin{textblock*}{18cm}(1.4cm,1.5cm) % {block width} (coords) \textcolor{GJBlue}{\footnotesize \\ Artificial Intelligence formulated this projection for compatibility purposes from the original article published at Global Journals. However, this technology is currently in beta. \emph{Therefore, kindly ignore odd layouts, missed formulae, text, tables, or figures.}} \end{textblock*} \let\tabcellsep& \section[{Introduction}]{Introduction}\par ourier analysis occurs in the modeling of timedependent phenomena that are exactly or approximately periodic. Examples of this include the digital processing of information such as speech; the analysis of natural phenomen such as earthquakes; in the study of vibrations of spherical, circular or rectangular structures; and in the processing of images. In a typical case, Fourier spectral methods write the solution to the partial differential equation as its Fourier series. Fourier series decomposes a periodic realvalued function of real argument into a sum of simple oscillating trigonometric functions (??????????, ??????????????), that can be recombined to obtain the original function. Substituting this series into the partial differential equation gives a system of ordinary differential equations for the time-dependent coefficients of the trigonometric terms in the series then we choose a timestepping method to solve those ordinary differential equations II. \section[{Fourier Series}]{Fourier Series}\par The Fourier series of a smooth and periodic real-valued function ð??"ð??"(??) ? [0; 2??] with period 2?? is L x n x f L a L n ) sin( ) ( 1 2 0 ? ? = .... 2 , 1 , 0 = n (2) dx L x n x f L b L n ) cos( ) ( 1 2 0 ? ? = .... 2 , 1 = n (3)\par Fourier series can be expressed neatly in complex form as follows ? ? = ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? + ? ? ? ? ? ? + + = 1 2 2 2 f(x) n L x in L x in n L x in L x in \section[{Discrete Fourier Transform}]{Discrete Fourier Transform}\par In many applications, particularly in analyzing of real situations, the function ( ) f x to be approximated is known only on a discrete set of "sampling points" of x. Hence, the integral (7) cannot be evaluated in a closed form and Fourier analysis cannot be applied directly. It then becomes necessary to replace continuous Fourier analysis by a discrete version of it. The linear discrete Fourier transform of a periodic (discrete) sequence of complex values ?? 0 , ?? 0 , ? , ?? ?? ?? ?1 with period ?? ?? ?? , is a sequence of periodic complex values ?? ? 0 , ?? ? 1 , ? , ?? ? ?? ?? ?1 defined by? ? = ? = 1 0 2 k 1 u \textasciitilde ? ? ? ? N j N k ij j e u N (8)\par The linear inverse transformation is? ? = = 1 0 2 k j u ? ? ? ? N k N k ij e (9)\par The most obvious application of discrete Fourier analysis consists in the numerical calculation of Fourier coefficients. Suppose we want to approximate a realvalued periodic function ( ) f x defined on the interval [0; 2??] that is sampled with an even number N ? of gridpoints jh = j x ? N 2L h = 1 N 0,1,...., j ? = ? (\textbf{10})\par by it is Fourier series. First we compute approximate values of the Fourier coefficients? ? = ? ? 1 0 2 k ) ( 1 c \textasciitilde ? ? ? ? N j N k ij e x f N (11)\par Because the discrete Fourier transform and its inverse exhibit periodicity with period ?? ?? , i.e. ?? ? ?? + ?? ?? = ?? ? ?? (this property results from the periodic nature of ?? 2???????? ?? ?? ), it makes no sense to use more than ?? ?? terms in the series, and it suffices to calculate one full period. The Fourier series formed with the approximate coefficients is? + ? = ? ? 2 / 1 2 / 2 ) ( f \textasciitilde ?? ? ? N N k N k ij k e c x (11')\par The function ð??"ð??" ?(??) not only approximates, but actually interpolates ð??"ð??"(??) at the sampling grid points ?? ??\par In matrix form, the discrete Fourier transform (8) can be written asj kj u M N ? 1 u \textasciitilde k= 1 ,... 1 , 0 , ? = ? N j k (12)\par Where? kj kj M = and ? ? ? N i e 2 ? = ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? = 1) - 1)(N - (N 1) - 2(N 1 - N 1) - 3(N 1) - 2(N 1 - N 6 4 2 3 2 1 1 1 1 1 1 1 1 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? M Similarly, the inverse discrete Fourier transform has the form k * j u ? kj M = 1 ,... 1 , 0 , ? = ? N j k (13) Where * * ? kj kj M = and ? ? ? N i e 2 * = where ð??"ð??" * is complex conjugate of ð??"ð??" ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? = 1) - 1)(N - (N * 1) - 2(N * 1 - N * 1) - 3(N * 1) - 2(N * 1 - N * 6 * 4 * 2 * 3 * 2 * * * 1 1 1 1 1 1 1 1 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? \section[{M}]{M}\par The FFT algorithm reduces the computational work required to carry out a discrete Fourier transform by reducing the number of multiplications and additions of (13) , computational time is reduced from ??(?? 2 ?? ) to ??(?? ?? log ?? ?? ).\par To apply spectral methods to a partial differential equation we need to evaluate derivatives of functions. Suppose that we have a periodic real-valued function ( ) f x ? [0; 2??] with period 2?? that is discretized I with an even number N ? of grid points, so that the grid size ? N 2L h = .The complex form of the Fourier series representation of ( )f x is /2 /2 1 f ( ) N i kx L k k N x c e ? ? ? = ? + ? ? ? ? At 2 N k ? =\par the above series gives a term , and since it cannot be differentiated, we should set its derivative to be zero at the grid points. The numerical derivatives of the function ( ) f x can be illustrated as a matrix multiplication.For the first derivative, we multiply the Fourier coefficients (11) by the corresponding differentiation matrix for an even number N ? of grid points. N N i Diag L ? ? ? ? ? ? ? ? = ? ? ? ? ? ? ? ? ? ? ? ? ? ?\par This matrix has non-zero elements only on the diagonal. For an odd number N ? of grid points the differentiation matrix corresponding to the first derivative is diagonal with elements.\par 0,1, 2,3...., 1, 0, 1 ,...., 3 2 1 2 2N N i L ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?\par Then, we compute an inverse discrete Fourier transform using (11 ? ) to return to the physical space and deduce the first derivative of f (x) on the grid. Similarly, taking the second derivative corresponds to the multiplication of the Fourier coefficients (11) by the corresponding differentiation matrix for an even number N ? of grid points. N N N i Diag L ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? = ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?\par In general, in case of an even number N ? of grid points approximating the m -th numerical derivatives of a grid function ( ) f x corresponds to the multiplication of the Fourier coefficients (11) by the corresponding differentiation matrix which is diagonal with elementsm ik L ? ? ? ? ? ? ? for 0,1, 2,3...., 1, , 1 ,...., 3 2 1 2 2 2 N N N k ? ? ? ? ? = ? ? ? ? ? ? ? ? ? ?\par with the exception that for odd derivatives we set the derivative of the highest mode \section[{Exponential Time Differencing}]{Exponential Time Differencing}\par The family of exponential time differencing schemes. This class of schemes is especially suited to semi-linear problems which can be split into a linear part which contains the stiffest part of the dynamics of the problem, and a nonlinear part, which varies more slowly than the linear part. Exponential time differencing schemes are time integration methods that can be efficiently combined with spatial approximations to provide accurate smooth solutions for stiff or highly oscillatory semi-linear partial differential equations.In this paper I will present the derivation of the explicit Exponential time differencing schemes for arbitrary order following the approach in [ ] [ ] [ ] t t = to t t t t n n ? + = = +1 to obtain. ? ? ? d t t u F t u t u t n n t c t c n n e e ? ? ? ? + + + + = 0 1 ) ), ( ( * ) ( ) ( (?? 1 )\par The next step is to derive approximations to the integral in equation (?? 1 ). This procedure does not introduce an unwanted fast time scale into the solution and the schemes can be generalized to arbitrary order.\par If we apply the Newton Backward Difference Formula, we can write a polynomial approximation to ) ), ( (? ? + + n n t t u F in the form ( ) ( ) ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? = ? = ? ? ? ? ? ? ? ? ? = ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? = ? + + ? ? ? n t G m t t u F k m m t n t G m m t t t G t t u F n k n k n m k m s m m n s m m n n n ), ( ) 1 ( * / ) 1 ( * / ) 1 ( , ) ), ( ( 0 1 0 1 0 ? ? ? ? (?? 2 ) 1 ,..., 1 ), 1 )...( 2 )( 1 ( ! ? = + ? ? ? = ? ? ? ? ? ? ? ? s m k m m m k m k note that 1 0 ! 0 = ? ? ? ? ? ? ? ? m If we substitute (?? 2 ) into (?? 1 ) we get ? ? d n t G m m t t u t u t n s m m t c t c n n e e ? ? ? ? ? ? ? ? ? ? ? ? = ? ? + ? ? ? ? ? ? ? ? ? ? ? ? + = 0 1 0 1 * / ) 1 ( * ) ( ) ( ) / ( * / * ) 1 ( ) ( ) ( 1 0 ) / 1 ( 1 0 1 t d n t G m m t t t u t u n t t c s m m t c n n e e ? ? ? ? ? ? ? ? ? ? ? ? ? ? = ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? = ? + ? ? ? (?? 3 )\par We will indicate the integral by) / ( / 1 0 ) / 1 ( t d m t e t t c m ? ? ? ? ? ? ? ? ? ? ? = ? ? ? ? ? ? ? ? (?? 4 )\par If we substitute (?? 2 ) and (?? 4 ) into (?? 3 ) we get ( )k n k n m k m s m m m t c n n t t u F k m t t u t u e ? ? = ? = ? + ? ? ? ? ? ? ? ? ? ? ? ? ? + = ), ( ) 1 ( * * ) 1 ( ) ( ) ( 0 1 0 1 ? (?? 5 )\par Which represent the general generating formula of the exponential time differencing schemes of order s The first-order exponential time differencing scheme is obtained by setting s=1c F u u n t c t c n n e e / ) 1 ( 1 ? + = ? ? +\par The second-order exponential time differencing scheme is obtained by setting s=2\{ \} ) /( ) 1 ( ) 1 2 ) 1 (( 2 1 1 t c F t c F t c t c u u n t c n t c t c n n e e e ? + ? + ? + ? ? ? + ? + = ? ? ? ? + \section[{?}]{?}\par By setting s=2 we get the fourth-order exponential time differencing scheme ) 6 /( ) ( 3 4 3 4 2 3 1 2 1 1 t c F F F F u u n n n n t c n n e ? ? ? + ? ? ? + = ? ? ? ? + 6 18 2611 24 ) 6 12 11 6 ( 2 2 3 3 2 2 3 3 1 ? ? ? ? ? ? ? + ? + ? + ? = ? ? t c? = ? + ? + ? ? ? ? ? ? ? . I V.\par On the Stability of Etdrk4 Method\par The stability analysis of the ETDRK4 method is as follows (see {\ref [21,} {\ref 24]} , , , , ? ? ? ? ? by the above expressions suffers from numerical instability for y close to zero. Because of that, for small y , instead of them, we will use their asymptotic expansions. ? Because c and ? may be complex, the stability region of the ETDRK4 method is four-dimensional and therefore quite difficult to represent. Unfortunately, we do not know any expression for ( , ) 1 r x y = we will only be able to plot it. The most common idea is to study it for each particular case; for example, assuming c to be fixed and real in {\ref [21]} or that both c and ? are pure imaginary numbers in {\ref [24]}. For dissipative PDEs with periodic boundary conditions, the scalars c that arise with a Fourier spectral method are negative. Let us take for example Burger's equation2 1 ( ) 2 t xx x u u u ? = ? [ , ] x ? ? ? ?\par where 0 1 ? < ? (3.4) Transforming it to the Fourier space gives In Figure {\ref .} 1 we draw the boundary stability regions in the complex plane x for y=0,-0.9,-5,-10,-18 .\par When the linear part is zero ( 0 y = ), we recognized the stability region of the fourth-order Runge-Kutta methods and, as y ?? ?\par , the region grows. Of course, these regions only give an indication of the stability of the method. In 0 y < , 1 y << the boundaries that are observed approach to ellipses whose parameters have been fitted numerically with the following result. x = , the intersection with the imaginary axis Im( )\par x increases as y , i.e., as 2\par ? .Since the boundary of stability grows faster than , the ETDRK4 method should have a very good behavior to solve Burger's equation, which confirms the results of paper \hyperref[b5]{[6]}.\par VI. \section[{Exponential Time Differencing Runge-Kutta Methods}]{Exponential Time Differencing Runge-Kutta Methods}\par Generally, for the one-step time-discretization methods and the Runge-Kutta methods, all the information required to start the integration is available.However, for the multi-step time-discretization methods this is not true.These methods require the evaluations of a certain number of starting values of the nonlinear term ) ), ( ( t t u F at the n -th and previous time steps to build the history required for the calculations.Therefore, it is desirable to construct exponential time differencing methods that are based on Runge-Kutta methods. \section[{Based in [ ]}]{Based in [ ]}\par 12 and [ ] = + ? ? ? + ? ? ? + ? ? + ? + + ? + + ? ? ? + ? ? + ? ? ? + ? ? (?? 12 )\par In general, the exponential time differencing Runge-Kutta method (?? 12 ) has classical order four, but Hochbruck and Ostermann {\ref [11]} showed that this method suffers from an order reduction. They also presented numerical experiments which show that the order reduction, predicted by their theory, may in fact arise in practical examples. In the worst case, this leads to an order reduction to order three for the Cox and Matthews method (?? 12 ) {\ref [12]} . However, for certain problems, such the numerical experiments conducted by {\ref Kassam and Trefethen [13]} , \hyperref[b5]{[6]} for solving various one-dimensional diffusion-type problems, and the numerical results obtained in for solving some dissipative and dispersive PDEs, the fourth-order convergence of the fourth-order Runge Kutta exponential time differencing method [12] is confirmed numerically.\par Finally, we note that as 0 c ? in the coefficients of the s -order exponential time differencing Runge-Kutta methods, the methods reduce to the corresponding order of the Runge-Kutta schemes.\par V. \section[{The Kuramoto-Sivashinsky Equation}]{The Kuramoto-Sivashinsky Equation}\par The Kuramoto-Sivashinsky equation,is one of the simplest PDEs capable of describing complex behavior in both time and space. This equation has been of mathematical interest because of its rich dynamical properties. In physical terms, this equation describes reaction diffusion problems, and the dynamics of viscousfuid films flowing along walls.\par Kuramoto-Sivashinsky equation in one space dimension can be written4 4\textbf{2 2 ) , ( ) , ( ) , ( ) , ( ) , ( x t}x u x t x u x t x u t x u t t x u ? ? ? ? ? ? ? ? ? = ? ? (14)\par Equation ( {\ref 14}) can be written in integral form if we introducet t x t x u ? ? = ) , ( ) , ( ? then 4 4 2 2 2 ) , ( ) , ( ) , (\textbf{2}1 ) , ( x t x x t x x t x t t x ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? = ? ? ? ? ? ? (15) or in form 0 2 / 2 2 4 = + + + ? ? ? u u u ut\par The Kuramoto-Sivashinsky equation with 2?? periodic boundary conditions in Fourier space can be written as follows {\ref 14}) can be written as fellows {\ref (} )? ? = = 1 0 ) ( ? ? N k L x ik k j e u x u? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? = â??" ? ? ? ? ? ? ? ? ? = ? ? â??" ? ? ? ? ? ? ? ? ? = ? ? ? ? ? ?= ? = , 1 4 , 1 2 = ? = ? ? ? ? ? ? ? ? ? ? ? ? i i ? N L h jh j x 2 , = = Equation (? k k k ik t u t t u k k 2 ) ( ) ( \textasciitilde 4 2 ? ? = ? ? (17) Where = = ? ? dx t x u L L L L x ik k e ? ? ) , (\textbf{2}1 2 [ ] 2 / 2 1 0 ) ( ) , ( 2 1 t u FFT e t x N N ijk N j j u = ? ? = ? ? ?\par In final form will be (\par )[ ] 2 4 2 ) ( 2 ) ( ) ( \textasciitilde t u FFT ik t u t t u k k k k ? ? = ? ? (18)\par Equation has strong dissipative dynamics, which arise from the fourth order dissipation The stiffness in the system (17) is due to the fact that the diagonal linear operator ( )k k 4 2 ?\par , with the elements, has some large negative real eigenvalues that represent decay, because of the strong dissipation, on a time scale much shorter than that typical of the nonlinear term. The nature of the solutions to the the Kuramoto-Sivashinsky equation varies with the system size of linear operator. For large size of linear operator, enough unstable Fourier modes exist to make the system chaotic. For small size of linear operator, insuffcient Fourier modes exist, causing the system to approach a steady state solution. In this case, the exponential time differencing methods integrate the system very much more accurately than other methods since the the exponential time differencing methods assume in their derivation that the solution varies slowly in time.\par VI. \section[{Numerical Result}]{Numerical Result}\par For the simulation tests, we choose two periodic initial conditions[ ] ? 4 , 0 , ) ( 2 cos 1 ? = ? ? ? ? ? ? x x u e x [ ] ? 4 , 0 ), sin( 4 . 2 ) cos( 6 . 0 2 sin 1 . 0 2 cos 7 . 1 ) ( 2 ? + + ? ? ? ? ? ? + ? ? ? ? ? ? = x x x x x x u\par When evaluating the coefficients of the exponential time differencing and the exponential time differencing Runge Kutta methods via the "Cauchy integral" approach , the solution appears more clear as a mesh plot and shows waves propagating, traveling periodically in time and persisting without change of shape. \section[{VII.}]{VII.} \section[{Conclusions}]{Conclusions}\par In this paper, the main objective of this study was for finding the solution of one dimensional semilinear fourth order hyperbolic Kuramoto-Sivashinsky equation, describing reaction diffusion problems, and the dynamics of viscous-fuid films flowing along walls. In order to achieve this, we applied Fourier spectral approximation for the spatial discretization. In addition, we evaluated the coeffcients of the exponential time differencing and the exponential time differencing -fourth order Runge Kutta methods via the "Cauchy integral" .Some typical examples have been demonstrated in order to illustrate the efficiency and accuracy of the exponential time differencing methods technique in this case. For the simulation tests, we chose periodic boundary conditions and applied Fourier spectral approximation for the spatial discretization. In addition, we evaluated the coefficients of the Exponential Time Differencing Runge-Kutta methods via the "Cauchy integral" approach. The equations can be used repeatedly with necessary adaptations of the initial conditions.\begin{figure}[htbp] \noindent\textbf{}\includegraphics[]{image-2.png} \caption{\label{fig_0}}\end{figure} \begin{figure}[htbp] \noindent\textbf{}\includegraphics[]{image-3.png} \caption{\label{fig_6}}\end{figure} \begin{figure}[htbp] \noindent\textbf{}\includegraphics[]{image-4.png} \caption{\label{fig_7}}\end{figure} \begin{figure}[htbp] \noindent\textbf{}\includegraphics[]{image-5.png} \caption{\label{fig_8}}\end{figure} \begin{figure}[htbp] \noindent\textbf{1362}\includegraphics[]{image-6.png} \caption{\label{fig_9}Figure 1 : 36 IFigure 2 :}\end{figure} \begin{figure}[htbp] \noindent\textbf{}\includegraphics[]{image-7.png} \caption{\label{fig_10}}\end{figure} \begin{figure}[htbp] \noindent\textbf{}\includegraphics[]{image-8.png} \caption{\label{fig_11}}\end{figure} \begin{figure}[htbp] \noindent\textbf{} \par \begin{longtable}{P{0.12323667711598746\textwidth}P{0.0006661442006269592\textwidth}P{0.020650470219435734\textwidth}P{0.0006661442006269592\textwidth}P{0.0006661442006269592\textwidth}P{0.008659874608150469\textwidth}P{0.017985893416927898\textwidth}P{0.04263322884012539\textwidth}P{0.0013322884012539184\textwidth}P{0.019984326018808778\textwidth}P{0.0006661442006269592\textwidth}P{0.006661442006269592\textwidth}P{0.0019984326018808777\textwidth}P{0.006661442006269592\textwidth}P{0.0019984326018808777\textwidth}P{0.003996865203761755\textwidth}P{0.004663009404388714\textwidth}P{0.0013322884012539184\textwidth}P{0.0006661442006269592\textwidth}P{0.006661442006269592\textwidth}P{0.0019984326018808777\textwidth}P{0.007327586206896552\textwidth}P{0.0026645768025078368\textwidth}P{0.01665360501567398\textwidth}P{0.0026645768025078368\textwidth}P{0.010658307210031347\textwidth}P{0.0019984326018808777\textwidth}P{0.0026645768025078368\textwidth}P{0.0026645768025078368\textwidth}P{0.009326018808777429\textwidth}P{0.0053291536050156735\textwidth}P{0.0006661442006269592\textwidth}P{0.0006661442006269592\textwidth}P{0.0006661442006269592\textwidth}P{0.04862852664576802\textwidth}P{0.12123824451410657\textwidth}P{0.09858934169278996\textwidth}P{0.017985893416927898\textwidth}P{0.019984326018808778\textwidth}P{0.019984326018808778\textwidth}P{0.011324451410658307\textwidth}P{0.014655172413793103\textwidth}P{0.021316614420062694\textwidth}P{0.007327586206896552\textwidth}P{0.008659874608150469\textwidth}P{0.015321316614420063\textwidth}P{0.008659874608150469\textwidth}P{0.004663009404388714\textwidth}P{0.021982758620689654\textwidth}P{0.008659874608150469\textwidth}P{0.01598746081504702\textwidth}P{0.012656739811912225\textwidth}P{0.009326018808777429\textwidth}P{0.004663009404388714\textwidth}P{0.0059952978056426335\textwidth}P{0.0026645768025078368\textwidth}P{0.0026645768025078368\textwidth}P{0.0053291536050156735\textwidth}P{0.0019984326018808777\textwidth}P{0.0019984326018808777\textwidth}} \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep ( ) du t cu t ( ) dt =\tabcellsep +\tabcellsep ?\tabcellsep ( ) u t\tabcellsep (3.2)\\ \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep and the fixed point 0 ( ) u t is stable if Re(\tabcellsep ) 0 + < . c ?\\ \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \multicolumn{4}{l}{( dt t du}\tabcellsep )\tabcellsep =\tabcellsep \multicolumn{3}{l}{cu}\tabcellsep \multicolumn{2}{l}{( t}\tabcellsep \multicolumn{2}{l}{)}\tabcellsep \multicolumn{3}{l}{+}\tabcellsep \multicolumn{2}{l}{F}\tabcellsep \multicolumn{4}{l}{( u}\tabcellsep \multicolumn{2}{l}{( t}\tabcellsep \multicolumn{3}{l}{),}\tabcellsep t\tabcellsep )\tabcellsep (3.1)\tabcellsep The application of the ETDRK4 method to (3.2) leads to a recurrence relation involving n u and 1 n u + .\\ \tabcellsep \tabcellsep \multicolumn{4}{l}{with}\tabcellsep \tabcellsep \multicolumn{28}{l}{( ( ), ) F u t t the nonlinear part, we suppose}\tabcellsep Introducing the previous notation x\tabcellsep =\tabcellsep ?\tabcellsep h\tabcellsep and y ch = ,\\ \multicolumn{35}{l}{that there exists a fixed point 0 u this means that}\tabcellsep and using the Mathematica © algebra obtain the following amplification factor\tabcellsep package, we\\ \multicolumn{35}{l}{0 cu F u t 0 ( , ) 0 + = . Linearizing about this fixed point, if ( ) u t is the perturbation of 0 u and 0 '( , ) F u t ? = then}\tabcellsep 1 + = n n r x y ( , ) u u\tabcellsep 0 = + ? ?\tabcellsep 1\tabcellsep x\tabcellsep +\tabcellsep ?\tabcellsep 2\tabcellsep x\tabcellsep 2\tabcellsep +\tabcellsep ?\tabcellsep 3\tabcellsep x\tabcellsep 3\tabcellsep +\tabcellsep ?\tabcellsep 4\tabcellsep x\tabcellsep 4\tabcellsep (3.3)\\ \multicolumn{2}{l}{where}\tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \\ \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep ?\tabcellsep 0\tabcellsep =\tabcellsep e\tabcellsep y\\ \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep y\tabcellsep 3\tabcellsep y\tabcellsep y\tabcellsep 3\tabcellsep y\\ \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \multicolumn{2}{l}{?}\tabcellsep \multicolumn{2}{l}{1}\tabcellsep \tabcellsep \multicolumn{2}{l}{=}\tabcellsep \multicolumn{12}{l}{2 3 e 4 8 3 y y ? +}\tabcellsep ?\tabcellsep 8 e y\tabcellsep 2 3\tabcellsep +\tabcellsep 4 e y\tabcellsep 2 3\tabcellsep y\tabcellsep 2 2 e y 1 4 2 y ? +\tabcellsep ?\tabcellsep 2 y e y 6\tabcellsep +\tabcellsep 4 y e\tabcellsep 2 2\tabcellsep ?\tabcellsep 2 2 y y e\\ \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \multicolumn{2}{l}{y}\tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \multicolumn{2}{l}{3}\tabcellsep \multicolumn{2}{l}{y}\tabcellsep y\tabcellsep 3\tabcellsep y\tabcellsep y\\ \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep ?\tabcellsep 2\tabcellsep \multicolumn{12}{l}{4 8 16 y y e 4 ? = +}\tabcellsep \multicolumn{2}{l}{2}\tabcellsep \multicolumn{2}{l}{?}\tabcellsep \multicolumn{6}{l}{4 e y 16}\tabcellsep \multicolumn{2}{l}{2}\tabcellsep \tabcellsep +\tabcellsep 8 e y\tabcellsep 2 4\tabcellsep y\tabcellsep 3 5 12 y y e 3 ? +\tabcellsep 2\tabcellsep ?\tabcellsep 3 e y 10\tabcellsep y\tabcellsep +\tabcellsep 4\tabcellsep e y\tabcellsep 2 3\tabcellsep ?\tabcellsep 2 3 y y e\tabcellsep ?\tabcellsep 2 2 e y 1 4 2 y +\tabcellsep ?\tabcellsep 3 3 , y e y\\ \tabcellsep \tabcellsep \multicolumn{2}{l}{?}\tabcellsep \multicolumn{2}{l}{3}\tabcellsep \multicolumn{9}{l}{5 4 16 e 5 y y = ?}\tabcellsep y 2\tabcellsep \tabcellsep \multicolumn{3}{l}{+}\tabcellsep \tabcellsep \multicolumn{4}{l}{5 e 16 y}\tabcellsep \multicolumn{2}{l}{y}\tabcellsep \tabcellsep \multicolumn{2}{l}{+}\tabcellsep \multicolumn{5}{l}{8 e y}\tabcellsep 3 2 y 5\tabcellsep ?\tabcellsep 2 5 e 20 y\tabcellsep y\tabcellsep +\tabcellsep 5 5 y e 8\tabcellsep y 2\tabcellsep ?\tabcellsep 4 2 y\tabcellsep +\tabcellsep 4 10 e y ?\tabcellsep y 2\tabcellsep +\tabcellsep 4 e y 16\tabcellsep y\tabcellsep ?\tabcellsep 3 2 y 4 e y 12\tabcellsep +\tabcellsep 6 y e\tabcellsep 2 4\tabcellsep y\tabcellsep ?\tabcellsep 2 y e\tabcellsep 5 4\tabcellsep y 2\\ \tabcellsep \tabcellsep \multicolumn{3}{l}{?}\tabcellsep \multicolumn{4}{l}{y 2 3 e 2 y}\tabcellsep \multicolumn{2}{l}{+}\tabcellsep \multicolumn{2}{l}{3 y e 4 y}\tabcellsep \multicolumn{3}{l}{?}\tabcellsep \multicolumn{4}{l}{2 e y}\tabcellsep \multicolumn{2}{l}{3 3}\tabcellsep y 2\tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \\ \tabcellsep \tabcellsep ?\tabcellsep 4\tabcellsep \tabcellsep \multicolumn{2}{l}{=}\tabcellsep \multicolumn{7}{l}{6 8 24 e 6 y y ?}\tabcellsep \multicolumn{2}{l}{y 2}\tabcellsep \multicolumn{3}{l}{+}\tabcellsep \multicolumn{5}{l}{6 e 16 y}\tabcellsep \multicolumn{2}{l}{y}\tabcellsep \tabcellsep \multicolumn{2}{l}{+}\tabcellsep \multicolumn{6}{l}{3 2 y 5 e 16 y}\tabcellsep ?\tabcellsep 2 6 e 24 y\tabcellsep y\tabcellsep +\tabcellsep 5 6 y e 8\tabcellsep y 2\tabcellsep 5 6 y + +\tabcellsep 5 18 e y ?\tabcellsep y 2\tabcellsep +\tabcellsep 5 e y 20\tabcellsep y\tabcellsep ?\tabcellsep 3 2 y 5 e y 12\tabcellsep +\tabcellsep 6 y e\tabcellsep 2 5\tabcellsep y\tabcellsep +\tabcellsep 2\tabcellsep e y\tabcellsep 5 5\tabcellsep y 2\\ \tabcellsep \tabcellsep \multicolumn{10}{l}{y 2 4 e 2 6 4 y y + ?}\tabcellsep +\tabcellsep \multicolumn{5}{l}{4 y e 6 y}\tabcellsep \tabcellsep \multicolumn{4}{l}{?}\tabcellsep \multicolumn{3}{l}{2 e y}\tabcellsep 3 4\tabcellsep \multicolumn{2}{l}{y 2}\tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \\ \tabcellsep \tabcellsep \multicolumn{2}{l}{An}\tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \multicolumn{9}{l}{important}\tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \multicolumn{10}{l}{remark:}\tabcellsep computing\\ 0\tabcellsep 1\tabcellsep 2\tabcellsep \tabcellsep \tabcellsep 3\tabcellsep \tabcellsep \tabcellsep 4\tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \tabcellsep \end{longtable} \par \caption{\label{tab_2}}\end{figure} \footnote{© 2014 Global Journals Inc. (US)} \backmatter \begin{bibitemlist}{1} \bibitem[Beylkin et al. ()]{b0}\label{b0} ‘A New Class of Time Discretization Schemes for the Solution of Nonlinear PDEs’. G Beylkin , J M Keiser , L Vozovoi . \textit{J. Comput. Phys} 1998. 147 p. . \bibitem[Norsett ()]{b3}\label{b3} ‘An A-Stable Modication of the Adams-Bashforth Methods’. S P Norsett . \textit{Conf. on Numerical Solution of Dierential Equations}, Lecture Notes in Math. (Berlin) 1969. 1969. Springer-Verlag. 109 p. . \bibitem[Burden and Faires ()]{b6}\label{b6} R L Burden , J D Faires . \textit{Numerical Analysis. Wadsworth Group}, 2001. (seventh edition) \bibitem[Hochbruck and Ostermann ()]{b7}\label{b7} ‘Explicit Exponential Runge-Kutta Methods for Semi-linear Parabolic Problems’. M Hochbruck , A Ostermann . \textit{SIAM J. Numer. Anal} 2005. 43 p. . \bibitem[Cox and Matthews ()]{b8}\label{b8} ‘Exponential Time Differencing for Stiff Systems’. S M Cox , P C Matthews . \textit{J. Comput. Phys} 2002. 176 p. . \bibitem[Klein ()]{b4}\label{b4} ‘Fourth Order Time-Stepping for Low Dispersion Korteweg-de Vries and Nonlinear Schr??dinger Equations’. C Klein . \textit{Electronic Transactions on Numerica Analysis} 2008. 29 p. . \bibitem[Kassam and Trefethen ()]{b5}\label{b5} ‘Fourth-Order Time Stepping for Stiff PDEs’. A K Kassam , L N Trefethen . \textit{SIAM J. Sci. Comput} 2005. 26 p. . \bibitem[Friedli ()]{b2}\label{b2} ‘Generalized Runge-Kutta Methods for the Solution of Stiff Dierential Equations’. Friedli . \textit{Numerical Treatment of Dierential Equations}, Lecture Notes in Mathematics R Burlirsch, R Grigorie, J Schr??der (ed.) (Berlin) 1978. Springer. 631 p. . \bibitem[Kassam ()]{b9}\label{b9} \textit{High Order Time stepping for Stiff Semi-Linear Partial Differential Equations}, A K Kassam . 2004. Oxford University (PhD thesis) \bibitem[Certaine ()]{b1}\label{b1} ‘The Solution of Ordinary Dierential Equations with Large Time Constants’. J Certaine . \textit{Mathematical Methods for Digital Computers}, A Ralston, H S Wilf (ed.) (New York) 1960. Wiley. p. . \end{bibitemlist} \end{document}