>    restart;

>    u := (x,y) -> 3*cos(3*Pi*y)*sin(Pi*x) - 5*cos(48*Pi*y)*sin(16*Pi*x) + 4*sin(12*Pi*y)*sin(3*Pi*x);
plot3d( u(x,y), x=0..2*Pi, y=0..0.5, labels =[`x`, `time`, `u(x,t)`], axes=FRAMED, style=patch);

u := proc (x, y) options operator, arrow; 3*cos(3*Pi*y)*sin(Pi*x)-5*cos(48*Pi*y)*sin(16*Pi*x)+4*sin(12*Pi*y)*sin(3*Pi*x) end proc

[Maple Plot]

>    u1 := (x,y) -> 3*cos(3*Pi*y)*sin(Pi*x);
plot3d( u1(x,y), x=0..2*Pi, y=0..0.5, labels =[`x`, `time`, `u(x,t)`], axes=FRAMED, style=patch);

u1 := proc (x, y) options operator, arrow; 3*cos(3*Pi*y)*sin(Pi*x) end proc

[Maple Plot]

>    u2 := (x,y) -> - 5*cos(48*Pi*y)*sin(16*Pi*x) ;
plot3d( u2(x,y), x=0..2*Pi, y=0..0.5, labels =[`x`, `time`, `u(x,t)`], axes=FRAMED, style=patch);

u2 := proc (x, y) options operator, arrow; -5*cos(48*Pi*y)*sin(16*Pi*x) end proc

[Maple Plot]

>   

>    u3 := (x,y) -> 4*sin(12*Pi*y)*sin(3*Pi*x);
plot3d( u3(x,y), x=0..2*Pi, y=0..0.5, labels =[`x`, `time`, `u(x,t)`], axes=FRAMED, style=patch);

u3 := proc (x, y) options operator, arrow; 4*sin(12*Pi*y)*sin(3*Pi*x) end proc

[Maple Plot]

>