The command I used is plot(plot::Implicit3d(x*1000-3.14*2^2*(1+y/1000)^2*544*(ln(5/z))^0.23*5/z, x=15..25, y=-15..10, z=0.5..4.5), Scaling = Constrained).
I'm sure that the value of the function falls into the ranges of x, y and z
Your command works fine in the latest version, R2014b. What version of Matlab/MuPAD are you using (type version into your command window)? The Symbolic Math toolbox gained many new features starting with R2012a. In particular, R2012a gained relational operators and assumptions on variables:
Other functions in the toolbox were likely updated to use these in this and subsequent releases.
Also, if you are using an older version, be wary of what you come across in Google searches as The MathWorks posts the documentation for the latest release. Make sure that you either use the help or doc included with the install or or use the archived online documentation specifically for your release: http://mathworks.com/help/doc-archives.html
Thanks for your answer, Andrew, unfortunately I failed both in MATLAB 2010b and version 2012a. I wonder if this error is caused by an older version, maybe I should try it in the latest version as your suggestion.