>    restart;

>    u1 :=evalf(sum(14*(1-(-1)^n )/(n*Pi*sinh(n*Pi*2/3))*sinh(n*Pi/3)*sin(n*Pi/2) ,n=1..1)+sum(-2*(-1)^n/(n*Pi*sinh(n*Pi*3/2))*sinh(n*Pi/2)*sin(n*Pi*3/4),n=1..1));

u1 := 2.803326196

>    u2 :=evalf(sum(14*(1-(-1)^n )/(n*Pi*sinh(n*Pi*2/3))*sinh(n*Pi/3)*sin(n*Pi/2) ,n=1..2)+sum(-2*(-1)^n/(n*Pi*sinh(n*Pi*3/2))*sinh(n*Pi/2)*sin(n*Pi*3/4),n=1..2));

u2 := 2.803919511

>    u5 :=evalf(sum(14*(1-(-1)^n )/(n*Pi*sinh(n*Pi*2/3))*sinh(n*Pi/3)*sin(n*Pi/2) ,n=1..5)+sum(-2*(-1)^n/(n*Pi*sinh(n*Pi*3/2))*sinh(n*Pi/2)*sin(n*Pi*3/4),n=1..5));

u5 := 2.685272622

>    u10 :=evalf(sum(14*(1-(-1)^n )/(n*Pi*sinh(n*Pi*2/3))*sinh(n*Pi/3)*sin(n*Pi/2) ,n=1..10)+sum(-2*(-1)^n/(n*Pi*sinh(n*Pi*3/2))*sinh(n*Pi/2)*sin(n*Pi*3/4),n=1..10));

u10 := 2.684518156

>    u50:=evalf(sum(14*(1-(-1)^n )/(n*Pi*sinh(n*Pi*2/3))*sinh(n*Pi/3)*sin(n*Pi/2) ,n=1..50)+sum(-2*(-1)^n/(n*Pi*sinh(n*Pi*3/2))*sinh(n*Pi/2)*sin(n*Pi*3/4),n=1..50));

u50 := 2.684510862

>    u := (x,y) -> sum(14*(1-(-1)^n )/(n*Pi*sinh(n*Pi*2/3))*sinh(n*Pi*x/3)*sin(n*Pi*y/3) ,n=1..50)+sum(-2*(-1)^n/(n*Pi*sinh(n*Pi*3/2))*sinh(n*Pi*x/2)*sin(n*Pi*y/2),n=1..50) + sinh(4*Pi*(2-x)/3)*sin(4*Pi*y/3)/sinh(8*Pi/3) ;

u := proc (x, y) options operator, arrow; sum(14*(1-(-1)^n)/n/Pi/sinh(2/3*n*Pi)*sinh(1/3*n*Pi*x)*sin(1/3*n*Pi*y),n = 1 .. 50)+sum(-2*(-1)^n/n/Pi/sinh(3/2*n*Pi)*sinh(1/2*n*Pi*x)*sin(1/2*n*Pi*y),n = 1 .....

>    plot3d( u(x,y), x=0..2, y=0..3, labels =[`x`, `y`, `u`], axes=FRAMED, style=patch);

[Maple Plot]

>