for文の修正について

Question

コード
#include<studio>
using namespace std;
int main() {
int i,k,a[2];
double xl,xr,xc;
double fl,fr,fc,ss;

k=0;
for(i=0;i<10;i++){
fl=sin(i)+1/sqrt(2);
fc=sin(i+1)+1/sqrt(2);
if(fl*fc<0){
a[k]=i;
k++;
}
}
for(k=0;k<2;k++){
ss=1.0;
xl=a[k];
xr=a[k]+1;
while(ss>1.0e-10){
xc=(xr+xl)/2;
fl=sin(xl)+1/sqrt(2);
fr=sin(xr)+1/sqrt(2);
fc=sin(xc)+1/sqrt(2);
if(fl*fc<0){
xr=xc;
}
else{
xl=xc;
}
ss=fabs(1-xr/xl);
}
}
}

0<x<2πの範囲で
f(x)=sinx+1/√2=0となるxを求めるコードを書いたのですが、二つ目のfor文が機能していない状態です。(解は２つなのですが片方しか出てこないです。)
同様にg(x)=x^3-5x^2/12-x/12+1/48=0となるようなコードも書いたのですが、これに関しては1つも正しい答えが出ない状態です。
どのように修正すればいいでしょうか？
(求めようとしている方法は二分法です。)

Comments

[cateye]

>機能していない状態・・・これコンパイルできないですよね？
・・・コンパイルしたままのコードを貼り付けて下さい。
studio・・・こんなファイル無い！
インデントが無茶苦茶
タグがCになってますがC++に変えて下さい。

[dit.]

同じ質問が2つ投稿されているようです。
回答がつかないうちに片方を削除依頼出しておいた方がよろしいかと。
https://teratail.com/questions/196849

[jimbe]

C++でしたら, タグもC＋＋にしていただけますか.
また, 出力が見当たりません. 解の一つが出ているのでしたら出力があるはずです. それが無いということはご提示の為にソースを変更していることになり, ソースと変更しているということは原因箇所も変更されている可能性があるということになります.
実際に動作されているコードそのままと, どの様に入力したらどの様に出力されるハズがどうなった…という情報の組み合わせが正確でなければ, 正確な回答は得難いかと思います.

Answer 1

#include <cstdio>
#include <cmath>

using namespace std;

int main()
{
    int    i, k, a[2];
    double xl, xr, xc;
    double fl, fr, fc, ss;

    k = 0;
    for (i = 0; i < 10; i++) {
        fl = sin(i) + 1 / sqrt(2);
        fc = sin(i + 1) + 1 / sqrt(2);
        if (fl * fc < 0) {
            a[k] = i;
            k++;
        }
    }
    for (k = 0; k < 2; k++) {
        ss = 1.0;
        xl = a[k];
        xr = a[k] + 1;
        while (ss > 1.0e-10) {
            xc = (xr + xl) / 2;
            fl = sin(xl) + 1 / sqrt(2);
            fr = sin(xr) + 1 / sqrt(2);
            fc = sin(xc) + 1 / sqrt(2);
            if (fl * fc < 0) {
                xr = xc;
            } else {
                xl = xc;
            }
            ss = fabs(1 - xr / xl);
        }
    }

    return 0;
}

これでコンパイルは通りますが・・・どうしたいんでしょう？

Comments

[cateye]

sqrt(2);は事前に計算しておいた方が効率いいですよｄ＾＾