}{2C_1-2C_2\cdot t}}\approx 1.0435$$ according to. This prediction is also shown on the figure as line (0).

From this figure we see that at high resolutions all schemes produce the same growth of amplitude of perturbations, which is in agreement with the analytical prediction. As the resolution decreases the results for all schemes for all schemes begin to deviate from the analytical value. Among the method of lines schemes, the RK4 is the most accurate and is able to reproduce the correct behavior of the perturbation up to the resolution of 16 grid points per wavelength. CFLN also requires only 16 grid points per wavelength to reproduce the correct behavior. The other two schemes requires from 32 to 64 grid points per wavelength. The major difference between the CFLN and the method of lines schemes is that at low resolutions, CFLN is amplifying the perturbations more than it should according to the analytical predictions, whereas the method of lines schemes are damping. The RK4 is the least damping and the most accurate of these schemes. The behavior of the schemes at low resolutions is fully consistent with the results of the Von Neumann stability analysis, presented in section \[sec:2.01\]. According to that analysis the CFLN scheme must be amplifying at all wavelengths whereas the method of lines schemes must be damping.

**Discussion and conclusion**
=============================

In this paper we studied the properties of four numerical schemes, CFLN, ICN, RK3 and RK4, applied to a non-linear scalar wave equation. This equation has a number of non-trivial analytic solutions whose properties, including stability, were studied and summarized in section \[sec:theory\]. We carried out the Von Neumann stability analysis of the schemes and studied their phase and amplitude errors. Finally we carried out numerical experiments and compared the results of those experiments with the perturbation analysis of the equation and the Von Neumann stability analysis of the schemes.

We find that all four schemes presented in this paper are stable and converge with the second order accuracy. The stability range for the schemes were determined and are presented in table \[tab:stableinterval\]. Those ranges are valid for both linear and non-linear solutions.

For non-linear schemes we checked that the amplification factor $|\lambda|$ is less than $1+O(\Delta t)$ for sufficiently small $\Delta t$.

With respect to the dissipation errors, we find that the CFLN has the least amount of dissipation
} { 2C_1 - 2C_2\cdot t}}\approx 1.0435$$ according to. This prediction is also shown on the figure as course (0). 

 From this name we see that at high resolutions all dodge produce the same emergence of amplitude of perturbation, which is in agreement with the analytical prediction. As the resolving power decreases the results for all schema for all schemes begin to deviate from the analytic value. Among the method acting of lines schemes, the RK4 is the about accurate and is able to reproduce the right demeanor of the perturbation up to the resolution of 16 grid points per wavelength. CFLN besides requires only 16 power system points per wavelength to regurgitate the correct behavior. The other two schemes requires from 32 to 64 grid point per wavelength. The major difference between the CFLN and the method of lines schemes is that at low resolution, CFLN is amplifying the perturbations more than it should according to the analytical predictions, whereas the method of lines schemes are damping. The RK4 is the least damping and the about accurate of these scheme. The behavior of the system at depleted resolutions is fully consistent with the results of the Von Neumann stability analysis, presented in incision \[sec:2.01\ ]. According to that analysis the CFLN scheme must be amplifying at all wavelengths whereas the method of line schemes must be damping. 

 * * Discussion and conclusion * * 
 = = = = = = = = = = = = = = = = = = = = = = = = = = = = = 

 In this paper we study the property of four numerical schemes, CFLN, ICN, RK3 and RK4, applied to a non - linear scalar wave equation. This equality has a number of non - trivial analytic solutions whose properties, including stability, were studied and summarized in department \[sec: theory\ ]. We carried out the Von Neumann stability psychoanalysis of the schemes and studied their phase and amplitude errors. Finally we behave out numerical experiment and compared the results of those experiment with the disruption analysis of the equation and the Von Neumann stability analysis of the schemes. 

 We find that all four schemes presented in this paper are static and converge with the second order accuracy. The stability scope for the schemes were determined and are presented in table \[tab: stableinterval\ ]. Those ranges are valid for both linear and non - linear solution. 

 For non - linear schemes we checked that the amplification factor $ |\lambda|$ is less than $ 1+O(\Delta t)$ for sufficiently small $ \Delta t$. 

 With respect to the dissipation error, we witness that the CFLN has the least amount of profligacy
}{2C_1-2C_2\cfot t}}\approx 1.0435$$ according tu. This prediction is anso shkwn on tfe figure as line (0).

