前提・実現したいこと
ガウスの消去法の実装を作っています。
最終的な結果が適切に出せませんでした。
該当のソースコード
C
1#include <stdio.h> 2#include <stdlib.h> 3 4void Gauss(); 5void PrintAnswer(); 6double *ans; 7static int n; 8 9int main(void) { 10 11 printf("Input the number of dimention\n"); 12 scanf_s(" %d", &n); 13 14 //ガウスの除去法 15 Gauss(); 16 17 return 0; 18} 19 20void Gauss() { 21 22 int i, j, k; 23 double f; 24 double **pA; 25 double *pB; 26 27 pA = (double**)malloc(sizeof(double)*(n+2)); 28 for (i = 0;i < n;i++) { 29 pA[i] = (double *)malloc((n + 1) * sizeof(double)); 30 } 31 pB = (double*)malloc(sizeof(double)*(n+2)); 32 33 printf("係数を入力してください\n"); 34 for (i = 0; i < n; i++) { 35 for (j = 0; j < n; j++) { 36 printf("a[%d][%d]:",i,j); 37 scanf_s(" %lf", &pA[i][j]); 38 } 39 } 40 41 printf("解を入力してください\n"); 42 43 for (i = 0; i < n; i++) { 44 printf("%d行目の解:",i); 45 scanf_s(" %lf", &pB[i]); 46 } 47 48 for (k = 0; k < n - 1; k++) 49 for (i = k + 1; i < n; i++) { 50 f = -pA[i][k] / pA[k][k]; 51 for (j = k + 1; j < n; j++) { 52 pA[i][j] += f * pA[k][j]; 53 } 54 pB[i] += f * pB[k]; 55 56 } 57 // backward substitution 58 for (i = n - 1; i >= 0; i--) { 59 for (j = i + 1; j < n; j++) 60 pB[i] -= pA[i][j] * pB[j]; 61 pB[i] = pB[i] / pA[i][i]; 62 } 63 64 PrintAnswer(pB); 65 66 free(pA); 67 free(pB); 68} 69 70void PrintAnswer(double *ans) { 71 72 int i; 73 74 printf("Answer:["); 75 for (i = 0;i < n;i++) { 76 printf("%lf", &ans[i]); 77 if (i != n - 1) { 78 printf(","); 79 } 80 } 81 printf("] "); 82}
あなたの回答
tips
プレビュー