Keyword: 複素多項式の根
概要
本サンプルは複素多項式の根を求めるC言語によるサンプルプログラムです。 本サンプルは以下に示される次数が5の複素多項式の根を求めて出力します。
※本サンプルはnAG Cライブラリに含まれる関数 nag_zeros_complex_poly() のExampleコードです。本サンプル及び関数の詳細情報は nag_zeros_complex_poly のマニュアルページをご参照ください。
ご相談やお問い合わせはこちらまで
入力データ
（本関数の詳細はnag_zeros_complex_poly のマニュアルページを参照）| このデータをダウンロード |
nag_zeros_complex_poly (c02afc) Example Program Data
5
5.0 6.0
30.0 20.0
-0.2 -6.0
50.0 100000.0
-2.0 40.0
10.0 1.0
- １行目はタイトル行で読み飛ばされます。
- ２行目に多項式の次数(n)を指定しています。
- ３から８行目に各項の係数の実部(a().re)と虚部(a().im)を指定しています。
出力結果
（本関数の詳細はnag_zeros_complex_poly のマニュアルページを参照）| この出力例をダウンロード |
nag_zeros_complex_poly (c02afc) Example Program Results Degree of polynomial = 5 Roots of polynomial z = -2.4328e+01 -4.8555e+00 z = 5.2487e+00 +2.2736e+01 z = 1.4653e+01 -1.6569e+01 z = -6.9264e-03 -7.4434e-03 z = 6.5264e-03 +7.4232e-03
- ３行目に入力された多項式の次数が出力されています。
- ７から１１行目に複素多項式の根（実部と虚部）が出力されています。
ソースコード
（本関数の詳細はnag_zeros_complex_poly のマニュアルページを参照）
※本サンプルソースコードはnAG数値計算ライブラリ（Windows, Linux, MAC等に対応）の関数を呼び出します。
サンプルのコンパイル及び実行方法
| このソースコードをダウンロード |
/* nag_zeros_complex_poly (c02afc) Example Program.
*
* CLL6I261D/CLL6I261DL Version.
*
* Copyright 2017 Numerical Algorithms Group.
*
* Mark 26.1, 2017.
*
*/
#include <nag.h>
#include <stdio.h>
#include <nag_stdlib.h>
#include <nagc02.h>
int main(void)
{
Nag_Boolean scale;
Complex *a = 0, *z = 0;
Integer exit_status = 0, i, n;
NagError fail;
INIT_FAIL(fail);
printf("nag_zeros_complex_poly (c02afc) Example Program Results\n");
/* Skip heading in data file */
scanf("%*[^\n]");
scanf("%ld", &n);
if (n > 0) {
if (!(a = nAG_ALLOC(n + 1, Complex)) || !(z = nAG_ALLOC(n, Complex)))
{
printf("Allocation failure\n");
exit_status = -1;
goto END;
}
}
else {
printf("Invalid n.\n");
exit_status = 1;
return exit_status;
}
scale = Nag_TRUE;
for (i = 0; i <= n; i++)
scanf("%lf%lf", &a[i].re, &a[i].im);
/* nag_zeros_complex_poly (c02afc).
* Zeros of a polynomial with complex coefficients
*/
nag_zeros_complex_poly(n, a, scale, z, &fail);
if (fail.code != NE_NOERROR) {
printf("Error from nag_zeros_complex_poly (c02afc).\n%s\n", fail.message);
exit_status = 1;
goto END;
}
printf("\nDegree of polynomial = %4ld\n\n", n);
printf("Roots of polynomial\n\n");
for (i = 0; i < n; ++i)
printf("z = %13.4e %+14.4e\n", z[i].re, z[i].im);
END:
nAG_FREE(a);
nAG_FREE(z);
return exit_status;
}