1 files changed, 15 insertions, 0 deletions
diff --git a/1670/CH9/EX9.16/9_16.sce b/1670/CH9/EX9.16/9_16.sce
new file mode 100755
index 000000000..3c43c7fcb
--- /dev/null
+++ b/1670/CH9/EX9.16/9_16.sce
@@ -0,0 +1,15 @@
+//Example 9.16
+//Particular Solution
+//Page no. 297
+clc;clear;close;
+
+x=poly(0,'x')
+n=2;
+f=x^(n)-2*x^(n-1)+x^(n-2);
+z=roots(f)
+disp(z,f)
+printf('\t\t\nC.F. = (c1+n*c2) (%g)',z(1))
+A=[1,0,0;8,1,0;12,4,1];
+B=[1;0;0];
+C=inv(A)*B;
+printf('\n\n\t  2\nP.I = (%g)n+(%g)n+(%g)',C(1),C(2),C(3))
+\ No newline at end of file