VPF::[C\C++] Задача решающая трансцендентные ур-я

Kot

Дата 10.10.2006, 01:09 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 3
Регистрация: 29.7.2006

Репутация: нет
Всего: нет

По алгоритмам вопросов нет - они простые. Проблемы возникают при использовании типа double: функции sqrt и acos не работают. Как быть? smile

Код


#include <conio.h>
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
 
//uravnenie
// acos(x)-sqrt(1-(0.3*(x*x*x)))=0
 
//pervaja proizvodnaya funkcii
//((0.9*(x*x))/(2*(sqrt(1-(0.3*(x*x*x))))))-(1/(sqrt(1-(x*x))));
 
//Metod Dihotomii
void findRootDihotomy()
{
int N;//Chislo iterazii
double x;
double f_a;
double A;
double B;
double E;
clrscr();
printf("Vvedite levuiu granicu:\n");
scanf("%f",&A);
printf("Vvedite pravuiu granicu:\n");
scanf("%f",&B);
printf("Vvedite pogreshnost':\n");
scanf("%f",&E);
f_a=acos(A)-sqrt(1-(0.3*(A*A*A)));
N=0;
do
{
x=(A+B)*.5;
N++;
if ((f_a-( acos(x)-sqrt(1-(0.3*(x*x*x))) ))<.0000000)
{
B=x;
}
else
{
A=x;
f_a=( acos(x)-sqrt(1-(0.3*(x*x*x))) );
}
}
while (abs(B-A)>E);
printf("Metod dihotomii:\n");
printf("x=%f\n",x);
printf("Chislo iterazii=%d\n",N);
}
 
//Metod N'iutona
 
void findRootNewton()
{
double A;
double B;
double E;
double x_p;//predidushee znachenie argumenta
double x_s;//tekushee znachenie argumenta
double f;
int N;//Chislo iterazii
int flag_f;
int count;
clrscr();
printf("Vvedite levuiu granicu:\n");
scanf("%f",&A);
printf("Vvedite pravuiu granicu:\n");
scanf("%f",&B);
printf("Vvedite pogreshnost':\n");
scanf("%f",&E);
N=0;
flag_f=0;

printf("Nachal'noe znachenie kornya:\n");
printf("1 - Znachenie levoi granicy\n");
printf("2 - Znachenie pravoi granicy\n");
scanf("%d",&count);
if (count==1)
x_p=A;
if (count==2)
x_p=B;

do
{
if(flag_f==1)
x_p=x_s;
f=((0.9*(x_p*x_p))/(2*(sqrt(1-(0.3*(x_p*x_p*x_p))))))-(1/(sqrt(1-(x_p*x_p))));
x_s=( x_p-( ( acos(x_p)-sqrt(1-(0.3*(x_p*x_p*x_p))) )/f ) );
N++;
}
while (E<abs(x_p-x_s));
printf("Metod N'iutona:\n");
printf("x=%f\n",x_s);
printf("Chislo iterazii=%d\n",N);
}
 
//Metod Iterazii
void findRootIter()
{
double A;
double B;
double E;
double x_p;//predidushee znachenie argumenta
double x_s;//tekushee znachenie argumenta
int N;//Chislo iterazii
int flag_f;
int count;
clrscr();
printf("Vvedite levuiu granicu:\n");
scanf("%f",&A);
printf("Vvedite pravuiu granicu:\n");
scanf("%f",&B);
printf("Vvedite pogreshnost':\n");
scanf("%f",&E);
N=0;
flag_f=0;

printf("Nachal'noe znachenie kornya:\n");
printf("1 - Znachenie levoi granicy\n");
printf("2 - Znachenie pravoi granicy\n");
scanf("%d",&count);
if (count==1)
x_p=A;
if (count==2)
x_p=B;

do
{
if(flag_f==1)
x_p=x_s;
x_s=( acos(x_p)-sqrt(1-(0.3*(x_p*x_p*x_p))) );
N++;
}
while (E<abs(x_p-x_s));
printf("Metod iterazii:\n");
printf("x=%f\n",x_s);
printf("Chislo iterazii=%d\n",N);
}
 
 
void main ()
{
int pos;
pos=0;
clrscr();
printf(" ---------------------------------------------------- \n");
printf("|                ...:::MENU:::...                    |\n");
printf(" ---------------------------------------------------- \n");
 
printf("Sdelaite vybor:\n");
 
//print punkty menu
printf("1 - Metod Dihotomii\n");
printf("2 - Metod N'utona\n");
printf("3 - Metod prostikh iterazii\n");
printf("4 - Vihod iz programmy\n");
printf(" ---------------------------------------------------- \n");
 
 
do
{
printf("Vybor:\n");
scanf ("%d",&pos);
switch(pos)
{
case 1: findRootDihotomy(); break;
case 2: findRootNewton(); break;
case 3: findRootIter(); break;
case 4: break;
}
}
while (pos!=1 && pos!=2 && pos!=3 && pos!=4);
}

Earnest

Дата 10.10.2006, 08:00 (ссылка)

(нет голосов)

Загрузка ...

Эксперт

Профиль
Группа: Экс. модератор
Сообщений: 5962
Регистрация: 17.6.2005
Где: Рязань

Репутация: 2
Всего: 183

Для домашних заданий, курсовых, существует "Центр Помощи"

Тема перенесена!

--------------------

...

ptr

Дата 10.10.2006, 13:06 (ссылка)

(нет голосов)

Загрузка ...

Опытный

Профиль
Группа: Участник
Сообщений: 661
Регистрация: 31.5.2006
Где: Новосибирск

Репутация: нет
Всего: 14

Цитата(Kot @ 10.10.2006, 05:09

)

Проблемы возникают при использовании типа double: функции sqrt и acos не работают.

В каком смысле не работают?

abs() надо поменять на fabs().

--------------------

Единственный способ определить границы возможного - это выйти за эти границы, в невозможное.
Артур Кларк.

0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| Центр помощи \| Следующая тема »