Найти корень уравнения методом половинного деления - C (СИ)

Узнай цену своей работы

Формулировка задачи:

В программе методом половинного деления необходимо найти корень уравнения. Загвоздка в том что при вводе буквы программа автоматически приравнивает все переменные друг к другу и в итоге перебрасывает действие к последнему Else. Так же сделал защиту от дурака на примере одной переменной j, но к сожалению если ввести например 1к защита уже не работает. Пожалуйста подскажите что исправить)
Листинг программы 
  1. #include <stdlib.h>
  2.  
  3. double func(double x,double a, double b, double c, double d)
  4. {
  5.         return   ((a*x+b)+1/(c*x+d));
  6. }
  7.  
  8. int main() {
  9. double func(double x,double a, double b, double c, double d);
  10.  
  11.     int n=0;
  12.         double a;    //  переменная
  13.         double b;    //  переменная
  14.         double c;    //  переменная
  15.         double d;    //  переменная
  16.         double j;    //  координата отрезка поиска корня
  17.         double g;    //  координата отрезка поиска корня
  18.         double u;    //  вспомогательная переменная
  19.         double eps;  //  абсолютная погрешность
  20.  
  21.         printf("j=");
  22.         if (scanf ("%d", &j) != 1) {
  23.             printf("Error.\n");
  24.  
  25.         }
  26.  
  27.         printf("g=");
  28.         scanf("%lf", &g);
  29.         printf("eps=");
  30.         scanf("%lf", &eps);
  31.         printf("a=");
  32.         scanf("%lf", &a);
  33.         printf("b=");
  34.         scanf("%lf", &b);
  35.         printf("c=");
  36.         scanf("%lf", &c);
  37.         printf("d=");
  38.         scanf("%lf", &d);
  39.  
  40. if (func(j,a,b,c,d)*func(g,a,b,c,d)<=0)
  41. {
  42.  
  43.         do
  44. {
  45.         u=(j+g)/2;
  46.         if (func(u,a,b,c,d)*func(j,a,b,c,d)<=0)
  47.             g=u;
  48.         else j=u;
  49.  
  50.                 n+=1;
  51.         printf("n=%d\n", n);
  52.         printf("u=%lf\n", u);
  53.         }
  54.         while (fabs(j-g)>=eps);}
  55.  
  56. else {
  57.         printf( "\n Neverno zadany nachalnye usloviya" );  // если ни одно из выше-перечисленных условий не подошло, то программа покажет этот вариант ответа
  58.     }
  59.                 getch();
  60.         return 0;
  61. }

Решение задачи: «Найти корень уравнения методом половинного деления»

textual
Листинг программы 
  1. #include <stdio.h>
  2. #include <stdlib.h>
  3. #include <ctype.h>
  4. #include <math.h>
  5.  
  6. double func(double x, double a, double b, double c, double d)
  7. {
  8.     return   ((a*x + b) + 1 / (c*x + d));
  9. }
  10.  
  11. int print_error() {
  12.     printf("\nIncorrect number...\n");
  13.     return 0;
  14. }
  15.  
  16. double get_number(char* s)
  17. {
  18.     char arr_int[10] = { 0 }, arr_fract[10] = { 0 };
  19.     char ch;
  20.     unsigned short fract = 0, i = 0, j = 0;
  21.     int sign = 1;
  22.     double d;
  23.  
  24.     printf("\nEnter %s: ", s);
  25.  
  26.     while ((ch = getchar()) != '\n') {
  27.  
  28.         if ((ch == '-') && (i == 0)) { sign = -1;}
  29.         else if (isdigit(ch) || (ch == '.')) {
  30.             
  31.             if ((ch == '.') && (i != 0)) fract = 1;
  32.             if ((ch == '.') && !( (i==1) && (sign == 1) ) )  fract = 1;
  33.             
  34.             if ( isdigit(ch) && fract == 1) { arr_fract[j] = ch; j++; }
  35.             if ( isdigit(ch) && fract == 0) { arr_int[i] = ch; i++; }
  36.  
  37.         }
  38.         else return print_error();
  39.  
  40.         
  41.  
  42.         
  43.  
  44.  
  45.     }
  46.     //printf("\n");
  47.     d = atol(arr_int) + atol(arr_fract)*pow(10, -1.0*j);
  48.  
  49.     return d * sign;
  50. }
  51.  
  52. int main() {
  53.  
  54.     int n = 0;
  55.     double a;
  56.     if ((a = get_number("a")) == 0) { system("pause"); return -1; }    //  переменная
  57.     double b;   
  58.     if ((b = get_number("b")) == 0) { system("pause"); return -1; }         //  переменная
  59.  
  60.     double c = get_number("c");    //  переменная
  61.  
  62.     double d = get_number("d");    //  переменная
  63.  
  64.     double j = get_number("j");    //  координата отрезка поиска корня
  65.  
  66.     double g = get_number("g");    //  координата отрезка поиска корня
  67.  
  68.     double u = get_number("u");    //  вспомогательная переменная
  69.  
  70.     double eps = get_number("eps");  //  абсолютная погрешность
  71.  
  72.     if (func(j, a, b, c, d)*func(g, a, b, c, d) > 0) {
  73.         while (fabs(j - g) >= eps) {
  74.             u = (j + g) / 2;
  75.             if (func(u, a, b, c, d)*func(j, a, b, c, d) <= 0)
  76.                 g = u;
  77.             else j = u;
  78.  
  79.             n += 1;
  80.             printf("n=%d\n", n);
  81.             printf("u=%lf\n", u);
  82.         }
  83.     }
  84.     else {
  85.         printf("\n Neverno zadany nachalnye usloviya\n");
  86.     }
  87.     
  88.                                             
  89.     system("pause");
  90.     return 0;
  91. }

Объяснение кода листинга программы

  1. Включаются необходимые заголовочные файлы
  2. Определяется функция func(double x, double a, double b, double c, double d), которая представляет собой выражение под корнем
  3. Определяется функция print_error(), которая выводит сообщение об ошибке и возвращает 0
  4. Определяется функция get_number(char* s), которая считывает число с плавающей точкой или знак минус и возвращает его в переменной типа double
  5. В функции main() считываются значения переменных a, b, c, d, j, g и u с помощью функции get_number(string)
  6. Определяются значения переменных eps и n
  7. Выполняется проверка на корректность начального условия (если функция func(j, a, b, c, d) умноженная на func(g, a, b, c, d) больше 0)
  8. Если условие выполняется, то в цикле while происходит бинарный поиск корня методом половинного деления
  9. Внутри цикла вычисляются значения u, j и g, а также выводится информация о количестве итераций и значениях переменных u, j и g
  10. Если условие в начале функции main() не выполняется, то выводится сообщение об ошибке
  11. В конце программы выполняется пауза и возвращается 0

ИИ поможет Вам:

  • решить любую задачу по программированию
  • объяснить код
  • расставить комментарии в коде
  • и т.д
Попробуйте бесплатно

Оцени полезность:

15   голосов , оценка 4.2 из 5

Нужна аналогичная работа?

Оформи быстрый заказ и узнай стоимость

Бесплатно
Оформите заказ и авторы начнут откликаться уже через 10 минут
Похожие ответы