From cb5c3be111bf34e90374ad4a525fe029e10f90ca Mon Sep 17 00:00:00 2001
From: LinlyBoi <libkyy@e.email>
Date: Sun, 14 May 2023 22:00:16 +0300
Subject: [PATCH] second derivative done :DDDDD

---
 newton_differentiation.m | 20 ++++++++++++++++++++
 1 file changed, 20 insertions(+)

diff --git a/newton_differentiation.m b/newton_differentiation.m
index 15cc2ec..2caa7ee 100644
--- a/newton_differentiation.m
+++ b/newton_differentiation.m
@@ -52,3 +52,23 @@ end
 function y = b_newt_diff2(Y,h)
 	y =  ( Y(1) - 2 * Y(2) + Y(3) ) / h^2
 end
+
+
+
+
+
+function y = newt_diff_fr(Y,X,x,h)
+	idx = find(X==x);
+	if (idx == length(X))
+		elems = [X(idx), X(idx-1), X(idx-2)]
+		y = b_newt_diff1(elems,h)
+
+	elseif (idx == 1)
+		elems = [X(idx), X(idx+1), X(idx+2)]
+		y = f_newt_diff1(elems,h)
+	else
+		elems = [X(idx-1), X(idx), X(idx+1)]
+		y = c_newt_diff1(elems,h)
+	endif
+	endfunction
+