Od[F_, x0_,n_,d_,cmin_, cmax_, cstep_,xmin_, xmax_,  ptsize_] := Module[{ mylist, c, x, i },
	mylist = {};
	For[c = cmin, c <= cmax, c = c + cstep,
		x = x0;
		For[i = 1, i < d, i = i + 1,
				x = F[x, c]];
		For[i = 1, i < n, i = i + 1,
				x = F[x, c];
				mylist = Append[mylist, {c, x}];
				];
			];
	ListPlot[mylist, PlotStyle->PointSize[ptsize], PlotRange->{xmin, xmax}]
			];

f[x_, k_] := k x (1 - x);

Od[f, 0.5,100,100,0.0, 4.0, 0.01,0, 1, 0.002]

Od[f, 0.5,200,100,2.8, 4.0, 0.005,0, 1, 0.002]

Od[f, 0.5,200,100,3.8, 3.875, 0.00025,0, 1, 0.002]

Od[#2 Sin[#1]&, N[Pi/2],100,100,1.0, N[Pi], 0.005,0, Pi, 0.002]

Od[#2 (#1-#1^3/3)&, 1.0,100,100,1.0, N[3Sqrt[3]/2], 0.005,0, Sqrt[3], 0.002]

Od[#2 Cos[#1]&, 1.0,100,100,-2.96,0.0, 0.02,-Pi, Pi, 0.006]

Od[#2 Cos[#1]&, 1.0,100,100,-2.96,0.0, 0.005,-Pi, Pi, 0.002]

g[x_, c_] :=x^2 + c;

Od[g,0,100,100,-2.0, 0.25, 0.01,-2, 2, 0.002]

Od[g,0,200,100,-1.5, -1.3, 0.001,-1.5, -0.9, 0.002]