Метод Ньютона

@Pain'ka · Регистрация: 23.11.2011

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

Привет всем! Найти 3 корня методом Ньютона: f(x)=tan(2.2082*x) - 3.2258*x.
Вроде программу написала по алгоритму, но что то не так. Подскажите, что, пожалуйста!!!!

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
#include <math.h>
#include <iostream>
#include <stdio.h>
#define q 2.2082
#define w 3.2258
#define n  10
 using namespace std;
 
double f(double x) {
 
        return  sin(x*q)/cos(q*x)-w*x;
}
 
double f1(double x) {
 
        return   q/(cos(x)*cos(x))-w;
}
 
double f2(double x) {
 
        return  (2*q*sin(x))/(cos(x)*cos(x)*cos(x));
}
 
int main() {
        double a,b,c,h,x1,x2,eps = 0.0000001;
        cout<<"a="; cin>>a;// левый конец отрезка
        cout<<"b="; cin>>b;// правый конец отрезка
       
        h=(b-a)/n;
        for(int i=1;i<=10;i++){
            x1=a+h*(i-1);
            x2=a+h*i;
            
        
        if(f(x1)*f2(x1)>0) c=x1;
        else c=x2;
        do {
                c=c-f(c)/f1(c);
             
        }
        while (fabs(f(c))>=eps);
                //cout<<"c="<<c<<"\n";
                printf("%e \n",c);
                }
               
        return 0;
}

@9-loh · 25.02.2012, 21:14

в интервале [a,b] точно есть корень?

нужна только одна производная(первая)

@Том Ардер · 25.02.2012, 21:17

Ошибки в производных:

Сообщение от Pain'ka

double f1(double x)
{
return q/(cos(x)*cos(x))-w;
}
double f2(double x)
{
return (2*q*sin(x))/(cos(x)*cos(x)*cos(x));
}

Аргумент у sin, cos д.б. qx

@Pain'ka · 25.02.2012, 21:25 **[ТС]**

Сообщение от 9-loh

в интервале [a,b] точно есть корень?
нужна только одна производная(первая)

условие выбора идет по второй производной,
интервал надо как то подобрать, а как я не знаю, примерный график накидала и по нему выбирала

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

Сообщение от Том Ардер

Аргумент у sin, cos д.б. qx

точно! спасибо!

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

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL3_image 8Observer8 10.02.2026 Содержание блога Библиотека SDL3_image содержит инструменты для расширенной работы с изображениями. Пошагово создадим проект для загрузки изображения формата PNG с альфа-каналом (с прозрачным. . .	Установка Qt-версии Lazarus IDE в Debian Trixie Xfce volvo 10.02.2026 В общем, достали меня глюки IDE Лазаруса, собранной с использованием набора виджетов Gtk2 (конкретно: если набирать текст в редакторе и вызвать подсказку через Ctrl+Space, то после закрытия окошка. . .	SDL3 для Web (WebAssembly): Работа со звуком через SDL3_mixer 8Observer8 08.02.2026 Содержание блога Пошагово создадим проект для загрузки звукового файла и воспроизведения звука с помощью библиотеки SDL3_mixer. Звук будет воспроизводиться по клику мышки по холсту на Desktop и по. . .	SDL3 для Web (WebAssembly): Основы отладки веб-приложений на SDL3 по USB и Wi-Fi, запущенных в браузере мобильных устройств 8Observer8 07.02.2026 Содержание блога Браузер Chrome имеет средства для отладки мобильных веб-приложений по USB. В этой пошаговой инструкции ограничимся работой с консолью. Вывод в консоль - это часть процесса. . .
SDL3 для Web (WebAssembly): Обработчик клика мыши в браузере ПК и касания экрана в браузере на мобильном устройстве 8Observer8 02.02.2026 Содержание блога Для начала пошагово создадим рабочий пример для подготовки к экспериментам в браузере ПК и в браузере мобильного устройства. Потом напишем обработчик клика мыши и обработчик. . .	Философия технологии iceja 01.02.2026 На мой взгляд у человека в технических проектах остается роль генерального директора. Все остальное нейронки делают уже лучше человека. Они не могут нести предпринимательские риски, не могут. . .	SDL3 для Web (WebAssembly): Вывод текста со шрифтом TTF с помощью SDL3_ttf 8Observer8 01.02.2026 Содержание блога В этой пошаговой инструкции создадим с нуля веб-приложение, которое выводит текст в окне браузера. Запустим на Android на локальном сервере. Загрузим Release на бесплатный. . .	SDL3 для Web (WebAssembly): Сборка C/C++ проекта из консоли 8Observer8 30.01.2026 Содержание блога Если вы откроете примеры для начинающих на официальном репозитории SDL3 в папке: examples, то вы увидите, что все примеры используют следующие четыре обязательные функции, а. . .

@9-loh 2 / 2 / 0 Регистрация: 18.12.2010 Сообщений: 31
	25.02.2012, 21:14
	в интервале [a,b] точно есть корень? нужна только одна производная(первая) 0