Difference between revisions of "File:Susinploat300.jpg"

From TORI
Jump to navigation Jump to search
(Importing image file)
 
 
Line 1: Line 1:
  +
[[Explicit plot]] of function [[SuSin]] in comparison to its approximations.
−
Importing image file
  
  +
  +
Black think solid curve shows $y= \mathrm{SuSin}(x)$, this is [[superfunction]] of [[sin]], id est, super sin.
  +
  +
The upper thin blue curves shows the leading term of its asymptotic expansion, $y=\sqrt{3/x}$, suggested in 2012 by [[Kursernas Hemsidor]] <ref>
  +
http://web.abo.fi/fak/mnf/mate/kurser/dynsyst/Iteration%20of%20sin(x).pdf Kursernas Hemsidor. 273027 Introduction to Dynamical Systems 2012. Derivation of Niklas Carlsson;s formula. Cited by the state for December 2013: Let function f(x) be sin(x). We want to evaluate, approximately, the value of the nth iterate of f(x) ... If the formula is correct, it will take 3⋅1010 .. steps to reat 0.00001 from 1. ..
  +
$f^n(x)\approx \sqrt{\frac{3}{n}}$, $x\approx 1$, $n$ lage ... How good is this formula? after 60000 iterations of the value 0.0071 and the discrepancy -4.7s-007 ..
  +
</ref>.
  +
  +
The red dashed curve shows the approximation suggested in 2012 by [[Thomas Curtright]] <ref> http://server.physics.miami.edu/~curtright/Schroeder.html
  +
Dr. Thomas Curtright. Continuous iterates continue to be interesting, after 150 years of study. (2012)
  +
<i> As a first illustration, we display the continuous iterates of the sine function, sin[t](x). Note that the maximum values at x = π/2 are approximately given by exp[(1-√t) ln(π/2)].
  +
</i></ref>
  +
  +
$y=\exp\Big( \big(1\!-\!\sqrt{x}\big)\ln(\pi/2)\Big)$
  +
  +
The lowest red thin curve shows the difference between the $SuSin(x)$ and the approximation by Thomas, scaled with factor 10.
  +
  +
==References==
  +
<references/>
  +
  +
==[[C++]] generator of curves==
  +
//Files [[ado.cin]], [[arcsin.cin]], and [[susin.cin]] should be loaded to working directory in order to compile the code below
  +
  +
//<poem><nomathjax><nowiki>
  +
#include <math.h>
  +
#include <stdio.h>
  +
#include <stdlib.h>
  +
#define DB double
  +
#define DO(x,y) for(x=0;x<y;x++)
  +
using namespace std;
  +
#include<complex>
  +
typedef complex<double> z_type;
  +
#define Re(x) x.real()
  +
#define Im(x) x.imag()
  +
#define I z_type(0.,1.)
  +
  +
#include "ado.cin"
  +
#include "arcsin.cin"
  +
#include "susin.cin"
  +
  +
