clear
normpforig=[];
normdforig=[];
normpf0=[];
normdf0=[];

iterorig=[];
iter0=[];
residualstd=[];
residualmean=[];
dims=[10:10:180];
for n=dims,
genthetaren;
 [normpf, normdf, X, y, theta, iter] = thetarend( adj,n,indsE);
normpforig=[normpforig normpf];
normdforig=[normdforig normdf];
iter0=[iter0 iter];
[normpf, normdf, X, y, theta, iter] = fthetaimp( adj,n,indsE,indsEc,indsEAL);
normpf0=[normpf0 normpf];
normdf0=[normdf0 normdf];
iterorig=[iterorig iter];
save testmc0 iter0 iterorig normpforig normdforig normpf0 normdf0
end

figure(1)
plot(dims,normpforig,'go',dims,normpforig,'y-')
figure(2)
plot(dims,normdforig,'go',dims,normdforig,'y-')
