-
Notifications
You must be signed in to change notification settings - Fork 2
Expand file tree
/
Copy pathPROJECT_LANG_1.rdata
More file actions
32 lines (27 loc) · 1.52 KB
/
Copy pathPROJECT_LANG_1.rdata
File metadata and controls
32 lines (27 loc) · 1.52 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
#!/usr/bin/env Rscript
# Project language file 1
# For: AI2001
# About
# I decided to make R the 1st and main project langugae file for this project (AI2001) as R is the perfect language to use for data mining and related Artificial Intelligence work, but will be backed up with Python, among other languages. It is getting its own project language file, starting here.
# Now time for a Mandelbrot (From Wikipedia/wiki/R/)
install.packages("caTools") # install external package
library(caTools) # external package providing write.gif function
jet.colors <- colorRampPalette(c("red", "blue", "#007FFF", "cyan", "#7FFF7F",
"yellow", "#FF7F00", "red", "#7F0000"))
dx <- 1500 # define width
dy <- 1400 # define height
C <- complex(real = rep(seq(-2.2, 1.0, length.out = dx), each = dy),
imag = rep(seq(-1.2, 1.2, length.out = dy), dx))
C <- matrix(C, dy, dx) # reshape as square matrix of complex numbers
Z <- 0 # initialize Z to zero
X <- array(0, c(dy, dx, 20)) # initialize output 3D array
for (k in 1:20) { # loop with 20 iterations
Z <- Z^2 + C # the central difference equation
X[, , k] <- exp(-abs(Z)) # capture results
}
write.gif(X, "Mandelbrot.gif", col = jet.colors, delay = 100)
# File info
# File type: R Source file (*.r *.rdata *.rds *.rda)
# File version: 1 (2023, Tuesday, February 21st at 9:30 pm PST)
# Line count (including blank lines and compiler line): 33
# End of script