課題でこういうのが出ました。cは習いたてで全然理解できません。内積はすんなり出せたのですが、直積を求める際 segmentation falut と表示されてしまいます。なぜか教えて下さい。
任意の次元数nのベクトルを2個入力させ,作成した関数を用いて内積と直積とを計算し,結果を表示するプログラムを作成せよ. ただし,次元数nの最大値は100とし,それぞれの関数は以下の基準に従うこと.
内積(inner product)
double inner_product ( const double vectA[], const double vectB[], int n)
引数として与えられた実数型配列vectA, vectBをそれぞれ n 次元ベクトルと見なし,内積を算出する.返却値は算出された内積とする.
直積(outer product)
void outer_product ( const double vectA[], const double vectB[], int n, double M[][])
引数として与えられた実数型配列vectA, vectBをそれぞれ n 次元ベクトルと見なし,直積を算出する.算出された行列(テンソル)は引数で与えられた実数型2次元配列 M に格納する. 上記宣言のの部分には正しい数値(オブジェクト形式マクロでも良い)を入れる必要がある.
#include <stdio.h>
double inner_product ( const double vectA[], const double vectB[], int n)
{
int i;
double answer;
for ( i = 1 ; i <= n; i ++ ){
answer += vectA[i] * vectB[i];
}
return answer;
}
void outer_product ( const double vectA[], const double vectB[], int n, double M[][100])
{
int i, j, k;
double vecta[0][100];
double vectb[100][0];
double answer[i][j];
for ( i = 1 ; i <= n ; i++ ) {
vecta[0][i] = vectA[i];
vectb[i][0] = vectB[i];
}
for ( i = 1 ; i <= n ; i++ ) {
for ( j = 0 ; j < n ; j++ ) {
M[i][j] = 0.0 ;
M[i][j] = vecta[0][i] * vectb[i][0];
printf("%.2f", M[i][j]);
}
}
}
int main(void)
{
int n = 100;
double vecta [100];
double vectb [100];
double vectc [100];
double answer[100][100];
double ANSWER[100][100];
int i = 0;
double naiseki;
int j = 0;
int k= 0;
naiseki = 0;
printf("ベクトルの次数は?");
scanf("%d", &n);
printf("ベクトルAを入力して下さい");
for (i = 1; i <= n; i ++){
scanf("%lf", &vecta[i]);
}
printf("ベクトルBを入力して下さい");
for (i = 1; i <= n; i ++){
scanf("%lf", &vectb[i]);
}
printf("ベクトルA = ");
for (i = 1; i <= n; i++){
printf("%.2f ", vecta[i]);
}
printf("\nベクトルB = ");
for (i = 1; i <= n; i++){
printf(" %.2f ", vectb[i]);
}
printf("\n");
naiseki = inner_product ( vecta, vectb, n );
printf("内積:%.2f\n", naiseki);
answer[0][0] = 0;
printf("直積:"); outer_product(vecta, vectb, n, answer);
return 0;
}