int main(){ int j,k,m,n; DB x,y, p,q, t; z_type z,c,d;
  +
DB x0=0.;
  +
DO(m,14){y=Re(susin(z_type(1.,1.e-9)+x0))-1.;
  +
x0+=4.*y;
  +
printf("%2d %19.16f %19.16f\n",m,x0,y);}
  +
//FILE *o;o=fopen("susinplot1.eps","w"); ado(o,1002,244);
  +
//FILE *o;o=fopen("04.eps","w"); ado(o,1002,348);
  +
FILE *o;o=fopen("susinploa.eps","w"); ado(o,1002,348);
  +
#define M(x,y) {fprintf(o,"%6.4f %6.4f M\n",0.+x,0.+y);}
  +
#define L(x,y) {fprintf(o,"%6.4f %6.4f L\n",0.+x,0.+y);}
  +
fprintf(o,"1 106 translate\n 100 100 scale\n");
  +
fprintf(o,"1 setlinejoin 2 setlinecap\n");
  +
for(m=0;m<11;m++){M(m,-1) L(m,2) }
  +
for(n=-1;n<3;n++){M( 0,n) L(10,n)}
  +
fprintf(o,".006 W 0 0 0 RGB S\n");
  +
M(0,M_PI/2.); L(10,M_PI/2)
  +
fprintf(o,".004 W 0 0 0 RGB S\n");
  +
fprintf(o,"1 setlinejoin 1 setlinecap\n");
  +
DO(m,100){ x=.5+.1*m; y=sqrt(3./x); if(m==0) M(x,y) else L(x,y) ; if ( x>10.) break;}
  +
fprintf(o,".006 W 0 0 1 RGB S\n");
  +
//M(0,M_PI/2.);
  +
fprintf(o,"1 setlinejoin 0 setlinecap\n");
  +
DO(m,300){ x=.0001+.04*m/(1+5./(.3+m)); y=exp((1.-sqrt(x))*log(M_PI/2)); if(m/2*2==m) M(x,y) else L(x,y) ; if ( x>10.) break;}
  +
fprintf(o,".02 W 1 0 0 RGB S\n");
  +
fprintf(o,"1 setlinejoin 2 setlinecap\n");
  +
M(0,M_PI/2.);
  +
DO(m,2002){ x=.005*(m+.3); z=z_type(x,1.e-8); c=susin(z); y=Re(c); L(x,y); printf("%8.5f %8.5f\n",x,y); }
  +
fprintf(o,".012 W 0 0 0 RGB S\n");
  +
M(0,0);
  +
DO(m,20022){ x=.005*(m+.3); z=z_type(x,1.e-8); c=susin(z);
  +
y=exp((1.-sqrt(x))*log(M_PI/2))- Re(c); y*=10; L(x,y); printf("%8.5f %8.5f\n",x,y);
  +
if(x>10) break; }
  +
fprintf(o,".006 W 1 0 0 RGB S\n");
  +
  +
fprintf(o,"showpage\n");
  +
fprintf(o,"%c%cTrailer\n",'%','%');
  +
fclose(o);
  +
system("epstopdf susinploa.eps");
  +
system( "open susinploa.pdf"); //for macintosh
  +
getchar(); system("killall Preview"); // For macintosh
  +
}
  +
</nowiki></nomathjax></poem>
  +
  +
==[[Latex]] generator of labels==
  +
  +
<poem><nomathjax><nowiki>
  +
\documentclass[12pt]{article}
  +
\usepackage{geometry}
  +
\usepackage{graphics}
  +
\usepackage{rotating}
  +
\newcommand \rot {\begin{rotate}}
  +
\newcommand \ero {\end{rotate}}
  +
\paperwidth 1026pt
  +
\paperheight 345pt
  +
\topmargin -109pt
  +
\oddsidemargin -90pt
  +
\newcommand \sx {\scalebox}
  +
\pagestyle{empty}
  +
\begin{document}
  +
\begin{picture}(1016,328)
  +
%\put(20,1){\includegraphics{susinplot1}}
  +
\put(20,1){\includegraphics{susinploa}}
  +
%\put(20,1){\includegraphics{04}}
  +
\put(2,312){\sx{2.4}{$y$}}
  +
\put(-1,257){\sx{2.8}{$\frac{\pi}{2}$}}
  +
\put(2,198){\sx{2.4}{$1$}}
  +
\put(2,98){\sx{2.4}{$0$}}
  +
\put(15,-16){\sx{2.4}{$0$}}
  +
\put(115,-16){\sx{2.4}{$1$}}
  +
\put(215,-16){\sx{2.4}{$2$}}
  +
\put(315,-16){\sx{2.4}{$3$}}
  +
\put(415,-16){\sx{2.4}{$4$}}
  +
\put(516,-16){\sx{2.4}{$5$}}
  +
\put(616,-16){\sx{2.4}{$6$}}
  +
\put(717,-16){\sx{2.4}{$7$}}
  +
\put(817,-16){\sx{2.4}{$8$}}
  +
\put(917,-16){\sx{2.4}{$9$}}
  +
\put(1010,-16){\sx{2.5}{$x$}}
  +
%\put(45,134){\sx{2.5}{$y\!=\!\mathrm{SuSin}(x)$}}
  +
\put(190,246){\sx{1.8}{\rot{-12}$y\!=\! \sqrt{3/x}$\ero}}
  +
\put(190,201){\sx{1.8}{\rot{-6}$y\!=\!\mathrm{SuSin}(x)$\ero}}
  +
\put(190,176){\sx{1.8}{\rot{-6}$y\!=\!\exp((1\!-\!\sqrt{x})\ln(\pi/2))$\ero}}
  +
  +
