home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1993 #1 / NN_1993_1.iso / spool / comp / lang / pascal / 7957 < prev    next >
Encoding:
Internet Message Format  |  1993-01-07  |  2.2 KB
  1. Path: sparky!uunet!munnari.oz.au!ariel.ucs.unimelb.EDU.AU!ucsvc.ucs.unimelb.edu.au!lugb!news
  2. Newsgroups: comp.lang.pascal
  3. Subject: Re: powers in TURBO pascal
  4. Message-ID: <1993Jan8.011629.16911@lugb.latrobe.edu.au>
  5. From: CHEAC@LURE.LATROBE.EDU.AU (CZERWINSKI,Alex)
  6. Date: Fri, 8 Jan 1993 01:16:29 GMT
  7. Sender: news@lugb.latrobe.edu.au (USENET News System)
  8. References: <30DEC92.02825818.0033@music.mus.polymtl.ca> <dmurdoch.262.725736007@mast.queensu.ca> <1992Dec30.174246.2523@news.uni-stuttgart.de>
  9. Organization: La Trobe University
  10. In-Reply-To: lehnert@ifsws1.sozialforschung.uni-stuttgart.de's message of Wed, 30 Dec 1992 17:42:46 GMT
  11. X-News-Reader: VMS NEWS 1.24
  12. Lines: 48
  13.  
  14. In <1992Dec30.174246.2523@news.uni-stuttgart.de> lehnert@ifsws1.sozialforschung.uni-stuttgart.de writes:
  15.  
  16. > Hi, everybody,
  17. > I suppose it's trivial, but ...
  18. > 
  19. > What is the TURBO pascal operator for "x to the power of y"?
  20. > 
  21. > I tried ^ (reserved for pointers), **, and others,
  22. > and consulted the Borland Reference Books;
  23. > Maybe somebody can help me.
  24. > 
  25. > Thanks,
  26. > Doris
  27. > -- 
  28. > Doris Lehnert            lehnert@ifsws1.sozialforschung.uni-stuttgart.de
  29. > Universitaet Stuttgart, Institut fuer Sozialforschung
  30.  
  31.  Here is an answer by Prof. Timo Salmi :
  32.  
  33.  
  34.  Q: I need a power function but there is none in Turbo Pascal.
  35.  
  36.  A: Pascals do not have an inbuilt power function. You have to write
  37. one yourself. The common, but non-general method is defining
  38.    function POWERFN (number, exponent : real) : real;
  39.      begin
  40.        powerfn := Exp(exponent*Ln(number));
  41.      end;
  42. To make it general use:
  43.    (* Generalized power function by Prof. Timo Salmi *)
  44.    function GENPOWFN (number, exponent : real) : real;
  45.    begin
  46.      if (exponent = 0.0) then
  47.        genpowfn := 1.0
  48.      else if number = 0.0 then
  49.        genpowfn := 0.0
  50.      else if abs(exponent*Ln(abs(number))) > 87.498 then
  51.        begin writeln ('Overflow in GENPOWFN expression'); halt; end
  52.      else if number > 0.0 then
  53.        genpowfn := Exp(exponent*Ln(number))
  54.      else if (number < 0.0) and (Frac(exponent) = 0.0) then
  55.        if Odd(Round(exponent)) then
  56.          genpowfn := -GENPOWFN (-number, exponent)
  57.        else
  58.          genpowfn :=  GENPOWFN (-number, exponent)
  59.      else
  60.        begin writeln ('Invalid GENPOWFN expression'); halt; end;
  61.    end;  (* genpowfn *)
  62.