summaryrefslogtreecommitdiff
path: root/24/CH11/EX11.1/Example11_1.sce
diff options
context:
space:
mode:
Diffstat (limited to '24/CH11/EX11.1/Example11_1.sce')
-rwxr-xr-x24/CH11/EX11.1/Example11_1.sce22
1 files changed, 22 insertions, 0 deletions
diff --git a/24/CH11/EX11.1/Example11_1.sce b/24/CH11/EX11.1/Example11_1.sce
new file mode 100755
index 000000000..3eee34e75
--- /dev/null
+++ b/24/CH11/EX11.1/Example11_1.sce
@@ -0,0 +1,22 @@
+//Given that
+t = poly(0, 't')
+A = -1.00-0.600*t+0.250*t^2
+
+//Sample Problem 11-1a
+printf("**Sample Problem 11-1a**\n")
+Ts = [-3:0.5:6]
+As = horner(A, Ts)
+xset('window', 1)
+xtitle( 'angular variable for the disk v/s time', 'time(sec)', 'Y-axis')
+plot(Ts, As, 'm-o')
+
+//Sample Problem 11-1b
+printf("\n**Sample Problem 11-1b**\n")
+To = roots(derivat(A))
+printf("At t=%fsec, theta approaches its minimum value equal to %f\n", To, horner(A, To))
+
+//Sample Problem 11-1c
+printf("\n**Sample Problem 11-1c**\n")
+Os = horner(derivat(A), Ts)
+plot(Ts, Os, 'r-+')
+legend('theta(rad)', 'omega(rad/s)') \ No newline at end of file