hold off
t0=clock;
k1=.03;g=50;n3=8;
a=(1+k1*(g-1))/g;
n1=1000;
m=(1-k1)*(1-a);
n=a*(1-k1);
b=k1*(1-a);
p1=a^(1/(1-k1))
p2=p1^2*(m+b)/n;
v=[p1 p2];x=zeros(1,n1-2);
p=[v x];w=zeros(1,n1);w(n3)=a^(n3/(1-k1));q=n3/((n3+1)*n);
v1=a^(n3/(1-k1));s=a^(k1/(1-k1));t=(1-a)/(1-k1);v2=n3*t*v1*s;v3=.5*v2*t*s*(n3+1+k1);
v4=[v1 v2 v3];
for i=3:n1,
suma=0;
k=fix(i/2);
l=1;
while (l<k)
suma=suma+p(l)*p(i-l);
l=l+1;
end
if i-2*k>0
	suma=suma+p(k)*p(i-k);
	p(i)=suma*(2*m+i*b)/((i-1)*n);
else
	suma=suma+p(k)^2/2;
	p(i)=suma*(2*m+i*b)/((i-1)*n);
end
end
duratad=etime(clock,t0)
resd=p;t1=clock;
%stem(resd)
x=(1:n1);
%semilogy(x,resd)
loglog(x,resd,'b')
xlabel('g');
ylabel('P(m/1) si P(m/8)');
hold on
pack
%sigma=0;
for l=n3+1:n1
sigma=0;%w(n3)=a^(n3/(1-k1));
	for j=1:l-n3
	sigma=sigma+p(j)*w(l-j);
	end
  w(l)=sigma*q*(l*b+m*(n3+1))/(l-n3);
end
durataf=etime(clock,t1);      
loglog(x,w,'r')	
duratag=etime(clock,t1);
cc=conv(p,p);ee=conv(cc,cc);ff=conv(ee,ee); dd=w(1,8:200)-ff(1,1:193);
duratae=etime(clock,t0);
mc=w(1,n3:n3+15);