• API
• FAQ
• Tools
• Archive
daily pastebin goal
43%
SHARE
TWEET

# Untitled

a guest Feb 21st, 2018 50 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
1. # Pkg.add("Decimals") # Not needed because this is on juliabox
2.
3. using Decimals
4.
5. # Defining Variables
6. a = Decimal(0.15)
7. b = Decimal(0.2)
8. c = Decimal(1)
9. d = Decimal(-2)
10. n = Decimal(-0.5)
11. t = 0
12. values = []
13. column = Dict()
14.
15. # Creating the first column
16. while n < Decimal(0.5)
17.     column[n] = a + b*n^1 + c*n^2 + d*n^3
18.     n = n + Decimal(0.02)
19. end
20.
21. # Add the first column to the array of values
22. push!(next, column)
23.
24. # Starting an infinite loop
25. while true
26.     n = Decimal(-0.5)
27.     column = Dict()
28.
29.     # Creating the n'th column
30.     while n < Decimal(0.5)
31.         index = string(n)
32.         lastIndex = string(n-Decimal(0.02))
33.         nextIndex = string(n+Decimal(0.02))
34.         column[index] = values[end][index] + (values[end][lastIndex]*phi - 2*values[end][n]*phi + values[end][nextIndex]*phi)
35.
36.         if values[end][n] == values[end][lastIndex]
37.             break
38.         end
39.
40.         n = nextIndex
41.     end
42.
43.     if values[end][string(n)] == values[end][string(n-Decimal(0.02))]
44.       break
45.     end
46.
47.     push!(values, column)
48.     t = t + 4
49. end
50.
51. println(done!)
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy.

Top