代码如下：
solution=NDSolve[{D[u[x,t],t,t]+D[u[x,t],t]+D[u[x,t],x,x,x,x]-D[u[x,t],x,x]+0.01/(1+u[x,t])^2==0,u[x,0]==0.5,Derivative[0, 1][u][x, 0] == 0,u[0,t]==0,u[1,t]==0},u,{x,0,1},{t,0,1}];
Plot3D[Evaluate[u[x,t]/.First[solution]],{x,0,1},{t,0,1}];