質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

新規登録して質問してみよう
ただいま回答率
85.50%
C

C言語は、1972年にAT&Tベル研究所の、デニス・リッチーが主体となって作成したプログラミング言語です。 B言語の後継言語として開発されたことからC言語と命名。そのため、表記法などはB言語やALGOLに近いとされています。 Cの拡張版であるC++言語とともに、現在世界中でもっとも普及されているプログラミング言語です。

Q&A

0回答

412閲覧

ガウスの消去法の結果がうまく出ない

tusanch6

総合スコア12

C

C言語は、1972年にAT&Tベル研究所の、デニス・リッチーが主体となって作成したプログラミング言語です。 B言語の後継言語として開発されたことからC言語と命名。そのため、表記法などはB言語やALGOLに近いとされています。 Cの拡張版であるC++言語とともに、現在世界中でもっとも普及されているプログラミング言語です。

0グッド

0クリップ

投稿2019/05/29 08:51

前提・実現したいこと

ガウスの消去法の実装を作っています。

最終的な結果が適切に出せませんでした。

該当のソースコード

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}

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

izmktr

2019/05/29 10:15

多分関係ないポイントだけど、mallocが3回でているのにfreeが2回しかでてこないところとか、pA = (double**)malloc(sizeof(double*)*(n+2));が正しいとか、n+2の2やn+1の1はどこからでてきたのか気になります
jimbe

2019/05/29 12:56

> 最終的な結果が適切に出せませんでした。 それで, 何をご質問なのでしょう.
guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

まだ回答がついていません

会員登録して回答してみよう

アカウントをお持ちの方は

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.50%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問