how to draw a temperature distribution on (x,y,z)?

I have 3 coordinates x,y,z. For each x,y,z I have a corresponding temperature t. How can I draw it in matlab to get a "smooth" graph as in the below matlab example - graph (b)?

http://www.google.co.uk/imgres?um=1&newwindow=1&sa=N&rlz=1C1GPCK_en-GBGB495&hl=en&biw=1280&bih=667&tbm=isch&tbnid=YFrJj663YEdmJM:&imgrefurl=http://www.emeraldinsight.com/journals.htm%3Fissn%3D0961-5539%26volume%3D18%26issue%3D7/8%26articleid%3D1752627%26show%3Dhtml&docid=upUAHvGEujXwOM&itg=1&imgurl=http://www.emeraldinsight.com/content_images/fig/1340180707031.png&w=1405&h=608&ei=oLW0UeTQBvPK0AX9oICwDA&zoom=1&ved=1t:3588,r:6,s:0,i:99&iact=rc&dur=4201&page=1&tbnh=122&tbnw=282&start=0&ndsp=15&tx=234&ty=35

Answer 1

[Iman Ansari]

t defines color:

[x y] = meshgrid(-10:0.1:10);
z =zeros(size(x));
z(x+2*y>-9 & 2*x+y<9 & x-y>-10 & x-y<15) = 1;

t = - x + y;
s=surf(x,y,z,t);
set(s,'EdgeColor','none')

[D]

thank you Iman what is z in your answer? in my problem z is the 3rd coordinate. Could you explain please? thanks

[Iman Ansari]

Here z is 3rd coordinate too, x,y and z define the shape and t defines its color:

[x y] = meshgrid(-10:0.1:10);
z = zeros(size(x));
z(x+2*y>-9 & 2*x+y<9 & x-y>-10 & x-y<15) = 1;

t = - x + y;
t(z==1) = 0;
s=surf(x,y,z,t);
set(s,'EdgeColor','none')

[D]

Thank you Iman.

I used function scatter3, like this: scatter3(x,y,z,[2000],t,'.'); grid off box off view(32,18); axis equal tight off vis3d; % azimuth 26 camproj perspective camlight; lighting gouraud; alpha(0.75); rotate3d on

And I got these:

https://lh4.googleusercontent.com/fd7y3UlG6ADff97F7tqSqRB9SASpNgnQ4xo2OhZj5Uc=w209-h207-p-no

https://lh6.googleusercontent.com/cby9SBv6Gx_GvZ9ycuaQBmgtXmCDlaRZzOf6UZU6Hj8=w112-h207-p-no

The problem is that the temperature points are big (big circles) which i want to smooth out between other circles. If I make the cricles smaller, the distances between these circles are even bigger, and the picture looses its appearance:

https://lh6.googleusercontent.com/ULZA8UR6CcJ30Ux9AtNwYx5c_nYonzMM4s_Y7pl34ak=w112-h207-p-no

How to overcome this? Can I maybe use interpolation, so that empty spaces between the circles could be filled out with average values?