%@AUTEUR:Guillaume Connan
prologues:=2;

verbatimtex
%&latex
\documentclass{article}
\begin{document}
etex

input courbes;
input geo;

color vert_e, turquoise, orange, vert_fonce, rose, vert_mer, bleu_ciel, or, rouge_v,bleu_m,bleu,bleu_f;
vert_e:=(0,0.790002,0.340007);
turquoise:=(0.250999,0.878399,0.815699);
orange:=(0.589999,0.269997,0.080004);
vert_fonce:=(0,1.4*0.392193,0);
rose:=(1.0, 0.752907, 0.796106);
bleu_ciel:=(1.2*0.529405,1.2*0.807794,1);%.2*0.921598);
or:=(1,0.843104,0);
rouge_v:=(0.829997,0.099994,0.119999);
bleu_m:=(0.7*0.529405,0.7*0.807794,0.7);%*0.921598);
bleu_f:=(0.211762,0.3231176,0.3686392);
bleu:=(0.529405,0.807794,1);


vardef f(expr x) =1/(sqrt x) enddef;

vardef trace (suffix g)(expr a,b,inc) =
  save i; numeric i;
  for i=a step inc until b:
    (i*ux,f(i)*uy) ..
  endfor (b*ux,f(b)*uy)
enddef;

% lim infinie en infini : Energie cinétique

% Déclarations des constantes %
numeric xmin, xmax, ymin, ymax;

ux:=1cm; uy:=0.9cm; xmin := -0.5; xmax := 8;
ymin := xmin; ymax := xmax-4;

% Définitions des axes et labels associés
vardef axes =
  drawarrow (ux*xmin,0) -- (ux*xmax,0); % axe des x
  drawarrow (0,uy*ymin) -- (0,uy*ymax); % axe des y
  label.lrt(btex $v$ etex,(xmax*ux,0)); % label de l'axe des x
  label.lft(btex $\mathcal{E}_c$ etex,(0,ymax*uy)); % label de l'axe des y
enddef;


beginfig(1);

  label.llft(btex $0$ etex,(0,0));
  
%Fonction, graphe
  vardef f(expr x) =.2*x*x enddef;
  
  vardef tracee (suffix g)(expr a,b,inc) =
    save i; numeric i;
    for i=a step inc until b:
      (i*ux,f(i)*uy) ..
    endfor (b*ux,f(b)*uy)
  enddef;
  
%Points sur le graphe
  pair A[],M[],O[];
  path z;
  u:=.7;%facteur d'échelle
  aa:=2;bb:=5;cc:=9;%numeros des points particuliers
  
  for i=1 upto 12: A[i]:=(i*u*ux,0);% A comme abscisse
    M[i]:=(i*u*ux,f(i*u)*uy); %M comme courbe !!
    O[i]:=(0,f(i*u)*uy); %O comme ordonnée
  endfor ;
  
  z:=O[bb]--O[bb]shifted (xmax*ux,0)--O[bb]shifted (xmax*ux,0)shifted(0,(ymax-f(5*u))*uy)--O[bb]shifted(0,(ymax-f(5*u))*uy)--cycle;
  fill z withcolor bleu_ciel;
  axes;

  draw A[bb]--M[bb]--O[bb] dashed evenly withcolor bleu_m;

  label.bot(btex $v_0$ etex,A[bb]);

%Tracé des droites
  
  pair C,D;
  
  C=O[bb]; %Point seuil
  D=1/3(O[cc]+2O[cc-1]); %Point seuil
  
  pickup pencircle scaled 2bp;
  
  draw C--C  shifted (xmax*ux,0) withcolor bleu_m;
  
  label.lft(btex $\mathcal{E}_{\mbox{seuil}}$ etex,C shifted (0,-.1*uy));
  
  draw tracee(g,0,4.472,.1)withcolor bleu_f;
endfig;
end