home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #20 / NN_1992_20.iso / spool / sci / math / symbolic / 2384 < prev    next >
Encoding:
Text File  |  1992-09-11  |  2.6 KB  |  76 lines
  1. Newsgroups: sci.math.symbolic
  2. Path: sparky!uunet!cs.utexas.edu!wupost!darwin.sura.net!Sirius.dfn.de!chx400!bernina!neptune!monagan
  3. From: monagan@inf.ethz.ch (Michael)
  4. Subject: Re: plots and eval
  5. Message-ID: <1992Sep11.135137.22089@neptune.inf.ethz.ch>
  6. Summary: plot3d functionality in Maple
  7. Sender: news@neptune.inf.ethz.ch (Mr News)
  8. Nntp-Posting-Host: rutishauser-gw.inf.ethz.ch
  9. Organization: Dept. Informatik, Swiss Federal Institute of Technology (ETH), Zurich, CH
  10. References: <1992Sep10.164017.3507@magnus.acs.ohio-state.edu>
  11. Date: Fri, 11 Sep 1992 13:51:37 GMT
  12. Lines: 62
  13.  
  14. In article <1992Sep10.164017.3507@magnus.acs.ohio-state.edu>, hmarvel@magnus.acs.ohio-state.edu (Howard P Marvel) writes:
  15. > Consider the following series of statements in Maple:
  16. > 
  17. > c := proc (alpha,beta) if signum(alpha+beta)=1 then alpha+ beta else 0 fi en\
  18. > d;
  19. > c(a,b);
  20. > a:=1;
  21. > b:=1;
  22. > c(a,b);
  23. > plot3d(c(x,y),x=0..1,y=0..1);
  24. > 
  25. > The first sets up a proc that returns alpha + beta if alpha + beta is
  26. > a positive real. However, if alpha and beta are variables, signum
  27. > cannot evaluate the expression, and the proc returns 0. The first call
  28. > c(a,b) returns 0 and the second returns 2. Now consider the plot
  29. > statement. This should yield values from 0 to 2, right? Instead,
  30. > signum fails and all plotted values are 0. My question is, is there a
  31. > way to force Maple to evaluate c for the values in question and to
  32. > plot the values?  
  33. > 
  34. > Thanks in advance for any suggestions.
  35.  
  36. Dear Howard,
  37.  
  38. The plot3d function has no special evaluation rules.  The arguments
  39. get evaluated, then function is executed.  So what happens when you do
  40.  
  41.     plot3d( c(x,y), x=0..1, y=0..1 );
  42.  
  43. is that c(x,y) is executed on symbolic x and y, and this returns 0,
  44. so it's just as if you had written  plot3d( 0, x=0..1, y=0..1 );
  45. The plot3d command has two ways of being called.  You can plot either
  46. a formula in x and y, or a function of two variables.  Consider
  47.  
  48.     plot3d( binomial(n,k), n=0..3, k=0..3 );
  49.  
  50. versus
  51.  
  52.     plot3d( binomial, 0..3, 0..3 );
  53.  
  54. The second way allows you to plot a function of two variables defined
  55. by a Maple procedure or operator expression.  So the simple solution
  56. to your problem is
  57.  
  58.     plot3d( c, -1..1, -1..1 );
  59.  
  60. The alternative is to use quotes to prevent c(x,y) being evaluated
  61. before the plot routine gets the expression, i.e.
  62.  
  63.     plot3d( 'c(x,y)', x=-1..1, y=-1..1 );
  64.  
  65. Some more examples to clarify.
  66.  
  67.     plot3d( 'binomial(x,y)*c(x,y)', x=-3..3, y=-3..3 );
  68.     plot3d( binomial*c, -3..3, -3..3 );
  69.  
  70. These are wrong
  71.  
  72.     plot3d( 'binomial(x,y)*c(x,y)', -3..3, -3..3 );
  73.     plot3d( binomial*c, x=-3..3, y=-3..3 );
  74.  
  75. Mike Monagan
  76.