newton-verfahren

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

float a[6];
float b[5];
float sw;
float lsg;

void eingabe(){
	for(int i = 0; i<6; i++)
	{
		printf("%i. Koeffizienten eingeben: ", i + 1);
		scanf("%f",&a[i]);
	}
}

void ableitung(){
	for(int i = 0; i <5; i++)
	{
		b[i] = a[i] * (5-i);
	}
}

float newton(float x){
	float aFktWert = 1;
	float bFktWert = 1;

	while(abs(aFktWert) > 0.01)
	{
		aFktWert = 0;
		bFktWert = 0;

		for(int i = 0; i<6; i++)
		{
			aFktWert += a[i]*pow(x,5-i);
		}

		for(int i = 0; i<5; i++)
		{
			bFktWert += b[i]*pow(x,4-i);
		}

		x = x - (aFktWert/bFktWert);
	}

	return x;
}

int main(void){
	eingabe();
	ableitung();

	printf("\nFunktion: f(x) = %fx^5+%fx^4+%fx^3+%fx^2+%fx+%f\n",a[0],a[1],a[2],a[3],a[4],a[5]);
	printf("Ableitung: f'(x) = %fx^4+%fx^3+%fx^2+%fx+%f\n",b[0],b[1],b[2],b[3],b[4]);

	printf("\nStartwert eingeben: ");
	scanf("%f",&sw);

	lsg = newton(sw);
	printf("\nNulstelle bei x = %.2f",lsg);
	return 0;
}