Julia set: Difference between revisions

{{task}}

 

;Task

*   [[Mandelbrot_set|Mandelbrot Set]]

<br><br>

 

=={{header|AWK}}==

{{trans|COBOL}}

# syntax: GAWK -f JULIA_SET.AWK [real imaginary]

BEGIN {

exit(0)

}

{{out}}

<pre>

 

=={{header|BASIC}}==

 

==={{header|Sinclair ZX81 BASIC}}===

70 FOR A=1 TO 50

80 LET R=X*X-Y*Y+C

100 IF R*R>1000 THEN GOTO 150

110 LET X=R

120 LET Y=I

130 NEXT A

150 NEXT H

 

==={{header|ZX Spectrum Basic}}===

{{trans|Sinclair ZX81 BASIC}}

Higher resolution is obtainable, if you have the time to wait for it.

20 LET cimag=0.156

30 FOR v=-16 TO 16

140 PLOT h+100,150-v

150 NEXT h

{{out}}

Screenshot [http://edmundgriffiths.com/spectrumjulia.jpg here].

 

=={{header|C++}}==

[[File:JuliaSetCpp.png|200px|thumb|right]]

#include <windows.h>

#include <string>

julia j; j.draw( c ); return 0;

}

 

 

 

 

}

}

}

}

 

}

 

=={{header|COBOL}}==

Plots—in ASCII or EBCDIC art—a Julia set for the function <i>f</i>(<i>z</i>) = <i>z</i>² + <i>c</i>, based on a value of <i>c</i> input by the user (real part then imaginary part, pressing the carriage return key after each). The sample output is for the inputs <tt>-0.8</tt> and <tt>0.156</tt>.

PROGRAM-ID. JULIA-SET-PROGRAM.

DATA DIVISION.

MOVE SPACE TO WS-PLOT-CHARACTER.

MOVE Z-REAL TO X.

{{out}}

<pre>

</pre>

 

=={{header|Elixir}}==

{{trans|AWK}}

def set(c_real, c_imag) do

IO.puts "#{c_real}, #{c_imag}"

c_real = if r=Enum.at(System.argv, 0), do: Float.parse(r) |> elem(0), else: -0.8

c_imag = if c=Enum.at(System.argv, 1), do: Float.parse(c) |> elem(0), else: 0.156

 

{{out}}

</pre>

 

=={{header|J}}==

[[File:example-j-julia.png|200px|thumb|right]]

pal2=: 255,~0,<.(254$1 0.8 0.6)*Hue 5r6*(i.%<:)254

g=: [: %: 0.3746j0.102863 0.132565j0.389103 _0.373935j_0.353777 1&p.

 

See also: [http://webbox.lafayette.edu/~reiterc/j/fvj4/index.html Fractals Visualization and J, 4th edition, Part 1] (by Clifford A. Reiter), Chapter 6

[[File:julia_set_java.png|200px|thumb|right]]

{{works with|Java|8}}

import java.awt.image.BufferedImage;

 

public class JuliaSet extends JPanel {

 

public JuliaSet() {

int w = getWidth();

int h = getHeight();

 

for (int x = 0; x < w; x++) {

for (int y = 0; y < h; y++) {

while (zx * zx + zy * zy < 4 && i > 0) {

zx = tmp;

i--;

}

image.setRGB(x, y, c);

}

});

}

 

=={{header|JavaScript}}==

take a look [http://paulo-jorente.de/tests/juliaset/ here].

var maxIterations = 450, minX = -.5, maxX = .5,

minY = -.5, maxY = .5, wid, hei, ctx,

drawFractal();

}

 

=={{header|jq}}==

{{trans|awk}}

 

# $re : -0.8

# $im : 0.156

else .x = .z_real | .y = .z_imag

end )

 

With the above program in a file called julia.jq, the following invocation of jq 1.5 produces the same output as shown in the awk entry on this page:

The following code creates the fractal as a ppm file named julia.ppm. There is no need of an external library to create this image since the ppm format is straightforward to generate.

 

function iter(z,c)

n = 0

writedlm(out, [f(width,height,a,b,i,j) for j = 1:height, i = 1:width], '\n')

end

 

We can then produce a 600x300 ppm image of the Julia set associated to the parameter ''-0.786+0.147i'' as follows.

 

[[File:JuliaSet_julia.png|250px|thumb|right]]

 

@inline function hsv2rgb(h, s, v)

 

 

end

 

end

img[y,x] = RGB{Float64}(r, g, b)

end

end

 

 

=={{header|Kotlin}}==

import java.awt.image.BufferedImage

import javax.swing.JFrame

 

init {

preferredSize = Dimension(800, 600)

}

 

setRenderingHint(RenderingHints.KEY_ANTIALIASING, RenderingHints.VALUE_ANTIALIAS_ON)

while (zx * zx + zy * zy < 4 && i > 0) {

val tmp = zx * zx - zy * zy + cX

i--

}

}

}

}

}

}

 

