4D plot around a cylinder

Hi, I have a formula for stress distribution around a borehole, sigtet(r,T)=(((sig1+sig2)*(1+R^2/a^2))/2)-(((sig1-sig2)*(1+(3*R^4)/a^4)*cos(T*2))/2)-sig12*(1+3*(R/a)^4)*sin(2*T)-delp*(R/a)^2

for a special depth I can enter sig1, sig2 and delp values and get the stress distribution around the wellbore by this code: R=10.5, sig1=100, sig2=35, sig12=0, delp=10,

    for ii=1:500 
    for jj=1:201
  a=R+(ii+1)/100; 
  T=((jj-1)/10)*pi;
  sigtet(ii,jj)=(((sig1+sig2)*(1+R^2/a^2))/2)-(((sig1-sig2)*(1+(3*R^4)/a^4)*cos(T*2))/2)-sig12*(1+3*(R/a)^4)*sin(2*T)-delp*(R/a)^2;
  x(ii,jj)=(a)*cos(T); 
  y(ii,jj)=(a)*sin(T);
  end 
  end 
  pcolor(x,y,sigtet) 
  shading interp 
but sig1 and sig2 and delp values are functions of depth (z), as u see for 2D plot we have stress values in different radius and angles, and the color shows stress intensity for a known depth, now I want to graph it in 3D plot which shows the stress distribution not only for different radius and angles but for different depths. so the new function will be:
sigtet(r,T,z)=(((sig1(z)+sig2(z))*(1+R^2/a^2))/2)-(((sig1(z)-sig2(z))*(1+(3*R^4)/a^4)*cos(T*2))/2)-sig12*(1+3*(R/a)^4)*sin(2*T)-delp(z)*(R/a)^2

Answer 1

[David Sanchez]

I am sorry to disappoint you, but the 4D plot is not yet available. You should choose 3 of your variables at a time and plot them together. Visit (type it)

doc plot3
doc mesh
doc surf

I am sure you will find a plotting function that suits your needs.

[Farid sjdclkjas]

thank you for your answer. here there are 3 variables, (a,T,z) and sigtet is a function of them, but I want to plot "sigtet" on a cylinder.is it possible?