diff options
author | Nguyễn Gia Phong <vn.mcsinyx@gmail.com> | 2020-02-16 14:26:55 +0700 |
---|---|---|
committer | Nguyễn Gia Phong <vn.mcsinyx@gmail.com> | 2020-02-16 14:26:55 +0700 |
commit | 82e6cf7d1046d6cee16f7e8b044ec33e7ec6c4b7 (patch) | |
tree | f7b7ae0bce69070c47a1b31a85bd2bc69dfecf09 /usth/MATH2.2/hw/cubic | |
parent | 29d1001e2e21eff289bff23412e284c8b2e44595 (diff) | |
download | cp-82e6cf7d1046d6cee16f7e8b044ec33e7ec6c4b7.tar.gz |
[usth] Numerical Method is MATH2.4
Diffstat (limited to 'usth/MATH2.2/hw/cubic')
-rwxr-xr-x | usth/MATH2.2/hw/cubic | 25 |
1 files changed, 0 insertions, 25 deletions
diff --git a/usth/MATH2.2/hw/cubic b/usth/MATH2.2/hw/cubic deleted file mode 100755 index ca24739..0000000 --- a/usth/MATH2.2/hw/cubic +++ /dev/null @@ -1,25 +0,0 @@ -#!/usr/bin/env octave -disp ("Solve quadratic equation a*x^3 + b*x^2 + c*x + d = 0"); -a = input ("a = "); -while (a == 0) - disp ("a must be nonzero"); - a = input ("a = "); -end -b = input ("b = "); -c = input ("c = "); -d = input ("d = "); - -% Using the trigonometric and hyperbolic solutions on Wikipedia: -% https://en.wikipedia.org/wiki/Cubic_function -p = (3*a*c - b^2) / (3*a^2); -q = (2*b^3 - 9*a*b*c + 27*a^2*d) / (27*a^3); -if (p == 0) - disp ("The given cubic equation has only one distinct solution:"); - x = cbrt(-q) - b/a/3 -else - r = 2 * sqrt(-p/3); - disp ("The three solutions of the given quadratic equation are:"); - for k = 0 : 2 - x = r * cos((acos(3*q/p/r) - 2*pi*k) / 3) - b/a/3 - end -end |