\documentclass[]{standalone}
\usepackage{tikz}
\usepackage{pgfplots}
\usetikzlibrary{calc} 
\pgfplotsset{compat=1.14} 

\usepackage{xcolor}

%Radius of rolling circle
\pgfmathparse{1} \let\r\pgfmathresult
   
%Number of curves to draw.
\pgfmathparse{1} \let\N\pgfmathresult
   
%Relative position of rollong point. Lambda=\pm 1: cycloid
\pgfmathparse{2.0} \let\LambdaMin\pgfmathresult
\pgfmathparse{2.0} \let\LambdaMax\pgfmathresult

%Steps of aLambda
\pgfmathparse{(\LambdaMax-\LambdaMin)/\N} \let\dLambda\pgfmathresult
  
%Variation of colors - color in [0,100].
\pgfmathparse{100/\N} \let\dcolor\pgfmathresult

%Fixed time - draw rolling
\pgfmathparse{4/3} \let\T\pgfmathresult

%Epsilon (not too small!)
\pgfmathparse{1} \let\Eps\pgfmathresult


\input{Curve.Styles.tikz.tex}
\input{Curve.tikz.tex}
\input{PQ.tikz.tex}
\input{Rolling.tikz.tex}

\input{Trochoid.tikz.tex}
\input{Trochoid.Rolling.tikz.tex}

\begin{document}
   %Multiply \T with pi, \T is t_0
   \pgfmathparse{pi*\T} \let\T\pgfmathresult
   %Second parametervalue, t
   \pgfmathparse{\T+\Eps} \let\TT\pgfmathresult

   
   \begin{tikzpicture}
   \begin{scope}   
       \clip (-\r-1,-2*\r) rectangle (5*pi,6*\r);
       \foreach \n in {0,1,...,\N}
       {
           %\Lambda and color, \color
           \pgfmathparse{\LambdaMin+\n*\dLambda} \let\Lambda\pgfmathresult
           \pgfmathparse{\n*\dcolor} \let\col\pgfmathresult
          
           \DrawCurve{\Lambda}{color=blue!\col!cyan,thin} 
       }
       
       %Draw cycloid
       \pgfmathparse{1} \let\Lambda\pgfmathresult
       \DrawCurve{\Lambda}{color=blue,thick}
       \DrawRolling{\T}
       
       \CalcR{\Lambda}{\T}
       \filldraw[rolling] (R) circle(1pt);
       \draw[rolling,-latex] (Rolling) -- (R);
            
       \CalcR{\LambdaMax}{\T}
       
       \coordinate (Frenet) at (0,0);
       \DrawFrenet{\LambdaMax}{\T}{(Frenet)}
       
       \filldraw[curves] (R) circle(1pt) node[right] {$t_0$};
      
       %Tangent      
       \CalcDr{\LambdaMax}{\T}
       %\draw[-latex,curves] (R)-- +(dR);
       
       \DrawTaylorII{\LambdaMax}{\T}{\TT}
      
       %Coordinate system
       \draw[-latex]
          (-\r,0) 
          --
          (4*pi+0.5*\r*\LambdaMax,0) node [right]
          {$x$};
       \draw[-latex]
          (0,-\r)
          --
          (0,1.5*\r*\LambdaMax)
          node [above] {$y$};
              
   \end{scope}   
   \end{tikzpicture}
\end{document}