Задачка на Си.

пятница, 15 апреля 2011

23:43

все записи пользователя в сообществе infinity235

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

double function(double c, int n);
double differentiation(double c, int n);

int main(void)
{
	int n; // Степень корня
	double c; // Число под корнем
	int i; // Счётчик попыток
	double state[100];
	
	printf("Введите число c>");
	scanf("%lf", &c);
	
	printf("Введите n - степень корня>");
	scanf("%d", &n);
	
	// Исходное предположение с/2
	state[0] = function(c/2, n);
	printf("Root: %f ", state[0]);
	printf("Diff: %f\n\n\n", differentiation(c/2, n));

	for(i=1; i<10; ++i)
	{
	state[i]=state[i-1]-function(state[i-1],n)/differentiation(state[i-1], n);
	printf("Root: %f ", state[i]);
	printf("Diff: %f\n", differentiation(state[i], n));
	}

	getch();
	return(0);
}


double function(double c, int n)
{
double fun;
fun = pow(c,(double)1/n);
return(fun);
}


double differentiation(double c, int n)
{
double dif;
dif = ((double)1/(double)n)*pow(c,((double)1/(double)n-1));
return(dif);
}

Неправильно работает цикл.
Либо где-то ошибка в условии, либо я неправильно понимаю задачу.

@темы: C++

URL

"Тридцать спиц сходятся в ступице, но как раз промеж... Я, наверное, схожу на эту постановку... http://bulvar.2... Тишина и покой. Цветущая магнолия - тяжелые крупные цветы...

Неужели не очевидно, что и то, и другое - всего лишь маск... Люмпен на рынок не смотрит, люмпен или наемный рабочий см... Я - создатель края. Не веришь? Нет, я не просто создатель...

Комментарии

16.04.2011 в 01:57

Joeman

А разве f(x) не равно x^n - c?

URL

16.04.2011 в 02:27

Joeman

Вот моё предложение решения:

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

double function(double x, double c, int n);
double differentiation(double x, double c, int n);

int main(void)
{
	int n; // Степень корня
	double c; // Число под корнем
	int i; // Счётчик попыток
	double state[100];
	setlocale(LC_ALL,"Russian");
	printf("Введите число c>");
	scanf("%lf", &c);
	double state0, state1;
	printf("Введите n - степень корня>");
	scanf("%d", &n);
	
	// Исходное предположение с/2
	state0 = c/2;
	state1 = 100;
	printf("Root: %f ", state[0]);
	printf("Diff: %f\n\n\n", differentiation(state0,c/2, n));

	//for(i=1; i<10; ++i)
	while (fabs(state1- state0) > 1.0e-6)
	{
	state0 = state1;
	state1=state0-function(state0,c,n)/differentiation(state0,c, n);
	//printf("Root: %f ", state1);
	//printf("Diff: %f\n", differentiation(state1,c, n));
	
	}
	printf("Корень степени %d из числа %f равен %f", n, c, state1);
	getch();
	return(0);
}


double function(double x, double c, int n)
{
double fun;
//fun = pow(c,(double)1/n);
fun = pow(x,n) - c;
return(fun);
}


double differentiation(double x, double c, int n)
{
double dif;
dif = n*pow(x,n-1);//((double)1/(double)n)*pow(c,((double)1/(double)n-1));
return(dif);
}

URL