Keyword: 複素多項式の根
概要
本サンプルは複素多項式の根を求めるC言語によるサンプルプログラムです。 本サンプルは以下に示される次数が5の複素多項式の根を求めて出力します。
※本サンプルはnAG Cライブラリに含まれる関数 nag_zeros_complex_poly() のExampleコードです。本サンプル及び関数の詳細情報は nag_zeros_complex_poly のマニュアルページをご参照ください。
ご相談やお問い合わせはこちらまで
入力データ
(本関数の詳細はnag_zeros_complex_poly のマニュアルページを参照)1 2 3 4 5 6 7 8
このデータをダウンロード |
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
- 1行目はタイトル行で読み飛ばされます。
- 2行目に多項式の次数(n)を指定しています。
- 3から8行目に各項の係数の実部(a().re)と虚部(a().im)を指定しています。
出力結果
(本関数の詳細はnag_zeros_complex_poly のマニュアルページを参照)1 2 3 4 5 6 7 8 9 10 11
この出力例をダウンロード |
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
- 3行目に入力された多項式の次数が出力されています。
- 7から11行目に複素多項式の根(実部と虚部)が出力されています。
ソースコード
(本関数の詳細はnag_zeros_complex_poly のマニュアルページを参照)
※本サンプルソースコードはnAG数値計算ライブラリ(Windows, Linux, MAC等に対応)の関数を呼び出します。
サンプルのコンパイル及び実行方法
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64
このソースコードをダウンロード |
/* 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; }