V1
V1
In[]:=
t=0;A=2;c1=1;c2=1;c=1/25;omega=2;pts=Array[List/*N,101{1,1},{{-5,5},{-4,6}}];{xi,eta}={pts[[All,All,1]],pts[[All,All,2]]};P=2Sinh[xi+omegat]/(3Cosh[xi+omegat])+(5Sinh[xi+omegat]/(6Cosh[xi+omegat])^3)+0.9;Q=-Sinh[eta]/Cosh[eta];x=xi-2.5Tanh[xi+omegat];y=eta;Px=-Sech[xi+omegat]^2;Qy=-Sech[eta]^2;u=(A-c1c2)PxQy/(1+c1P+c2Q+cPQ)^2;pols=BlockMap[{{#[[1,1]],#[[1,2]],#[[2,1]]},{#[[2,1]],#[[1,2]],#[[2,2]]}}&,Transpose[{x,y,u},{3,1,2}],{2,2},{1,1}];pols=Catenate[Catenate[pols]];Graphics3D[Polygon[pols],BoxRatios{1,1,1}]
Out[]=
V2
V2
In[]:=
ClearAll[funcu]funcu[xi_,eta_]:=Module[{t=0,A=2,c1=1,c2=1,c=1/25,omega=2,P,Q,x,y,Px,Qy},P=2Sinh[xi+omegat]/(3Cosh[xi+omegat])+(5Sinh[xi+omegat]/(6Cosh[xi+omegat])^3)+0.9;Q=-Sinh[eta]/Cosh[eta];x=xi-2.5Tanh[xi+omegat];y=eta;Px=-Sech[xi+omegat]^2;Qy=-Sech[eta]^2;{x,y,(A-c1c2)PxQy/(1+c1P+c2Q+cPQ)^2}]ParametricPlot3D[funcu[xi,eta],{xi,-5,5},{eta,-4,6},BoxRatios{1,1,1},PlotRangeAll,PlotPoints80,ColorFunction->Function[{x,y,z,ϕ,θ},ColorData["BlueGreenYellow"][z]]]
Out[]=
