for i=1:101
   for j=1:101
       V(i)=3+0.05*(i-1);
       Mmax(j)=0.09-0.0005*(j-1); 
       t=p(101*(i-1)+j);
       [V,Mmax]=meshgrid(3:0.05:8,0.09:-0.0005:0.04);
       surf(V(i),Mmax(j),t);
    end
end 
Warning: Z must not be a scalar or vector, not rendering surface.