figure hold on for ij = 1:10 for jt = 1:100 u(4,ij,jt)= 2.*ij.^2+5.* jt.^3; plot3(ij,jt,u(4,ij,jt)) end end xlabel('Z');ylabel('T'); zlabel('u4');