Archimedean spiral: Difference between revisions
m
syntax highlighting fixup automation
(Added Quackery) |
Thundergnat (talk | contribs) m (syntax highlighting fixup automation) |
||
Line 18:
=={{header|Action!}}==
Action! does not provide trigonometric functions. Therefore a simple implementation for Sin and Cos function has been provided.
<
0 4 9 13 18 22 27 31 36 40 44 49 53 58 62 66 71 75 79 83
88 92 96 100 104 108 112 116 120 124 128 132 136 139 143
Line 68:
DO UNTIL CH#$FF OD
CH=$FF
RETURN</
{{out}}
[https://gitlab.com/amarok8bit/action-rosetta-code/-/raw/master/images/Archimedean_spiral.png Screenshot from Atari 8-bit computer]
Line 74:
=={{header|Ada}}==
{{libheader|SDLAda}}
<
with SDL.Video.Windows.Makers;
Line 145:
Window.Finalize;
SDL.Finalise;
end Archimedean_Spiral;</
=={{header|ALGOL W}}==
{{Trans|AWK}}
This version doubles the characters horiontally to give a slightly more rounded shape.
<
% Translation of AWK which was a trnslation of Applesoft Basic program %
integer procedure max ( integer x, y ) ; begin if x > y then x else y end;
Line 191:
write()
end for_i
end.</
{{out}}
<pre>
Line 240:
Uses Dyalog's [https://sharpplot.com/ SharpPlot] integration, which works on all supported platforms.
<
InitCauseway ⍬ ⍝ initialise current namespace
sp←⎕NEW Causeway.SharpPlot
sp.DrawPolarChart {⍵(360|⍵)}⌽⍳720
View sp</
[https://i.imgur.com/hZDqjjM.png See the plot on imgur.]
Line 250:
=={{header|AutoHotkey}}==
Requires [https://github.com/tariqporter/Gdip GDIP]
<
{
MsgBox, 48, gdiplus error!, Gdiplus failed to start. Please ensure you have gdiplus on your system
Line 300:
Gdip_Shutdown(pToken)
ExitApp
Return</
=={{header|AWK}}==
<
# syntax: GAWK -f ARCHIMEDEAN_SPIRAL.AWK
# converted from Applesoft BASIC
Line 339:
function max(x,y) { return((x > y) ? x : y) }
function min(x,y) { return((x < y) ? x : y) }
</syntaxhighlight>
{{out}}
<pre>
Line 386:
==={{header|AmigaBASIC}}===
{{trans|Locomotive Basic}}
<
b=1.5
pi=3.141592
Line 394:
r=a+b*t
LINE -(320+2*r*SIN(t),100+r*COS(t))
NEXT</
==={{header|Applesoft BASIC}}===
<
120 LET W = H + H / 2
130 HGR2
Line 416:
280 HPLOT X,Y
290 NEXT
</syntaxhighlight>
==={{header|BASIC256}}===
<
# Basic-256 ver 1.1.4
# Archimedean Spiral
Line 452:
imgsave "spiral-Basic-256.png", "PNG"
</syntaxhighlight>
Line 459:
==={{header|Commodore BASIC}}===
Commodore BASIC 2.0 lacks in-built graphics capability. This implementation is written for Commodore BASIC 7.0 that was built into the Commodore 128 computer. Should also work for Commodore BASIC 3.5.
<
2 REM USING COMMODORE BASIC 7.0
3 REM OF THE COMMODORE 128
Line 473:
90 X0 = X : Y0 = Y
100 NEXT T
110 GOTO 110</
==={{header|FreeBASIC}}===
<
' compile with: fbc -s gui
Line 501:
Print : Print "hit any key to end program"
Sleep
End</
==={{header|GW-BASIC}}===
<
20 B = 1
30 SCREEN 1
Line 515:
100 IF INKEY$="" THEN GOTO 100
110 SCREEN 2:SCREEN 0
120 END</
==={{header|IS-BASIC}}===
<
110 OPTION ANGLE DEGREES
120 PLOT 640,360,ANGLE 90;
130 FOR I=2 TO 33.2 STEP .05
140 PLOT FORWARD I,LEFT 5;
150 NEXT</
==={{header|Locomotive Basic}}===
{{trans|Commodore BASIC}}
<
20 mode 2:rad:move 320,200
30 for t=0 to 40*pi step 0.2
Line 533:
50 draw r*sin(t)+320,r*cos(t)+200
60 next
70 while inkey$="":wend</
==={{header|Run BASIC}}===
<
'runs in Run Basic
'Run Basic website http://www.runbasic.com
Line 562:
print "Thank you and Goodbye"
end
End</
==={{header|QBasic}}===
<
WINDOW (-2.67, -2!)-(2.67, 2!)
PI = 4 * ATN(1)
Line 576:
Y = (A + B * T) * SIN(T)
LINE -(X, Y)
NEXT</
==={{header|Sinclair ZX81 BASIC}}===
{{trans|Applesoft BASIC}}
Works with the unexpanded (1k RAM) ZX81. The output is quite blocky, but identifiably a spiral.
<
20 LET B=0.7
30 FOR T=0 TO 7*PI STEP 0.05
40 LET R=A+B*T
50 PLOT R*COS T+32,R*SIN T+22
60 NEXT T</
{{out}}
Screenshot [http://edmundgriffiths.com/zx81archspiral.jpg here].
Line 593:
The BQN online REPL supports some basic plotting functionality through <code>•Plot</code>. This is used to create a spiral plotting function:
<
When called with argument 200, it is similar to the given example diagram.
Line 603:
=={{header|C}}==
Interactive code which asks the parameters a and b as inputs, the number of cycles and the division steps. Requires the [http://www.cs.colorado.edu/~main/bgi/cs1300/ WinBGIm] library.
<syntaxhighlight lang=C>
#include<graphics.h>
#include<stdio.h>
Line 636:
closegraph();
}
</syntaxhighlight>
=={{header|C sharp|C#}}==
<
using System.Linq;
using System.Drawing;
Line 679:
}
}
</syntaxhighlight>
=={{header|C++}}==
[[File:SpiralCpp.png|200px|thumb|right]]
<
#include <windows.h>
#include <string>
Line 794:
spiral s; s.draw( 16, 8 ); return 0;
}
</syntaxhighlight>
=={{header|Clojure}}==
{{Works with| Incanter}}
<
(use '(incanter core stats charts io))
Line 814:
(view (parametric-plot arq-spiral 0 (* 10 Math/PI)))
</syntaxhighlight>
<pre> </pre>
Line 821:
Common Lisp doesn't provide native graphical output. Libraries or bitmapped output could be used instead, but for this solution, the output is accomplished with character printing.
<
(let* ((min-x (apply #'min (mapcar #'car coords)))
(min-y (apply #'min (mapcar #'cdr coords)))
Line 891:
</syntaxhighlight>
=={{header|FOCAL}}==
<
1.2 S B=2
1.3 S N=250
Line 910:
4.1 S X2=R*FSIN(.2*(T+1))
4.2 S Y2=R*FCOS(.2*(T+1))</
This program uses FOCAL-11 on a DEC GT40 vector graphics terminal.
Line 918:
{{Works with|Frege|3.23.888}}
<
import Java.IO
Line 1,006:
drawSpiral g
f <- File.new "SpiralFrege.png"
void $ ImageIO.write buffy "png" f</
Output is [http://funwithsoftware.org/images/2016-SpiralFrege.png here] due to [[User talk:Short Circuit#Is file uploading blocked forever?|Is file uploading blocked forever?]]
Line 1,013:
{{works with|go|1.9}}
Creates a PNG file using only built-in packages.
<
import (
Line 1,060:
log.Fatal(err)
}
}</
=={{header|Haskell}}==
Line 1,067:
{{libheader|Juicy.Pixels}}
{{libheader|Rasterific}}
<
-- stack --resolver lts-7.0 --install-ghc runghc --package Rasterific --package JuicyPixels
Line 1,093:
polyline points
writePng "SpiralHaskell.png" img</
Output is [http://funwithsoftware.org/images/2016-SpiralHaskell.png here] due to [[User talk:Short Circuit#Is file uploading blocked forever?|Is file uploading blocked forever?]]
Line 1,099:
=={{header|J}}==
[[File:Archimedian spiral j.png|200px|thumb|right]]
<
'aspect 1' plot (*^)j.0.01*i.1400</
<div style="clear:both"></div>
Line 1,107:
[[File:archimedian_spiral.png|300px|thumb|right]]
{{works with|Java|8}}
<
import static java.lang.Math.*;
import javax.swing.*;
Line 1,190:
});
}
}</
=={{header|JavaScript}}==
Line 1,196:
{{Works with|Chrome}}
[[File:ASjs.png|200px|right|thumb|Output ASjs.png]]
<
<!-- ArchiSpiral.html -->
<html>
Line 1,221:
}
</script></body></html>
</syntaxhighlight>
{{Output}}
<pre>
Line 1,231:
Assumes the same HTML canvas embedding as above, but is functionally composed.
Defines and logs a set of points, before rendering them to canvas.
<
<head>
<title>Archimedean spiral</title>
Line 1,239:
<h3>Archimedean spiral</h3></p>
<canvas id="spiral" width="640" height="640" style="border: 2px outset;"></canvas>
<script></
<
const
ai = 0.05,
Line 1,272:
Array.from({
length: 1 + n - m
}, (_, i) => m + i);</
<
=={{header|jq}}==
Line 1,279:
'''Works with gojq, the Go implementation of jq'''
====SVG version====
<
def pi: 1 | atan * 4;
Line 1,304:
spiral(0; 10; 0.025)
</syntaxhighlight>
{{out}}
Line 1,312:
====ASCII Art Version====
{{trans|awk}}
<
def min($x;$y): if $x <= $y then $x else $y end;
def max($x;$y): if $x <= $y then $y else $x end;
Line 1,342:
| "\(.)\n" ;
spiral(1; 1; 0.02; 96)</
{{out}}
As for [[#awk|awk]].
Line 1,349:
{{works with|Julia|0.6}}
<
spiral(θ, a=0, b=1) = @. b * θ * cos(θ + a), b * θ * sin(θ + a)
x, y = spiral(1:0.1:10)
println(lineplot(x, y))</
{{out}}
Line 1,379:
=={{header|Kotlin}}==
{{trans|Java}}
<
import java.awt.*
Line 1,454:
f.isVisible = true
}
}</
=={{header|Lua}}==
{{libheader|LÖVE}}
{{works with|LÖVE|11.3}}
<
a=1
b=2
Line 1,480:
end
end
</syntaxhighlight>
=={{header|M2000 Interpreter}}==
<
module Archimedean_spiral {
smooth on ' enable GDI+
Line 1,510:
}
Archimedean_spiral
</syntaxhighlight>
=={{header|Maple}}==
<
plots[polarplot](1+2*theta, theta = 0 .. 6*Pi)
</syntaxhighlight>
=={{header|Mathematica}}/{{header|Wolfram Language}}==
The built-in function PolarPlot easily creates the desired plot
<
=={{header|MATLAB}}==
<
b = 1;
turns = 2;
theta = 0:0.1:2*turns*pi;
polarplot(theta, a + b*theta);</
=={{header|Nim}}==
{{libheader|gintro}}
<
import gintro/[glib, gobject, gtk, gio, cairo]
Line 1,603:
let app = newApplication(Application, "Rosetta.spiral")
discard app.connect("activate", activate)
discard app.run()</
=={{header|PARI/GP}}==
Line 1,611:
[[File:ArchiSpiral2.png|right|thumb|Output ArchiSpiral2.png]]
<
\\ The Archimedean spiral
\\ ArchiSpiral() - Where: lps is a number of loops, c is a direction 0/1
Line 1,631:
ArchiSpiral(640,5,1); \\ArchiSpiral2.png
}
</
{{Output}}
Line 1,644:
=={{header|Perl}}==
{{trans|Raku}}
<
use constant PI => 3.14159265;
Line 1,657:
$img->write(file => 'Archimedean-spiral.png');
</syntaxhighlight>
=={{header|Phix}}==
Line 1,664:
{{libheader|Phix/online}}
You can run this online [http://phix.x10.mx/p2js/Archimedean_spiral.htm here].
<!--<
<span style="color: #000080;font-style:italic;">--
-- demo\rosetta\Archimedean_spiral.exw
Line 1,713:
<span style="color: #000000;">main</span><span style="color: #0000FF;">()</span>
<!--</
=={{header|Processing}}==
Line 1,720:
====with points====
When drawn with points the rotation must be very small, and initially the animation is very slow. This is because the points will move further and further apart as the radius increases.
<
float theta;
float rotation;
Line 1,739:
// check restart
if (x>width/2.0) frameCount=-1;
}</
====with points, rotated====
Rotates the canvas matrix using the built-in rotate() and draws a simple point, rather than computing rotated coordinates with sin()/cos().
<
float rotation;
Line 1,760:
// check restart
if (theta>width/2.0) frameCount=-1;
}</
====with points, vector====
Rotates a vector object of increasing magnitude using the built-in PVector and draws its point, rather than computing rotated coordinates with sin()/cos().
<
float rotation;
Line 1,782:
// check restart
if (pv.mag()>width/2.0) frameCount=-1;
}</
====with line segments====
Draw each new line segments anchored to the previous point in order to keep the spiral visually connected no matter how much the radius expands.
<
float theta;
float rotation;
Line 1,807:
// check restart
if (px>width/2.0) frameCount=-1;
}</
====with line segments, rotated====
Uses the built-in rotate() and screenX() to rotate the frame of reference and then recover the rotated screen position of each next point. Draw each new line segments anchored to the previous point in order to keep the spiral visually connected no matter how much the radius expands.
<
float theta;
float rotation;
Line 1,836:
py = y;
if (theta>width/2.0) frameCount=-1; // start over
}</
==={{header|Processing Python mode}}===
====with points====
When drawn with points the rotation must be very small, and initially the animation is very slow. This is because the points will move further and further apart as the radius increases.
<
rotation = 0.1
Line 1,858:
if x > width / 2.0:
background(255)
theta = 0</
=={{header|PureBasic}}==
<
#XCENTER = 640/2
#YCENTER = 480/2
Line 1,882:
Repeat : Event = WaitWindowEvent() : Until Event = #PB_Event_CloseWindow
EndIf
End</
=={{header|Python}}==
Using the '''turtle''' module.
<
from math import *
color("blue")
Line 1,897:
goto(x, y)
up()
done()</
=={{header|Quackery}}==
<
turtle
0 n->v
Line 1,908:
1 20 v+
1 36 turn ]
2drop</
{{out}}
Line 1,915:
=={{header|R}}==
<
plot((1 + s) * exp(1i * s), type="l"))</
=={{header|Racket}}==
[[File:archemedian-spiral-racket.png]]
<
(require plot
racket/math)
Line 1,942:
;; writes to a file so hopefully, I can post it to RC...
(plot-file (list (archemedian-spiral-renderer2d 0.0 24 4))
"images/archemidian-spiral-racket.png")</
=={{header|Raku}}==
Line 1,948:
{{works with|Rakudo|2018.10}}
<
my ($w, $h) = (400, 400);
Line 1,958:
}
$png.write: 'Archimedean-spiral-perl6.png';</
=={{header|REXX}}==
Line 1,964:
Note: the value of <big><big> ''a'' </big></big> doesn't mean that much as the plot is automatically centered.
<
parse arg cy a b inc chr . /*obtain optional arguments from the CL*/
if cy=='' | cy=="," then cy= 3 /*Not specified? Then use the default.*/
Line 2,005:
if x=pi * .5 then return 1; if x==pi*1.5 then return -1
if abs(x)=pi | x=0 then return 0; q= x*x; z= x
do k=2 by 2 until p=z; p= z; _= -_ *q/(k*k+k); z= z+_; end; return z</
{{out|output|text= when using the following inputs: <tt> 13 , 5 , db </tt>}}
Line 2,205:
=={{header|Ring}}==
<
/*
+---------------------------------------------------------------------------------------------------------
Line 2,310:
return
</syntaxhighlight>
=={{header|Ruby}}==
Line 2,316:
{{libheader|JRubyArt}}
JRubyArt is an implementation of Processing in ruby, that uses JRuby to provide the interoperability with the java libraries.
<
INCR = 0.1
attr_reader :x, :theta
Line 2,337:
size(300, 300)
end
</syntaxhighlight>
=={{header|Rust}}==
<
extern crate bmp;
Line 2,371:
// Save the image
let _ = img.save("archimedean_spiral.bmp").unwrap_or_else(|e| panic!("Failed to save: {}", e));
}</
=={{header|SAS}}==
<
h=constant('pi')/40;
do i=0 to 400;
Line 2,387:
proc sgplot;
series x=x y=y;
run;</
=={{header|Scala}}==
===Java Swing Interoperability===
<
object ArchimedeanSpiral extends App {
Line 2,458:
)
}</
=={{header|Scheme}}==
{{libheader|Scheme/PsTk}}
<
(import (scheme base)
(scheme complex)
Line 2,498:
(draw-spiral canvas))
(tk-event-loop tk))
</syntaxhighlight>
=={{header|Scilab}}==
Line 2,518:
scf(2);
plot2d(x,y);</
=={{header|Seed7}}==
<
include "draw.s7i";
include "keybd.s7i";
Line 2,548:
DRAW_FLUSH;
ignore(getc(KEYBOARD));
end func;</
=={{header|Sidef}}==
{{trans|Raku}}
<
define π = Num.pi
Line 2,566:
}
img.write(file => 'Archimedean_spiral.png')</
Output image: [https://github.com/trizen/rc/blob/master/img/archimedean-spiral-sidef.png Archimedean spiral]
=={{header|Stata}}==
<
scalar h=_pi/40
set obs 400
Line 2,576:
gen x=(1+t)*cos(t)
gen y=(1+t)*sin(t)
line y x</
=={{header|Tcl}}==
This creates a little Tk GUI where you can interactively enter values for `a` and `b`. The spiral will be re-drawn automatically thanks to `trace`:
<
# create widgets
Line 2,645:
update ;# lay out widgets before trying to draw
draw
vwait forever ;# go into event loop until window is closed</
=={{header|VBA}}==
<
Dim chrt As Chart
Set chrt = ActiveSheet.Shapes.AddChart.Chart
Line 2,670:
Next i
plot_coordinate_pairs x, y
End Sub</
=={{header|Wren}}==
{{trans|Sidef}}
{{libheader|DOME}}
<
import "dome" for Window
Line 2,702:
static draw(dt) {}
}</
=={{header|XPL0}}==
Looks a lot like the C++ image.
<
[SetVid($12); \set 640x480 graphics
A:= 0.0; B:= 3.0; T:= 0.0;
Line 2,715:
T:= T + 0.03; \increase angle (Theta)
until T >= 314.159; \50 revs
]</
=={{header|Yabasic}}==
{{trans|Sinclair_ZX81_BASIC}}
<
10 LET A=1.5
20 LET B=0.7
Line 2,725:
40 LET R=A+B*T
50 LINE TO R*COS(T),R*SIN(T)
60 NEXT T</
=={{header|zkl}}==
[[File:ArchimedeanSpiral.zk.jpg|250px|thumb|right]]
Uses the PPM class from http://rosettacode.org/wiki/Bitmap/Bresenham%27s_line_algorithm#zkl
<
w,h:=640,640; centerX,centerY:=w/2,h/2;
bitmap:=PPM(w+1,h+1,0xFF|FF|FF); // White background
Line 2,741:
}
bitmap.writeJPGFile("archimedeanSpiral.jpg");
}(0,5,7);</
|