f(r,x):=(1-r^2)/(1-2*r*cos(x)+r^2);
valores:[0.1,0.2,0.3,0.9];
funciones:makelist(f(valores[k],x),k,1,length(valores));
plot2d(funciones,[x,-%pi,%pi]);