lepsi vypisovani, funkce.h

ukol-1/Makefile

@@ -1,6 +1,6 @@
 program=integraly
 
-HEAD=integraly.h random.h
+HEAD=integraly.h random.h funkce.h
 SRC=main.c integraly.c random.c
 OBJ=main.o integraly.o random.o
 

ukol-1/funkce.h

@@ -0,0 +1,23 @@
+/* jednorozmerna nezaporna funkce: R -> R+ */
+struct funkce_1d {
+	double (*f)(double);  /* fce vracejici hodnotu */
+	double a, b;  /* interval, na kterem chceme pocitat funkci (integral) */
+	double If;  /* skutecna (presna) hodnota integralu z f predpocitana analyticky */
+
+	double c;  /* funkce je na [a; b] omezena touto konstantou */
+
+	double (*g)(double);  /* na [a; b] plati: f(x) >= g(x) */
+	double Ig;  /* integral z funkce g na intervalu [a; b] */
+
+	char *jmeno;  /* symbolicky zapis funkce f */
+	char *jmeno_g;  /* symbolicky zapis funkce g */
+};
+
+/* n-rozmerna nezaporna funkce: R^dim -> R+ */
+struct funkce_nd {
+	unsigned int dim;  /* dimenze definicniho oboru */
+	double (*f)(double *);  /* fce vracejici hodnotu, rozmer pole musi byt presne dim */
+	double *a, *b;  /* vicerozmerny intergal, na kterem pocitame (integral) */
+	double c;  /* konstanta omezujici f na [a; b] */
+	char *jmeno;  /* pro hezci vypisy */
+};

ukol-1/integraly.h

@@ -1,26 +1,4 @@
-/* jednorozmerna nezaporna funkce: R -> R+ */
-struct funkce_1d {
-	double (*f)(double);  /* fce vracejici hodnotu */
-	double a, b;  /* interval, na kterem chceme pocitat funkci (integral) */
-	double If;  /* skutecna (presna) hodnota integralu z f predpocitana analyticky */
-
-	double c;  /* funkce je na [a; b] omezena touto konstantou */
-
-	double (*g)(double);  /* na [a; b] plati: f(x) >= g(x) */
-	double Ig;  /* integral z funkce g na intervalu [a; b] */
-
-	char *jmeno;  /* symbolicky zapis funkce f */
-	char *jmeno_g;  /* symbolicky zapis funkce g */
-};
-
-/* n-rozmerna nezaporna funkce: R^dim -> R+ */
-struct funkce_nd {
-	unsigned int dim;  /* dimenze definicniho oboru */
-	double (*f)(double *);  /* fce vracejici hodnotu, rozmer pole musi byt presne dim */
-	double *a, *b;  /* vicerozmerny intergal, na kterem pocitame (integral) */
-	double c;  /* konstanta omezujici f na [a; b] */
-	char *jmeno;  /* pro hezci vypisy */
-};
+#include "funkce.h"
 
 struct vysledek {
 	double I;  /* hodnota integralu */

ukol-1/main.c

@@ -21,7 +21,7 @@ double g_1(double x)
 struct funkce_1d func_1 = {
 	f_1,  /* fce vracejici hodnotu */
 	0.0, 4.0,  /* interval, na kterem chceme pocitat funkci (integral) */
-	25.0/3.0,  /* skutecna (presna) hodnota integralu z f predpocitana analyticky */
+	28.0/3.0,  /* skutecna (presna) hodnota integralu z f predpocitana analyticky */
 	3.0,  /* funkce je na [a; b] omezena touto konstantou */
 	g_1,  /* na [a; b] plati: f(x) >= g(x) */
 	8.0,  /* integral z funkce g na intervalu [a; b] */
@@ -43,9 +43,9 @@ struct funkce_1d func_2 = {
 	M_E,  /* skutecna (presna) hodnota integralu z f predpocitana analyticky */
 	5.0,  /* funkce je na [a; b] omezena touto konstantou */
 	g_2,  /* na [a; b] plati: f(x) >= g(x) */
-	0.0,  /* integral z funkce g na intervalu [a; b] */
+	21.0/8.0,  /* integral z funkce g na intervalu [a; b] */
 	"e^x + 2x",  /* symbolicky zapis funkce f */
-	"???",  /* symbolicky zapis funkce g */
+	"3.65x + 0.8",  /* symbolicky zapis funkce g */
 };
 
 
@@ -57,21 +57,27 @@ int main()
 	unsigned int pocty_kroku[POCTU_KROKU] = {10, 100, 10000, 1000000, 10000000}; /* tyto pocty kroku budeme zkouset */
 	unsigned int i, j, N;
 	char *nazev_funkce, *nazev_metody;
+	double delta;
 
 	srand(42);  /* nastavime nahodny generator na stejny seed pro reproducibilitu vysledku */
-	fprintf(stdout, "funkce\t\"pocet kroku\"\tmetoda\tintegral\todchylka\n");
 
 	/* vsechny funkce, pocty kroku a metody pro jednorozmerne funkce */
 	for(i = 0; i < FUNKCI_1D; i++) {
 		func = testovane_1d_funkce[i];
 		nazev_funkce = func->jmeno;
+		fprintf(stderr, "Jdeme integrovat funkci f(x) = %s na intervalu [%f; %f]\n", nazev_funkce, func->a, func->b);
+		fprintf(stderr, "Na intervalu by melo platit: a) f(x) <= %f  b) g(x) = %s <= f(x)  c) I(g(x)) = %f\n", func->c, func->jmeno_g, func->Ig);
+		fprintf(stderr, "Analyticky spocitana hodnota integralu = %f\n", func->If);
 
+		fprintf(stdout, "\"funkce\"\t\"pocet kroku\"\t\"metoda\"\t\"integral\"\t\"odhad odchylky\"\t\"skutecna odchylka\"\n");
 		for(j = 0; j < POCTU_KROKU; j++) {
 			N = pocty_kroku[j];
 			nazev_metody = "vzorkovani plochy";
 			result = integral_plocha(func, N);
-			fprintf(stdout, "\"%s\"\t%i\t\"%s\"\t%f\t%f\n", nazev_funkce, N, nazev_metody, result.I, result.s);
+			delta = fabs(result.I - func->If);
+			fprintf(stdout, "\"%s\"\t%i\t\"%s\"\t%f\t%f\t%f\n", nazev_funkce, N, nazev_metody, result.I, result.s, delta);
 		}
+		fprintf(stderr, "\n");
 	}
 
 	return 0;

ukol-1/random.h

@@ -1,5 +1,5 @@
 /* Jednorozmerna nahodna hodnota mezi a a b */
 double rovnom_1d(double a, double b);
 
-/* 2-rozmerna nahodna hodnota, hodnotu zapise na adresu out */
+/* n-rozmerna nahodna hodnota, hodnotu zapise na adresu out */
 void rvonom_nd(double *a, double *b, unsigned int dim, double *out);