質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

ただいまの
回答率

89.64%

ニュートン法のプログラミング

解決済

回答 1

投稿

  • 評価
  • クリップ 1
  • VIEW 241

challenge326

score -12

ニュートン法のは一般に、初期値に近い解を高速に計算する計算方法である。
区間[a; b] 内にある解を求める為に、次のような方法を考える。

  1. 区間[a; b] 内の1025 個点xn = a + (b 􀀀 a)  n=1024, (0  n  1024) 、そのな
    かでjf(xn)j が最小であるものを返す関数maekeisan を作る。
  2. 関数maekeisan で得られた点を初期値として、ニュートン法を計算する。
    問い:/* houteisiki05.c */ 以下のプログラムを完成させよ。

include <stdio.h>

include <math.h>

define EPS 1e-8 /* 打ち切り誤差*/

define LIMIT 50 /* 打ち切り回数*/

double f(double x){ return sin(x); }
double df(double x){ return cos(x); }
double maekeisan(double (*f)(double), double a, double b){
double x,saisho,xx;
int n;
saisho=fabs(f(a)); 【A 】;
for(n=0;【B 】;n++){
x=a+(b-a)*n/1024;
if (fabs(f(x))<saisho) {xx=x; saisho=【C 】}
}
return xx;
}
double newton2(double (*f)(double), double (*g)(double), double a, double b){
double x0,x1;
int k;
x0=【D 】;
printf("初期値=%lf \n",x0);
for (k=1;k<=LIMIT;k++) {
x1=x0-f(x0)/df(x0);
if (fabs(x1-x0)<EPS) return x1;
x0=x1;
}
if (k>LIMIT) { printf("収束しない\n"); return 9999;}
}
main(){
printf("x=%lf \n",newton2(f,df,3,7));
}
実行結果
初期値=3.140625
x=3.141593

  • 気になる質問をクリップする

    クリップした質問は、後からいつでもマイページで確認できます。

    またクリップした質問に回答があった際、通知やメールを受け取ることができます。

    クリップを取り消します

  • 良い質問の評価を上げる

    以下のような質問は評価を上げましょう

    • 質問内容が明確
    • 自分も答えを知りたい
    • 質問者以外のユーザにも役立つ

    評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

    質問の評価を上げたことを取り消します

  • 評価を下げられる数の上限に達しました

    評価を下げることができません

    • 1日5回まで評価を下げられます
    • 1日に1ユーザに対して2回まで評価を下げられます

    質問の評価を下げる

    teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。

    • プログラミングに関係のない質問
    • やってほしいことだけを記載した丸投げの質問
    • 問題・課題が含まれていない質問
    • 意図的に内容が抹消された質問
    • 過去に投稿した質問と同じ内容の質問
    • 広告と受け取られるような投稿

    評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。

    質問の評価を下げたことを取り消します

    この機能は開放されていません

    評価を下げる条件を満たしてません

    評価を下げる理由を選択してください

    詳細な説明はこちら

    上記に当てはまらず、質問内容が明確になっていない質問には「情報の追加・修正依頼」機能からコメントをしてください。

    質問の評価を下げる機能の利用条件

    この機能を利用するためには、以下の事項を行う必要があります。

質問への追記・修正、ベストアンサー選択の依頼

  • m.ts10806

    2020/01/16 15:04

    宿題は先生に聞いてください。

    キャンセル

  • cateye

    2020/01/16 15:23 編集

    三つも立て続けに質問して、回答があってもプログラム組めるんですか?
    誰もソースなんか提供しませんよ・・・
    https://teratail.com/help/avoid-asking

    キャンセル

  • 退会済みユーザー

    2020/01/23 01:52

    複数のユーザーから「やってほしいことだけを記載した丸投げの質問」という意見がありました
    「質問を編集する」ボタンから編集を行い、調査したこと・試したことを記入していただくと、回答が得られやすくなります。

回答 1

checkベストアンサー

-3

回答

#include<iostream>
#include<stdio.h>
#include<math.h>
using namespace std;

# define EPS 1e-8 // truncation error
# define LIMIT 50 // truncation number

double f(double x){return sin(x);}
double df(double x){return cos(x);}

double pre(double(*f)(double), double a, double b){
    double x, mim, xx;
    int n;

    mim = fabs(f(a));
    xx = fabs(f(b));

    for(n = 0; mim < xx; n++){
        x = a + (b - a) * n / 1024;
        if(fabs(f(x)) < mim){
            xx = x;
            mim = fabs(f(x));
        }
    }
    return xx;
}

double newton2(double(*f)(double), double(*g)(double), double a, double b){
    double x0, x1;
    int k;

    x0 = pre(f, a, b);
    cout<<"Initialize value: "<<x0<<endl;

    for(k = 1; k <= LIMIT; k++){
        x1 = x0 - f(x0)/ df(x0);

        if(fabs(x1 - x0) < EPS)
            return x1;

        x0 = x1;
    }

    if(k > LIMIT){
        cout<<"Couldn't converge"<<endl;
        return 9999;
    }

}

int main()
{
    cout<<"x: "<<newton2(f, df, 3, 7)<<endl;

    return 0;
}

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

15分調べてもわからないことは、teratailで質問しよう!

  • ただいまの回答率 89.64%
  • 質問をまとめることで、思考を整理して素早く解決
  • テンプレート機能で、簡単に質問をまとめられる