defaultCloseOperation = JFrame.EXIT_ON_CLOSE

title = "Julia Set"

isResizable = false

pack()

setLocationRelativeTo(null)

isVisible = true

}

 

Mathematica provides built-in functions for Julia sets.

Generate the set of points for the -0.77 +0.22 I Julia set with step sizes of 0.01

Visualize the same Julia set

 

=={{header|Perl}}==

[[File:JuliaSet.perl.png|250px|thumb|right]]

 

my($w, $h, $zoom) = (800, 600, 1);

}

 

 

 

 

}

 

 

=={{header|Processing}}==

float cX = -0.7;

float cY = 0.27015;

float zx, zy;

float maxIter = 300;

 

void draw() {

i -= 1;

}

}

}

noLoop();

}

 

 

 

 

 

=={{header|Python}}==

{{trans|zkl}}

 

 

if __name__ == "__main__":

i -= 1

# convert byte to RGB (3 bytes), kinda magic to get nice colors

 

=={{header|Racket}}==

[[File:JuliaSet_racket.png|250px|thumb|right]]

;; Based on Mandelbrot code (GPL) from:

;; https://github.com/hebr3/Mandelbrot-Set-Racket/blob/master/Mandelbrot.v6.rkt

 

(save-image M-Image "julias.png")

 

=={{header|REXX}}==

:::: which is a

{{trans|COBOL}}

parse arg real imag fine . /*obtain optional arguments from the CL*/

if real=='' | real=="," then real= -0.8 /*Not specified? Then use the default.*/

end /*h*/

if $\='' then say strip($, 'T') /*only display a plot line if non-blank*/

This REXX program makes use of &nbsp; '''scrsize''' &nbsp; REXX program (or BIF) which is used to determine the screen size of the terminal (console),

<br>and the plot size is adjusted according.

(Shown at &nbsp; ^'''1'''/_'''6''' &nbsp; size.)

 

<pre style="font-size:17%">

■■

■■

</pre>

 

=={{header|Ruby}}==

{{trans|AWK}}

puts Complex(c_real, c_imag)

-1.0.step(1.0, 0.04) do |v|

end

 

 

{{out}}

#

</pre>

 

=={{header|Rust}}==

 

use image::{ImageBuffer, Pixel, Rgb};

let _ = img.save("output.png");

 

 

=={{header|Sidef}}==

[[File:JuliaSet_sidef.png|250px|thumb|right]]

 

var (w, h) = (640, 480)

}

 

 

This version generates an ASCII representation:

 

var maxIter = 40

}

print "\n"

{{out}}

##

#######

#######

##

 

</pre>

 

=={{header|zkl}}==

{{trans|Java}}

[[File:JuliaSet.zkl.jpg|250px|thumb|right]]

w,h,zoom:=800,600, 1;

bitmap:=PPM(w,h,0xFF|FF|FF); // White background

 

bitmap.writeJPGFile("juliaSet.jpg",True);