\put(190,110){\sx{1.8}{\rot{-8}$y\!=\! 10\Big(\exp((1\!-\!\sqrt{x})\ln(\pi/2))-\mathrm{SuSin}(x)\Big)$\ero}}
  +
\end{picture}
  +
\end{document}
  +
</nowiki></nomathjax></poem>
  +
  +
[[Category:Explicit plot]]
  +
[[Category:SuSin]]
  +
[[Category:Superfunction]]
  +
[[Category:Super sin]]
  +
[[Category:C++]]
  +
[[Category:Latex]]
  +
[[Category:Book]]

Latest revision as of 08:53, 1 December 2018

Explicit plot of function SuSin in comparison to its approximations.

Black think solid curve shows $y= \mathrm{SuSin}(x)$, this is superfunction of sin, id est, super sin.

The upper thin blue curves shows the leading term of its asymptotic expansion, $y=\sqrt{3/x}$, suggested in 2012 by Kursernas Hemsidor [1].

The red dashed curve shows the approximation suggested in 2012 by Thomas Curtright [2]

$y=\exp\Big( \big(1\!-\!\sqrt{x}\big)\ln(\pi/2)\Big)$

The lowest red thin curve shows the difference between the $SuSin(x)$ and the approximation by Thomas, scaled with factor 10.

References

  1. ↑ http://web.abo.fi/fak/mnf/mate/kurser/dynsyst/Iteration%20of%20sin(x).pdf Kursernas Hemsidor. 273027 Introduction to Dynamical Systems 2012. Derivation of Niklas Carlsson;s formula. Cited by the state for December 2013: Let function f(x) be sin(x). We want to evaluate, approximately, the value of the nth iterate of f(x) ... If the formula is correct, it will take 3⋅1010 .. steps to reat 0.00001 from 1. .. $f^n(x)\approx \sqrt{\frac{3}{n}}$, $x\approx 1$, $n$ lage ... How good is this formula? after 60000 iterations of the value 0.0071 and the discrepancy -4.7s-007 ..
  2. ↑ http://server.physics.miami.edu/~curtright/Schroeder.html Dr. Thomas Curtright. Continuous iterates continue to be interesting, after 150 years of study. (2012) As a first illustration, we display the continuous iterates of the sine function, sin[t](x). Note that the maximum values at x = π/2 are approximately given by exp[(1-√t) ln(π/2)].

C++ generator of curves

//Files ado.cin, arcsin.cin, and susin.cin should be loaded to working directory in order to compile the code below

//


#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#define DB double
#define DO(x,y) for(x=0;x<y;x++)
using namespace std;
#include<complex>
typedef complex<double> z_type;
#define Re(x) x.real()
#define Im(x) x.imag()
#define I z_type(0.,1.)

#include "ado.cin"
#include "arcsin.cin"
#include "susin.cin"

int main(){ int j,k,m,n; DB x,y, p,q, t; z_type z,c,d;
DB x0=0.;
DO(m,14){y=Re(susin(z_type(1.,1.e-9)+x0))-1.;
        x0+=4.*y;
        printf("%2d %19.16f %19.16f\n",m,x0,y);}
//FILE *o;o=fopen("susinplot1.eps","w"); ado(o,1002,244);
//FILE *o;o=fopen("04.eps","w"); ado(o,1002,348);
FILE *o;o=fopen("susinploa.eps","w"); ado(o,1002,348);
#define M(x,y) {fprintf(o,"%6.4f %6.4f M\n",0.+x,0.+y);}
#define L(x,y) {fprintf(o,"%6.4f %6.4f L\n",0.+x,0.+y);}
fprintf(o,"1 106 translate\n 100 100 scale\n");
fprintf(o,"1 setlinejoin 2 setlinecap\n");
for(m=0;m<11;m++){M(m,-1) L(m,2) }
for(n=-1;n<3;n++){M( 0,n) L(10,n)}
fprintf(o,".006 W 0 0 0 RGB S\n");
M(0,M_PI/2.); L(10,M_PI/2)
fprintf(o,".004 W 0 0 0 RGB S\n");
fprintf(o,"1 setlinejoin 1 setlinecap\n");
DO(m,100){ x=.5+.1*m; y=sqrt(3./x); if(m==0) M(x,y) else L(x,y) ; if ( x>10.) break;}
fprintf(o,".006 W 0 0 1 RGB S\n");
//M(0,M_PI/2.);
fprintf(o,"1 setlinejoin 0 setlinecap\n");
DO(m,300){ x=.0001+.04*m/(1+5./(.3+m)); y=exp((1.-sqrt(x))*log(M_PI/2)); if(m/2*2==m) M(x,y) else L(x,y) ; if ( x>10.) break;}
fprintf(o,".02 W 1 0 0 RGB S\n");
fprintf(o,"1 setlinejoin 2 setlinecap\n");
M(0,M_PI/2.);
DO(m,2002){ x=.005*(m+.3); z=z_type(x,1.e-8); c=susin(z); y=Re(c); L(x,y); printf("%8.5f %8.5f\n",x,y); }
fprintf(o,".012 W 0 0 0 RGB S\n");
M(0,0);
DO(m,20022){ x=.005*(m+.3); z=z_type(x,1.e-8); c=susin(z);
            y=exp((1.-sqrt(x))*log(M_PI/2))- Re(c); y*=10; L(x,y); printf("%8.5f %8.5f\n",x,y);
if(x>10) break; }
fprintf(o,".006 W 1 0 0 RGB S\n");

fprintf(o,"showpage\n");
fprintf(o,"%c%cTrailer\n",'%','%');
fclose(o);
      system("epstopdf susinploa.eps");
      system( "open susinploa.pdf"); //for macintosh
      getchar(); system("killall Preview"); // For macintosh
}

