evaltime.m
View Options

			N = {}
			nstart = 7;
			nstop = 30;
			step = 1;

			ind = nstart;
			while ind <= nstop
			N = [N, {ind}];
			ind = ind + step;
			end

			%%%%%%%%%%%%%%%%% Compute

			W = length(N);

			disp(sprintf('Starting symulation of %d values', W))
			result = zeros(W,2);
			for i=1:W
			n = N{i};
			disp(sprintf('Evaluating value: %d', n))
			result(i,1) = n;
			tic;
			plotherm(n, {-6, 6}, 0);
			result(i,2) = toc;
			disp(sprintf('Elapsed time: %f\n', result(i,2)))
			end

			%%%%%%%%%%%%%%%%%%% Plotting

			figure
			grid on
			xlabel('n')
			ylabel('Elapsed time')
			semilogy(result(:,1), result(:,2))

			out = polyfit(result(:,1), log(result(:,2)), 1);

			disp(sprintf('t = %g * exp(%g * n)', exp(out(2)), out(1)))

compphys/exo1/evaltime.m33bba04159d6master