関数とループに関する質問(ｃ言語）

Question

前提・実現したいこと

ソース内で指定したｔのときのｘをif文の条件でｋにいれる関数を作っていましたが、実行してもなにも出力されませんでした。なぜソースの内容が出力されないのか教えていただけると嬉しいです。

該当のソースコード

#define _USE_MATH_DEFINES
#include <stdio.h>
#include <math.h>

#define Num 10

double calucuation(double k[Num])
{
    int N;
    int c;
    double x;
    int i=100;
    double delta = 1.0e-3;
    double t = (double)i*delta;
    
        for(N=1;N<=Num;N++){
            x = 0.0;
        for (c=1;c<=N;c++) {
            double f = 2*c-1;
            double a = 4/(M_PI*(2*c-1));
            x += a * sin(2.0*M_PI*f*t);
        }
        if(0.0<=t && t<0.5){
          k[N] = x - 1;  
        }else{
          k[N] = x + 1;
        } 
        } 
}

int main(){

    double k[Num];
    int N;

    FILE *fp; 
    fp = fopen("2026_acc.txt","w");

    for (N=1;N<=Num;N++) {
        k[N] = N;
        calucuation(k);
        printf("%d: %lf\n",N,k[N]);
        fprintf(fp,"N 近似精度\n");
        fprintf(fp,"%d %lf\n",N,k[N]);
    }
    fclose(fp);
    return 0;
}

Comments

[tiitoi]

とりあえず calucuation() が返り値が void になっていないのはおかしいのではないでしょうか

[-_.can777]

コメントありがとうございます！具体的にはどういうことでしょうか？

[tiitoi]

calucuation(double k[Num]) が値を返さない関数ならば、
double calucuation(double k[Num])
ではなく、
void calucuation(double k[Num])
とする必要があるかと思います。

Answer 1

double k[Num] と宣言した場合は、有効なインデックスは 0, 1, 2, ..., Num - 1 なので、ループの終了条件を N <= Num としてしまうと、範囲外の k[Num] にアクセスしてしまいます。結果が出力されないのは、これが原因でプログラムが落ちてしまっているのではないでしょうか。

二箇所ある以下を次のように修正すれば結果が出力されました。

- for (N = 1; N <= Num; N++) {
+ for (N = 1; N < Num; N++) {

あと、返り値を返さない関数は void としないと、コンパイラによってはコンパイルがまず通らないです。

- double calucuation(double k[Num])
+ void calucuation(double k[Num])