Nullstellen-Suche

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114

#include <stdio.h>
#include <math.h>

int fn;

double f(double x)
{
	if(fn<=5 && fn>=1)
	{
		switch(fn)
		{
			case 1:{
						return sin(x);
						break;
				   }
			case 2:{
						return exp(x)-5*x+1;
						break;
				   }
			case 3:{
						return x*x-2;
						break;
				   }
			case 4:{
						return (double)1/7+1/x;
						break;
				   }
			case 5:{
						return 2*x-2;
				   }
		}
	}
	return 0;
}

struct rueckgabe
{
	int		anzit;
	int		ok;
	double	x0;
};

struct rueckgabe nullstelle(double a,double b,double g)
{
	struct rueckgabe merker;
	double c=0;

	merker.anzit=0;
	merker.ok=1;

	do{

		c=(a+b)/2;

		if (f(a)>0 && f(b)<0)
		{
			if(f(c)>0)
			{
				a=c;
			}
			else
			{
				b=c;
			}
		}else if(f(a)<0 && f(b)>0)
		{
			if(f(c)<0)
			{
				a=c;
			}
			else
			{
				b=c;
			}
		}else if(f(a)==0)
		{
			c=b=a;  // x0=a
		}
		else if(f(b)==0)
		{
			c=a=b;  // x0=b
		}
		else merker.ok=0;

		merker.anzit++;

	}while (g<fabs(f(a)-f(b)) && merker.ok!=0);
	
	merker.x0=c;//oder mit linearer Interpolation: a-f(a)*(b-a)/(f(b)-f(a));
	return merker;
}

void main(void)
{
	double a,b,g;
	struct rueckgabe merker;
	
	printf("Waehlen Sie eine Funktion (1-5): ");
	scanf("%d",&fn);

	printf("\nUntere Intervallgrenze    a = ");
	scanf("%lf",&a);
	printf("Obere Intervallgrenze     b = ");
	scanf("%lf",&b);

	printf("Die absolute Genauigkeit  g = ");
	scanf("%lf",&g);
	
	merker=nullstelle(a,b,g);
	if(merker.ok!=0)
	printf("\nErgebnis                  E = %.10lf\nAnzahl der Iterationen    I = %d\nFunktionswert fuer x0     F = %.10lf\n\n",merker.x0,merker.anzit,f(merker.x0));
	else
	printf("\nKein Vorzeichenwechsel an den Intervallgrenzen!\n\n");
}