[SciPy-user] first order pde

Marcos Capistran marcos.capistran@gmail....
Fri Sep 28 08:23:17 CDT 2007


Hi Issa,

You must make sure that your numerical method satisfies the
Courant-Friedrichs-Lewy condition. On the other hand, (implicit) Euler
is a low order method, it may not be a good choice if you need a reliable
solution after long time integration.

good look,

On 9/28/07, Issa Karambal <issa@aims.ac.za> wrote:
> Hi all,
>
> I am having  problem to solve numerically a hyperbolic equation like
> u_x=-pi*u_t, u(x,0)=exp(sin(x)), u(0,t)=u(2pi,t)=1. I used the method of
> lines; so I discretized the spatial domain using the centered finite
> difference method and thereafter I used the implicit Euler method, but I
> am still having problem after long time integration.
> If anyone has an idea on how I can solve numerically my problem?
>
> THANKS,
>
> issa
>
> _______________________________________________
> SciPy-user mailing list
> SciPy-user@scipy.org
> http://projects.scipy.org/mailman/listinfo/scipy-user
>

-- 
Marcos Aurelio Capistrán Ocampo
CIMAT
A. P. 402
Jalisco S/N, Valenciana
Guanajuato, GTO 36240
Tel: (473) 73 2 71 55 Ext. 49640


More information about the SciPy-user mailing list