im 65 years old

Author: Duo-Star

lib/MathForest/Algebra/Trunk/Fertile/QNum.dart

@@ -1,5 +1,6 @@
 //
 import 'dart:math';
+import 'DNum.dart';
 /*
 
  */
@@ -10,15 +11,24 @@ class QNum {
   dynamic n3;
   dynamic n4;
 
-  QNum([dynamic n1, dynamic n2, dynamic n3, dynamic n4]):
-        n1 = n1 ?? -1,
-        n2 = n2 ?? 1,
-        n3 = n3 ?? 2,
-        n4 = n4 ?? 3;
+  QNum([dynamic n1, dynamic n2, dynamic n3, dynamic n4])
+    : n1 = n1 ?? -1,
+      n2 = n2 ?? 1,
+      n3 = n3 ?? 2,
+      n4 = n4 ?? 3;
+
+
+
+  static harmonic(DNum dn, num t) {
+    num a = dn.n1;
+    num b = dn.n2;
+    num m = (t * b - a) / (t - 1);
+    num n = (t * b + a) / (t + 1);
+    return QNum(a, m, b, n);
+  }
 
   @override
   String toString() {
     return 'QNum(${n1.toString()}, ${n2.toString()}, ${n3.toString()}, ${n4.toString()})';
   }
-
-}
+}

lib/MathForest/Geometry/D2/Conic/Circle.dart

@@ -26,8 +26,12 @@ class Circle {
 
   Circle.new2P(Vector p1, Vector p2) : p = p1, r = (p2 - p1).len;
 
+  Circle.newDiameter(DPoint dp) : p = dp.mid, r = dp.len/2;
+
+
   num get area => pi * r * r;
   num get cir => 2 * pi * r;
+  DPoint get F => DPoint.overlap(p);
 
   Vector indexPoint(num theta) => p + Vector(r) * cos(theta) + Vector(0, r) * sin(theta);
   DPoint indexDPoint(DNum theta) => DPoint(indexPoint(theta.n1), indexPoint(theta.n2));
@@ -38,6 +42,15 @@ class Circle {
     indexPoint(theta.n4),
   );
 
+  num thetaClosestP(Vector P) {
+    Vector dp = P-p;
+    return atan2(dp.y, dp.x);
+  }
+
+  Vector closestP(Vector P) {
+    return p + (P-p).unit * r;
+  }
+
   num disP(Vector p0){
     return funcs.abs((p0 - p).len - r);
   }