Latex generator of labels


\documentclass[12pt]{article}
\usepackage{geometry}
\usepackage{graphics}
\usepackage{rotating}
\newcommand \rot {\begin{rotate}}
\newcommand \ero {\end{rotate}}
\paperwidth 1026pt
\paperheight 345pt
\topmargin -109pt
\oddsidemargin -90pt
\newcommand \sx {\scalebox}
\pagestyle{empty}
\begin{document}
\begin{picture}(1016,328)
%\put(20,1){\includegraphics{susinplot1}}
\put(20,1){\includegraphics{susinploa}}
%\put(20,1){\includegraphics{04}}
\put(2,312){\sx{2.4}{$y$}}
\put(-1,257){\sx{2.8}{$\frac{\pi}{2}$}}
\put(2,198){\sx{2.4}{$1$}}
\put(2,98){\sx{2.4}{$0$}}
\put(15,-16){\sx{2.4}{$0$}}
\put(115,-16){\sx{2.4}{$1$}}
\put(215,-16){\sx{2.4}{$2$}}
\put(315,-16){\sx{2.4}{$3$}}
\put(415,-16){\sx{2.4}{$4$}}
\put(516,-16){\sx{2.4}{$5$}}
\put(616,-16){\sx{2.4}{$6$}}
\put(717,-16){\sx{2.4}{$7$}}
\put(817,-16){\sx{2.4}{$8$}}
\put(917,-16){\sx{2.4}{$9$}}
\put(1010,-16){\sx{2.5}{$x$}}
%\put(45,134){\sx{2.5}{$y\!=\!\mathrm{SuSin}(x)$}}
\put(190,246){\sx{1.8}{\rot{-12}$y\!=\! \sqrt{3/x}$\ero}}
\put(190,201){\sx{1.8}{\rot{-6}$y\!=\!\mathrm{SuSin}(x)$\ero}}
\put(190,176){\sx{1.8}{\rot{-6}$y\!=\!\exp((1\!-\!\sqrt{x})\ln(\pi/2))$\ero}}

\put(190,110){\sx{1.8}{\rot{-8}$y\!=\! 10\Big(\exp((1\!-\!\sqrt{x})\ln(\pi/2))-\mathrm{SuSin}(x)\Big)$\ero}}
\end{picture}
\end{document}

File history

Click on a date/time to view the file as it appeared at that time.

Date/TimeThumbnailDimensionsUserComment
current06:14, 1 December 2018Thumbnail for version as of 06:14, 1 December 20184,259 × 1,432 (460 KB)Maintenance script (talk | contribs)Importing image file

You cannot overwrite this file.

The following page uses this file:

Metadata

This file contains additional information, probably added from the digital camera or scanner used to create or digitize it.

If the file has been modified from its original state, some details may not fully reflect the modified file.

Retrieved from "https://mizugadro.mydns.jp/t/index.php?title=File:Susinploat300.jpg&oldid=21747"