y=x^3-4x+1という関数の計算方法

y=x^3-4x+1という関数において、a<=x<=bの範囲でn分割し、yの最小値を求めるプログラムを作成しています。for文を使って、計算値の比較を行いたいのですが、関数の計算方法が分かりません。


int main(int argc,char* argv[]){
   float a = atof(argv[1]);
   float b = atof(argv[2]);
   int n = atoi(argv[3]);
   float min = 0;

   for(int i=a;i<=b;i+=n){
    
   }
}

変数x,yを定義してyに上記の関数を代入し、xにiを代入してfor文を回そうとしましたがうまくいきませんでした。

maisumakun

2023/04/28 06:46

> 関数の計算方法が分かりません。どこまで書いてみて、何に詰まっている状況でしょうか？

can110

2023/04/28 06:49

その関数の「何が」分からないのでしょうか？ - 関数そのもの（x=-1のときのyの値の計算方法）が分からない。 - その関数をC言語での演算子を使った計算式として表現する方法が分からない。また、y=xという関数だった場合、プログラムは作成できますか？

opopopop

2023/04/28 06:49

書いて見たのが上記のところまでです。y=x^3-4x+1という関数の定義方法がわかりません。

1T2R3M4

2023/04/28 06:52 編集

>うまくいきませんでした。うまくいかなかったコードを提示していただけませんか。また、xが2の時のyは計算できますか。

opopopop

2023/04/28 06:50

c言語において関数をどのように表現すれば良いのかがわからないです。

maisumakun

2023/04/28 06:50

> 変数x,yを定義してyに上記の関数を代入し実際に書いたコードと、どのような事態になったのかをご提示ください。

opopopop

2023/04/28 07:19

#include<stdio.h> #include<math.h> #include<stdlib.h> int main(int argc,char* argv[]){ float a = atof(argv[1]); float b = atof(argv[2]); int n = atoi(argv[3]); float min = 0; float x; for(int i=a;i<=b;i+=n){ float y = (i*i*i)-(4*i)+1; if(min>y){ min = y; x = i; } else{ min = min; x = x; } } printf("[x,y]=[%.4f,%.4f]\n",x,min); } このような形で書きました。a,bにたとえば2.3,8.9と入力した時にxに0が出力されてしまいます。a,bはxの範囲ですので0が入るのはおかしいのですがどこがおかしいのかわかりません。ちなみにa,bに1.0,2.0を入れた時はうまくいきました。

行動規範の内容に同意します

回答3件

たとえばa=0,b=2,n=3の場合、すなわち0から2までを3分割する場合、xは0, 0.67, 1.33, 2の4つをとるとのことですが

C++
1for(int i=a;i<=b;i+=n){

そのループの中でiがどのような値をとるのかprintなりで確認してみてください。
そして、なぜi+=nとしたのか？を考えてみてください。
（型の問題はひとまずおいときます）

投稿2023/04/28 07:37

編集2023/04/28 08:03

can110

総合スコア38262

opopopop

2023/04/28 07:50

範囲を3分割ですので、0から2まで3分割すると、0,0.67,1.33,1.99といったようになると思います。 i+=nの部分は修正いたしました。

行動規範の内容に同意します

ベストアンサー

a,bはxの範囲ですので0が入るのはおかしいのですがどこがおかしいのかわかりません。

xを初期化していませんので、ループの最初でelseに進んだ場合、初期化前のxに対してx = xが行われてしまいますが、これは未定義の操作です（何が起きても文句を言えません）。

投稿2023/04/28 07:22

maisumakun

総合スコア145183

opopopop

2023/04/28 07:48

ご指摘いただいた内容と他の方の指摘を参照し、コードを書き直したのですが、実行しても何も表示されなくなってしまいました。 #include<stdio.h> #include<math.h> #include<stdlib.h> int main(int argc,char* argv[]){ if(argc != 4){ fprintf(stderr,"引数の数が違います\n"); exit(1); } float a = atof(argv[1]); float b = atof(argv[2]); int n = atoi(argv[3]); float min; float x = a; float delta = (b-a)/n; for(float i=a;i<=b;i+delta){ float y = (i*i*i)-(4*i)+1; if(min>y){ min = y; x = i; } else{ min = min; x = x; } } printf("[x,y]=[%.4f,%.4f]\n",x,min); return 0; }

maisumakun

2023/04/28 07:50

i+deltaでiの値は変わりませんので、同じiでループが続いて脱出不能となっている状況です。

opopopop

2023/04/28 07:53

i+=deltaに修正したらうまくいきました！ありがとうございます。

行動規範の内容に同意します

C++
1#include <iostream>
2
3float fun(double x) {
4  return x*x*x - 4*x + 1;
5}
6
7int main(int argc,char* argv[]){
8   float a = atof(argv[1]);
9   float b = atof(argv[2]);
10   int n = atoi(argv[3]);
11   bool at_first = true;
12
13   float min;
14   float delta = (b - a)/n;
15   for(int i = 0; i < n; ++i){
16     float x = a + delta*i;
17     float y = fun(x);
18     if ( at_first ) { min = y; at_first = false; } // 最初の値なら問答無用にminに代入
19     else if ( y < min ) { min = y; } // さもなくば 記録更新
20   }
21   std::cout << min << std::endl;
22}