помогите разобраться(Си)

[Nikolay625]

Вот есть программа для нахождения комплексных корней алгебраического уравнения методом Ньютона-Рафсона. Как тут реализовано нахождение производной?

#include <stdio.h>
#include <math.h>
#include <complex.h>
typedef float complex	Float;
void Newton_Raphson(int n, Float q[n + 1], Float p[n], float eps);
Float hor(int n, Float a[n + 1], Float x, Float* b)
{
	int i;
	Float val = a[n];
	if (b) {
		b[n - 1] = val;
	}

	for (i = n - 1; i >= 0; --i) {
		val *= x;
		val += a[i];
		if (b && i) {
			b[i - 1] = val;
		}
	}
	return val;
}
Float dhor(int n, Float a[n + 1], Float x)
{
	int i;
	Float val = n * a[n];
	for (i = n - 1; i > 0; --i) {
		val *= x;
		val += i * a[i];
	}
	return val;
}
void newton(int n, Float a[n + 1], Float* px, Float x0, float eps)
{
	Float x = x0, dx, f_val, df_val;
	int niter = 1;
	do {
		f_val = hor(n, a, x, NULL);
		df_val = dhor(n, a, x);
		if (df_val != 0) {
			dx = f_val / df_val;
		} else {
			break;
		}
		x -= dx;
	} while (++niter <= 100 && cabsf(dx) > eps);
	*px = x;
}
void Newton_Raphson(int n, Float q[n + 1], Float p[n], float eps)
{
	Float x, x0;
	int i;
	float err;

	while (n) {
		x0 = 1.0f;
		newton(n, q, &x, x0, eps);
		err = cabsf(hor(n, q, x, NULL));
		if (err > eps) {
			newton(n, q, &x, x0 * I, eps);
			err = cabsf(hor(n, q, x, NULL));
		}
		Float b[n];
		hor(n, q, x, b);
		p[n - 1] = x;
		for (i = 0; i < n; ++i) {
			q[i] = b[i];
		}
		--n;
	}
}
int main(int argc, char* argv[])
{
	int n;
	int i;
	printf("\nstepen uravnenia n=");
	scanf("%d", &n);
	Float q[n + 1];
	Float p[n];
	printf("\nvektor koefficientov:\n");
	for (i = 0; i <= n; i++) {
		float x;
		scanf("%f", &x);
		q[i] = x;
	}
	Newton_Raphson(n, q, p, 1e-6);
	printf("\nKorni uravnenija");
	printf("\nDeystv\t        Mnim\n");
	for (i = 0; i < n; i++) {
		printf("%g\t%10g\n", crealf(p[i]), cimagf(p[i]));
	}
	printf("\n");
	getch ();
	return 0;
}

[Stilet]

Почему-то у меня складывается вчепятление что нет тут никакой функции, и производная тут не вычисляется. Просто вводятся уже посчитанные на калькуляторе иксы и все.

[Nikolay625]

а за функции hor и dhor?

[Stilet]

Х.з... Метод хорд что ли...
Это надо у автора сего творения спрашивать )

[Nikolay625]

а метод Хорнера вам ни о чем не говорит? пытался гуглить, но чет никак

[Stilet]

Неа. Ни о чем )
Я не математик.

[Nikolay625]

5 автоматом гайс)))