generated from seanpm2001/Git-Template_V8
-
-
Notifications
You must be signed in to change notification settings - Fork 2
/
PROJECT_LANG_2.r
29 lines (27 loc) · 1.38 KB
/
PROJECT_LANG_2.r
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
#!/usr/bin/env Rscript
# Project language file 2
# For: DIFFicult
# About
# I decided to make R the 2nd project language file for this project (DIFFicult) as I wanted to use the language to program the UI for this project, as part of software diversity reasons.
# 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)
# File version: 1 (2022, Monday, January 3rd at 4:36 pm)
# Line count (including blank lines and compiler line): 30
# End of script