> | 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)); |
> | 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)); |
> | 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)); |
> | 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)); |
> | 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)); |
> | 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) ; |
> | plot3d( u(x,y), x=0..2, y=0..3, labels =[`x`, `y`, `u`], axes=FRAMED, style=patch); |
> |