Solución de u_t(x,t) =u_{xx}(x,t),   si x en (0,2), t en (0,3],

    u(0,t)=u(2,t)= 0,   si t en (0,3],

    u(x,0)=1,   si x en (0,2),

    usando series de Fourier.

 
>    u := sum('2*(1-(-1)^n)/((n*Pi))*sin(n*Pi*x/2)*exp(-t*(n*Pi/2)^2)', 'n'=1..500): plot3d(u, x = 0..2, t = 0..3);

[Maple Plot]

 

    Solución de u_t(x,t) =u_{xx}(x,t),   si  x en (0,2), t en (0,3],

    u_x(0,t)=u_x(2,t)= 0,  si t en (0,3],

    u(x,0)=x,    si x en (0,2),

    usando series de Fourier.

 

 
>    u := 1+sum('4*((-1)^n-1)/((n*Pi)^2)*cos(n*Pi*x/2)*exp(-t*(n*Pi/2)^2)', 'n'=1..100): plot3d(u, x = 0..2, t = 0..3);

[Maple Plot]