From this higuee we see that at high resouutions apl schemws pcoduce the same jdowth on ampmltude if perturbatioks, which is in agreement figh the analytical prediction. As the rqsolutipn decreases the reslles fkg cll schemes for all schemes begjn to dtviate from the anslytical value. Among the mfthof of lines schemes, the RK4 is jge iist accurate and is able to reprodoce the correct behavior of the oertuxbation up jo tjg resolution of 16 drid points ivr wavenength. VFLN also requlres mnlt 16 grid points per watelength to reproducg the corracc behavior. The other rwi schgmes sequkees frkm 32 fo 64 grld 'oints per savelength. Rhe major differenct bqneeen the CFLH and ehq method of lines schemes is that at lof rssolutions, CFLN is amplufying the perturbatilns more ehan it should according to the analytical predicdions, dhexcwr tje method of lines schemes are damping. The RK4 js tme least damping and the moxt avsurate of there schzjea. The behavior of hhe schgmes ar low resjlutoons is fully consistent wirh the resulns od the Von Neumann dtability aualysix, prexented in section \[sec:2.01\]. Aecordihg to that wnalysis ffe CFLN scheme mjst ba amplifying at all waveleggths wheceas che methud og linef schemes lust na damping.

**Discussioj and cmnclusion**
=============================

Ij this paper we studied the pro'xrties of fout ngmegical schzmes, CNLN, ICN, RK3 and WK4, applied to c non-linzar scxlar wave vquation. Vhis equatiog has a numbes of non-triviel analytyc silutuons whure properties, oncluding stability, were studied and xumozrized in sectiin \[wec:theory\]. We catridd jun tie Vog Neumann statiligy xmalysks of the sgheoes snd studied their phdse znd amplitude errots. Finally we carryed out numerocal experiments ajd cokpaced thr rgsults of those experiments witg the perhurnation analysys on thg equation cnd the Von Neumann stability analysis oh the schemes.

We find thqt all four schemes kresented in vhis pwper are vtable and converge qith the second ovder accuracy. The stabjlity sange for the schemes were determined and are presented in table \[tab:stableinteeval\]. Tiofe ranges ade vslid xox bpth ligear end non-linear solltions.

For non-linear schemes we ciecked thad che amplification factor $|\lambca|$ is less than $1+O(\Delta t)$ for sufficientmy small $\Delta t$.

With respect to the disxipation errors, we find bhat the XFLN gas thv lxast amount of dissilation
}{2C_1-2C_2\cdot t}}\approx 1.0435$$ according to. This prediction shown the figure line (0). From at resolutions all schemes the same growth amplitude of perturbations, which is in with the analytical prediction. As the resolution decreases the results for all schemes all schemes begin to deviate from the analytical value. Among the method of schemes, RK4 the accurate and is able to reproduce the correct behavior of the perturbation up to the resolution 16 grid points per wavelength. CFLN also requires 16 grid points per to reproduce the correct behavior. other schemes requires 32 64 points per wavelength. major difference between the CFLN and the method of lines schemes is that at low resolutions, CFLN amplifying the than it according the predictions, whereas the lines schemes are damping. The RK4 damping and the most accurate of these schemes. behavior of schemes at low resolutions is fully with the results of the Von Neumann stability presented in section \[sec:2.01\]. According to that analysis the CFLN scheme must be amplifying at whereas the method of schemes must be **Discussion conclusion** In paper we the properties of four numerical schemes, CFLN, ICN, RK3 and RK4, to a non-linear scalar wave equation. This equation has a non-trivial solutions whose properties, stability, were studied and in \[sec:theory\]. We carried out Neumann analysis and their and amplitude errors. Finally carried out numerical experiments and the results of those of the equation and the Von Neumann stability of the schemes. We find that all schemes presented in this paper are stable and converge with the second accuracy. The for the schemes were determined and are presented table \[tab:stableinterval\]. Those ranges valid for both linear and non-linear solutions. For non-linear we that the factor $|\lambda|$ is than $1+O(\Delta t)$ sufficiently small $\Delta respect to dissipation we CFLN has the least amount of
}{2C_1-2C_2\cdot t}}\approx 1.0435$$ according to. THis predictIon is AlsO shOwN on tHe fiGure as line (0).

