Метод половинного деления

@Karichan · Регистрация: 06.10.2015

Студворк — интернет-сервис помощи студентам

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

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
//---------------------------------------------------------------------------
 
#pragma hdrstop
#include "stdio.h"
#include "stdlib.h"
#include "math."
#include <conio.h>
 
//---------------------------------------------------------------------------
 
#pragma argsused
double uravnenie (double );  //прототип функции
double met_polovinnogo_deleniya(double a, double b, double eps);
 
int main(int argc, char* argv[])
{
    double x_min=-2.0, x_max=1.0, h = 0.05;
    int steps =   (x_max - x_min)/h + 1;
    double * arg;
    arg = (double *)malloc(steps*sizeof(double)); 
    if(arg==NULL)
    {
        printf("Error in memory allocation");
        return -1;
    }
    double * func;
    func = (double *)malloc(steps*sizeof(double)); 
    if(func==NULL)
    {
        printf("Error in memory allocation");
        return -1;
    }
    int i=0;
    for(double x=x_min;x<=x_max; x+=h, i++)
    {
         arg[i] = x;
         func[i] = uravnenie(x);
    }
    printf ("X\tY\n", arg[i], func[i]);
 
    for(int i=0; i<steps; i++)
    {
        printf ("%lf\t%lf\n", arg[i], func[i]);
    }
     double eps = 0.001;
    double arg1 =  met_polovinnogo_deleniya(-0.95, -0.9,eps);
    printf(" arg1 = %lf \n", arg1);
    system("pause");
    return 0;
}
//---------------------------------------------------------------------------
 
dou double uravnenie (double a)
{
    double res;
    res = x*x-sin(5*x);
    return res;
}
double met_polovinnogo_deleniya(double a, double b, double eps)
{
    double res=0;
      double z;
                while (fabs(a-b)>=eps);
    {
        z=(a+b)/2;
        if (f(z)*f(a)<=0) b=z;
        else a=z;
 
        printf(" x = %lf \n", x);
        }
        system("pause");
    return res;
}

@zss · 19.10.2015, 13:48

Сообщение от Karichan

while (fabs(a-b)>=eps);

Точка с запятой в конце означает ПУСТОЕ тело цикла.
Пользуйтесь отладчикоом для нахождения таких ошибок.
https://www.cyberforum.ru/cpp-... 62479.html

@Karichan · 19.10.2015, 14:17 **[ТС]**

Visual Studio, к сожалению, нет возможности пользоваться.

Добавлено через 22 минуты
Как-то странно элементы выделились.
Немного поправила. Была путаница в обозначениях.

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
//---------------------------------------------------------------------------
 
#pragma hdrstop
#include "stdio.h"
#include "stdlib.h"
#include "math.h"
#include <conio.h>
 
//---------------------------------------------------------------------------
 
#pragma argsused
double uravnenie (double );  //прототип функции
double met_polovinnogo_deleniya(double a, double b, double eps);
 
int main(int argc, char* argv[])
{
    double x_min=-2.0, x_max=1.0, h = 0.05;
    int steps =   (x_max - x_min)/h + 1;
    double * arg;
    arg = (double *)malloc(steps*sizeof(double)); //выделить память под массив аргументов функции
    if(arg==NULL)
    {
        printf("Error in memory allocation");
        return -1;
    }
    double * func;
    func = (double *)malloc(steps*sizeof(double)); //выделить память под массив аргументов функции
    if(func==NULL)
    {
        printf("Error in memory allocation");
        return -1;
    }
    int i=0;
    for(double x=x_min;x<=x_max; x+=h, i++)
    {
         arg[i] = x;
         func[i] = uravnenie(x);
    }
    printf ("X\tY\n", arg[i], func[i]);
 
    for(int i=0; i<steps; i++)
    {
        printf ("%lf\t%lf\n", arg[i], func[i]);
    }
     double eps = 0.001;
    double arg1 =  met_polovinnogo_deleniya(-0.95, -0.9,eps);
    printf(" arg1 = %lf \n", arg1);
    system("pause");
    return 0;
}
//---------------------------------------------------------------------------
 
double uravnenie (double a)
{
    double res;
    res = a*a-sin(5*a);
    return res;
}
double met_polovinnogo_deleniya(double a, double b, double eps)
{
    double res=0;
      double z;
                while (fabs(a-b)>=eps)
    {
        z=(a+b)/2;
        if (res(z)*res(a)<=0) b=z;
        else a=z;
 
        printf(" a = %lf \n", a);
        }
        system("pause");
    return res;
}

Добавлено через 3 минуты

C++
1
 if (res(z)*res(a)<=0) b=z;

Исправила на

C++
1
if (uravnenie(z)*uravnenie(a)<=0) b=z;

Вроде бы, теперь всё работает более-менее.

Новые блоги и статьи Все статьи Все блоги /
Семь CDC на одном интерфейсе: 5 U[S]ARTов, 1 CAN и 1 SSI Eddy_Em 18.02.2026 Постепенно допиливаю свою "многоинтерфейсную плату". Выглядит вот так: https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11617&stc=1&d=1771445347 Основана на STM32F303RBT6. На борту пять. . .	Символьное дифференцирование igorrr37 13.02.2026 / * Программа принимает математическое выражение в виде строки и выдаёт его производную в виде строки и вычисляет значение производной при заданном х Логарифм записывается как: (x-2)log(x^2+2) -. . .	Камера Toupcam IUA500KMA Eddy_Em 12.02.2026 Т. к. у всяких "хикроботов" слишком уж мелкий пиксель, для подсмотра в ESPriF они вообще плохо годятся: уже 14 величину можно рассмотреть еле-еле лишь на экспозициях под 3 секунды (а то и больше),. . .	И ясному Солнцу zbw 12.02.2026 И ясному Солнцу, и светлой Луне. В мире покоя нет и люди не могут жить в тишине. А жить им немного лет.
«Знание-Сила» zbw 12.02.2026 «Знание-Сила» «Время-Деньги» «Деньги -Пуля»	SDL3 для Web (WebAssembly): Подключение Box2D v3, физика и отрисовка коллайдеров 8Observer8 12.02.2026 Содержание блога Box2D - это библиотека для 2D физики для анимаций и игр. С её помощью можно определять были ли коллизии между конкретными объектами и вызывать обработчики событий столкновения. . . .	SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL_LoadPNG (без SDL3_image) 8Observer8 11.02.2026 Содержание блога Библиотека SDL3 содержит встроенные инструменты для базовой работы с изображениями - без использования библиотеки SDL3_image. Пошагово создадим проект для загрузки изображения. . .	SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL3_image 8Observer8 10.02.2026 Содержание блога Библиотека SDL3_image содержит инструменты для расширенной работы с изображениями. Пошагово создадим проект для загрузки изображения формата PNG с альфа-каналом (с прозрачным. . .