Image:Radiodrome-params-colour.png

From Wikipedia, the free encyclopedia

Image
File history
File links

Radiodrome-params-colour.png‎ (17KB, MIME type: image/png)

This is a file from the Wikimedia Commons. The description on its description page there is shown below.

This Math image (or all images in this article or category) should be recreated using vector graphics as an SVG file. This has several advantages; see Commons:Images for cleanup for more information. If an SVG form of this image is already available, please upload it. After uploading an SVG, replace this template with template {{SupersededSVG|new image name.svg}} in this image.

1 Summary
2 Licensing
3 Used
4 Source Code

[edit] Summary

Description	EN: Curve of pursuit,different parameters, colour - DE: Radiodrome, verschiedene Parameter, Farbe -
Source	self by W!B: - using Euler (Radiodrome_params.en, #SourceCode), StarOffice Draw
Date	created 3. Nov. 2005
Author	W!B:
Permission	Public Domain

[edit] Licensing

I, the author of this work, hereby release it into the public domain. This applies worldwide.

In case this is not legally possible:
I grant anyone the right to use this work for any purpose, without any conditions, unless such conditions are required by law.

[edit] Used

de:Radiodrome

[edit] Source Code

% RADIODROME_PARAMS.en
% SourceCode for EULER by W#B: 29.10.2005, Public Domain
%      EULER is a program for quick and interactive computations with real
%      and complex numbers, matrices, or with intervals.
%      Euler is freeware under the GNU general public license.
%      The Euler homepage is at http://www.rene-grothmann.de
% ==================================================================
>function radiodrome(y, k, d0=1)
## x   - vektor over vektor y, y(i) <> 0!!
## d0 .. normal distance of point of return (1/x'=0)
## k  .. parameter of radiodrome (factor of velocities)
$yy=abs(y)/d0;
$if abs(k)==1; x=sign(k)*((yy^2-1)/4-log(yy)/2);
$ else; k1=1+k; k2=1-k;  x=1/2*((yy^k1-1)/k1-(yy^k2-1)/k2);  endif;
$return x*d0;
$endfunction
>function radiodrome0(k, d0=1)
## x   - value x at y=0
##     if no zeropoint exists, returns x=0
##      (as k=0 will give correct x=0, test for both)
## d0 .. normal distance of point of return (1/x'=0)
$if abs(k)>=1;  x=0;  "radiodrome0: no zeropoint";
$ else; k1=1+k; k2=1-k;  x=k/(k1*k2);  endif;
$return x*d0;
$endfunction
% ------------------------------------------------------------------
>clg; reset;
>xmin=0; xmax=1.6; ymin=0; ymax=1.4; 
>delta0=0.000001; delta1=0.001; delta2=0.02;
>setplot(xmin,xmax,ymin,ymax); keepsquare(1); scaling(0);  
>y=ymin+delta0:delta1:ymax+delta2; x=radiodrome(y, 1);
>plotarea(x,y); hold on; markerstyle("x");
>color(1);  linewidth(1); xplot();
>color(10); linewidth(2); plot(x,y);
>color(12); k=1.4;x=radiodrome(y, k); plot(x,y);
>color(5);  k=0.7;x=radiodrome(y, k); plot(x,y); mark (radiodrome0(k), 0);
>color(14); k=0.5;x=radiodrome(y, k); plot(x,y); mark (radiodrome0(k), 0);
>color(11); k=0.3;x=radiodrome(y, k); plot(x,y); mark (radiodrome0(k), 0);
>"done"
% ------------------------------------------------------------------