From THis fIgure we see that at high reSolutIoNS all SChEmes pRoduce tHE sAME grOwTh Of aMpLItUde of PerTurbatiOns, which is In aGrEement with thE AnAlytical prEdiCtion. As the reSolUtion dEcReaSEs the ResUlts fOr all sCHemes fOr all scheMeS Begin tO Deviate FROm The aNalytical value. AmoNG tHE method of lines SchemeS, tHE Rk4 IS thE moSt accurate AnD is abLE to reprODuCE THe cORrect behavior Of the perturBAtiOn up to ThE reSOlutioN of 16 grId POinTs per waveleNgth. cFLN also rEquireS Only 16 griD Points pEr waveLenGth To rePRoDuCe tHe COrrECt BehAVioR. The otheR tWo SchemEs reQUIREs frOm 32 tO 64 griD poinTs per wavelengTh. THe maJOr dIfferEnce bEtweEn The CFlN and tHe metHoD of lines schemes Is thAt at low reSolUtIonS, CfLN is AMplifyIng The PerturbAtions mORe tHaN IT ShOuld according to the AnALYtIcal predIctionS, WhErEAs the metHoD of LineS SChemeS are DAmPing. The Rk4 is the LEaSt Damping AnD the moSt AccUraTe of tHEse sChemes. the behavIor of THe schemes at low REsolutions is fULlY COnSIsteNt wIth the resulTs of THe VoN NeuMAnN stABilitY analYsIS, pREsented in section \[sec:2.01\]. acCordinG to thAt analysis the cFLN scheme MUST be ampliFyinG At ALl wavelengths wHereaS the method OF lines scHemes Must be daMping.

**DiscUSSion and cOncLusIon**
=============================

in tHIS pAper we studied THE proPeRties of FouR numeriCal SchEmeS, CFlN, iCN, RK3 and Rk4, applied To A nOn-LiNeaR scalAR wave equAtIon. thIs eQuatiON has a nUmber Of noN-tRiVIal AnalytiC SoLUTionS wHoSe prOpeRtIes, inCludINg sTabilitY, were studIed ANd suMmArIzed in sEction \[sec:theoRy\]. we carried oUt The von NeuMANn stabilIty analysis of the schemes ANd studiEd tHeir pHase And amplitUde Errors. finALly we cArried Out nuMeRicAL ExperIMEnTs aNd Compared thE REsuLts of ThOse eXperimeNts with the perturbaTIon Analysis of the EquAtioN ANd The vOn nEumAnN StaBILity analysis of tHe schemes.

WE fINd That all fouR SchEmEs preseNted in tHis paPEr are stAble and coNverge witH tHe seCONd oRder accuraCy. The staBility ranGE for tHE sChemeS weRe deteRmIneD and aRe presENteD in taBle \[tab:StAbleinTervaL\]. THose rangEs are valid for both linear And non-LineaR soLutions.

FoR noN-LinEar schemeS we cHecked that The AmpLificAtiON factOr $|\laMBdA|$ is LEss thAn $1+O(\DELta t)$ for suFFiCieNTLy Small $\Delta t$.

wITH reSpect To tHE dissiPatiOn errors, we find thaT The CFLN has the lEast AMOunT of DIssiPaTion
}{2C_1-2C_2\cdot t}}\appro x 1.0435$$  acco rdi ngto . Th is p rediction is a l so s hown on the figure asline(0 ) .

F r om  this  figure we s eeth at  at  h i gh  reso lut ions al l schemespro du ce the sameg ro wth of amp lit ude of pertu rba tions,  w hic h  is i n a greem ent wi t h theanalytica lp redict i on. Ast h ereso lution decreasest he results for al l sche me s  f o r  al l s chemes beg in  to d e viate f r om t h e a n alytical valu e. Among th e  me thod o flin e s sche mes,th e  RK 4 is the mo st a ccurate a nd isa ble tor eproduc e thecor rec t be h av io r o ft hep er tur b ati on up to  t he  reso luti o n o f 16  gr id p oints  per waveleng th.  CFL N  al so re quire s on ly  16 g rid po intspe r wavelength to  rep roduce th e c or rec tbehav i or. Th e o the r two s chemesr equ ir e s fr om 32 to 64 grid p oi n t sper wave length .  T he major di ff ere nceb e tween  the CF LN and t he met h od  o f lines  s chemes  i s t hat  at l o w re soluti ons, CFL N isa mplifying thep erturbationsm or e th a n it  sh ould accord ingt o th e an a ly tic a l pre dicti on s ,w hereas the method o flinesschem es are dampin g. The RK4 i s  the lea st d a mp i ng and the mos t acc urate of t h ese sche mes.The beha vior of t h e  schemes  at  lo w r eso l u ti ons is fullyc o nsis te nt with  th e resul tsofthe  Vo nNeumann s tability  a na ly si s,prese n ted in s ec tio n\[s ec:2. 0 1\]. A ccord ingto  t h atanalysi s  t h e  CFL Nsc heme  mu st  be a mpli f yin g at al l wavelen gth s  whe re as  the me thod of lines  s chemes mus tbedampin g . 

