%% numerik 1, aufgabe 4, Bibiana Gluski, Marina Kasimir, Thomas Dettbarn
un_a=3*sqrt(3); 		%% zunaechst einmal die startwerte festlegen
un_b=3*sqrt(3);			%% den auch
Un_a=6*sqrt(3);
Un_b=6*sqrt(3);
n=3;

%% jetzt werden die formeln aus dem aufgabenzettel angewendet, ganz ganz
%% oft, die _a sind die oberen, die _b sind die unteren.
erg=fopen('Pi-Berechnung.txt','w+');
fprintf(erg,'Ergebnisse der pi-Berechnung\n');
fprintf(erg,'	n	un_a		u2n_a		U2n_a		Un_a			un_b		u2n_b		U2n_b		Un_b\n');
for i=1:10
	u2n_a=sqrt(4*n*(2*n-sqrt(4*n*n-un_a*un_a)));
	u2n_b=un_b*sqrt(4*n/(2*n+sqrt(4*n*n-un_b*un_b)));

	U2n_a=4*n*(sqrt(4*n*n+Un_a*Un_a)-2*n)/Un_a;
	U2n_b=4*n*Un_b/(sqrt(4*n*n+Un_b*Un_b)+2*n);	
%%	y=[n;un_a;u2n_a;U2n_a;Un_a;un_b,u2n_b,U2n_b,Un_b];
	fprintf(erg,'	%i	%2.10f	%2.10f	%2.10f	%2.10f		%2.10f	%2.10f	%2.10f	%2.10f\n',n,un_a,u2n_a,U2n_a,Un_a,un_b,u2n_b,U2n_b,Un_b);		%% einmal ausgeben

	n=n*2;
	un_a=u2n_a;		%% im naechsten schritt ist un ja u2n
	un_b=u2n_b;		
	Un_a=U2n_a;
	Un_b=U2n_b;
end
fclose(erg);		%% und fertig. schluss.
