Power function

[eoredson]

This attached file (by STxAxTIC) is called sxmath with tweaks for powers and fractional powers to any digit:

Erik.

[STxAxTIC]

Say, thanks for the plug! I'll have to drag that project back out into public display one of these days...

[_vince]

series experiment

Code: QB64: [Select]

1. deflng a-z
2. screen _newimage(800,600,12)
3.  
4. dim x as double, y as double
5. dim t as double
6.  
7. do
8.         t = t + 0.1
9.         cls
10.         for xx=1 to 800
11.                 x = 5*xx/800
12.                 y = x^t
13.                 yy = y
14.                 pset (xx, 600-yy), 7
15.  
16.                 y = ppow#(x, t, 2)
17.                 yy = y
18.                 pset (xx, 600-yy), 14
19.  
20.                 y = ppow#(x, t, 3)
21.                 yy = y
22.                 pset (xx, 600-yy), 12
23.  
24.                 y = ppow#(x, t, 4)
25.                 yy = y
26.                 pset (xx, 600-yy), 10
27.  
28.                 y = ppow#(x, t, 5)
29.                 yy = y
30.                 pset (xx, 600-yy), 11
31.         next
32.         color 7:locate 1,1:?"x^"+ltrim$(str$(t))
33.         color 14:locate 2,1:?"n=2"
34.         color 12:locate 3,1:?"n=3"
35.         color 10:locate 4,1:?"n=4"
36.         color 11:locate 5,1:?"n=5"
37.         if t>10 then t=0
38.         _display
39.         _limit 10
40. loop until _keyhit=27
41. system
42.  
43. function ppow#(c as double, x as double, n)
44.         dim cc as double, xx as double, p as double
45.         cc = exp(c*log(c))
46.         p = 1
47.         xx = 1
48.         for i=1 to n
49.                 xx = xx*log(c)*(x-c)/fact(i)
50.                 p = p + xx
51.         next
52.         ppow# = cc*p
53. end function
54.  
55. function fact(n)
56.         fact = 1
57.         for i=0 to n-1
58.                 fact = fact*(n-i)
59.         next
60. end function
61.