**Disc ussion and conclusion** ======= === ===== ==== ========= =
In thi s p a per we  studi ed th epro p e rties o ffou rnumericals c hem es, C FL N, I CN, RK3  and RK4, appliedt o a  non-linear s cal ar w a v eequ a ti o n.Th i s e q u ation has a num ber of non -t r iv ial analyt i c s ol utionswhose p roper t ies, in cluding s tability,  w eres t udi ed and sum marizedin sectio n  \[se c :t heory \].  We ca rr ied  outthe Vo n  Ne umann  stabi li ty ana lysis  o f the sc hemes and studied their  phase  andamp litude er ror s . F inally we  car ried out n ume ric al ex per i ments  and co mpa r ed th e re s ults of t h os e e x p er iments with t h e p ertur bat i on ana lysi s of the equation and the Von Ne uman n sta bil i ty a na lysis of the s che me s . 

We fin dthat all fo ur schem es prese nted i n this  papera r es tableandcon verge wit h t he secondor de r  accur acy.  T he sta bility rang e for the schemeswered e termi n edand a re  presen t ed i n table \[ tab:stablei nterva l\].  Thos e range sare va lid  f or both li n ear and n on-li near so lu tion s.
For n on-l i n ear s chem es  we  checkedt h at th ea mpl ific ation  f acto r $|\lamb d a|$ is l ess than $1 +O (\D e l ta t)$ fo r sufficient lysmall $ \Delta t$. 
Wit h  r e spect  to th e diss ipation err or s, we f ind t hat the C FLN has t h e l ea st a mount of  d issi pa tio n
}{2C_1-2C_2\cdot t}}\approx_1.0435$$ according_to. This prediction is_also shown_on_the figure_as_line (0).

From this_figure we see_that at high resolutions_all schemes produce_the_same growth of amplitude of perturbations, which is in agreement with the analytical prediction._As_the resolution_decreases_the_results for all schemes for_all schemes begin to deviate_from the_analytical value. Among the method of lines schemes,_the_RK4 is the_most accurate and is able to reproduce the correct_behavior of the perturbation up to_the resolution of_16_grid_points per wavelength. CFLN_also requires only 16 grid points_per wavelength to reproduce the correct_behavior. The other two schemes requires from_32 to 64 grid points per_wavelength. The major difference between_the CFLN_and the method of lines_schemes is that_at low_resolutions, CFLN is_amplifying the perturbations more than it_should according to_the analytical predictions, whereas the method_of_lines schemes are_damping._The_RK4 is_the least damping_and_the most_accurate_of these schemes. The behavior of_the_schemes at low resolutions is fully consistent_with the results of_the_Von Neumann stability analysis,_presented in section \[sec:2.01\]. According_to that analysis the CFLN scheme_must be_amplifying at_all wavelengths whereas the method of lines schemes must be damping.

**Discussion_and conclusion**
=============================

In this paper we studied_the properties of four_numerical schemes,_CFLN,_ICN, RK3 and_RK4,_applied to_a non-linear scalar wave equation. This equation_has a_number of non-trivial analytic solutions whose_properties, including stability, were_studied_and summarized in section \[sec:theory\]. We_carried out the Von Neumann stability_analysis of the schemes and_studied_their_phase and amplitude errors. Finally_we carried out numerical experiments and_compared the results_of those experiments with the perturbation analysis_of_the equation and the Von Neumann_stability_analysis of the schemes.

We find that_all_four_schemes presented in this paper_are stable and converge with the_second order accuracy. The stability range for the schemes_were determined and_are presented in table \[tab:stableinterval\]._Those_ranges_are valid for both linear and non-linear solutions.

For non-linear schemes_we checked_that the amplification_factor $|\lambda|$ is less than $1+O(\Delta t)$ for sufficiently small_$\Delta t$.

With respect to the dissipation errors,_we find that the CFLN has the least amount of dissipation