Feig:=proc(R) local k,r,tmp,ligne; tmp:=floor(R/0.01); res:=[0$tmp]; ligne:=[0$20]; for r from 0 to R by 1/100 do tmp:=(0.2); for k from 0 to 100 do tmp:=r*tmp*(1-tmp); od;; for k from 0 to 20 do ligne[k+1]:=point(r,tmp); tmp:=r*tmp*(1-tmp); od;; array_sto(ligne,res[r*100+1]); od;; RETURN(res); end:;