2nd resolution Wave_t_ex

[Lin Ji]

Hi, Christophe,
    I did not see the big warning about the second scheme in your last email.
Now, I read about it. But, I still don't quite understand how to assemble the
source. I do express TimeFct[] and dfdt[] in the same order as in your
example:

  yo[]=1/(ct*cy*Sqrt[2*3.1416])*Exp[-(($Time-t0)^2/(ct^2)+(Y[]-y0)^2/(cy^2)
       )];
  yo2[]=Cos[gr*($Time-t0)]*yo[];
  TimeFct[] = yo[];
  // TimeFct[] = ( (X[]<.03)  )? yo[] : 0 ;

  dfdt[] = TimeFct[] ;

  SaveFct[] = !($TimeStep % 2) ;

instead of 'TimeFct[] = dfdt[];'. It seems I have to set 'TimeFct[] =
dfdt[]/dfdt[t==t0]' in order to make 'dfdt[t==t0]*TimeFct[]' the correct
source as you explained in the email. Does that mean dfdt[t==t0] can not be
zero? Please help.

Lin

-- 
Dept. of Math., RPI
518-276,2184(work)
518-271,4486(home)