[usth/MATH2.2] Numerical Methods

The future starts now.

1 files changed, 14 insertions, 0 deletions

diff --git a/usth/MATH2.2/hw/quadratic b/usth/MATH2.2/hw/quadratic
new file mode 100755
index 0000000..b6ea91a
--- /dev/null
+++ b/usth/MATH2.2/hw/quadratic
@@ -0,0 +1,14 @@
+#!/usr/bin/env octave
+disp ("Solve quadratic equation a*x^2 + b*x + c = 0");
+a = input ("a = ");
+while (a == 0)
+  disp ("a must be nonzero");
+  a = input ("a = ");
+end
+b = input ("b = ");
+c = input ("c = ");
+
+sd = sqrt(b*b - a*c*4);
+disp ("The two solutions of the given quadratic equation are:");
+x0 = (-b - sd) / a / 2
+x1 = (-b + sd) / a / 2