Ця математична ілюстрація має бути конвертована у векторний формат SVG. Це дає декілька переваг, докладніше про які Ви можете прочитати на сторінці Commons:Media for cleanup. Якщо Вам вже зараз доступна векторна версія даного зображення, завантажте її, будь ласка. Після завантаження замініть цей шаблон на такий: {{vector version available|Назва_завантаженого_файлу.svg}}.

Опис файлу

Ліцензування

Я, власник авторських прав на цю роботу, передаю роботу в суспільне надбання. Застосовується по всьому світу. У деяких країнах це не може бути юридично можливо, в такому випадку: Я даю кожному право на використання цієї роботи для будь-яких цілей, без будь-яких умов, якщо такі умови не вимагаються за законом.

Used

• de:Radiodrome

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"
% ------------------------------------------------------------------