added rotating pendulum

wjb11005 · Jan 31, 2018 · a85741e · a85741e
1 parent f514f06
commit a85741e
Show file tree

Hide file tree

Showing 4 changed files with 349 additions and 0 deletions.
diff --git a/ch_03/pend_rot.m b/ch_03/pend_rot.m
@@ -0,0 +1,5 @@
+function dy=pend_rot(t,y,l,w,a)
+  dy=zeros(size(y));
+  dy(1)=y(2);
+  dy(2)=-9.81/l*sin(y(1))+w^2*sin(y(1))*cos(y(1))+a*w^2/l*cos(y(1));
+end
diff --git a/ch_03/rotating_pendulum.m b/ch_03/rotating_pendulum.m
@@ -0,0 +1,10 @@
+
+l=0.3;
+a=1;
+w=10;
+
+my_ode=@(t,y) pend_rot(t,y,l,w,a)
+
+[t,y]=ode23(my_ode,[0 2],[pi/6,0]);
+
+plot(t,y(:,1))
diff --git a/octave-workspace b/octave-workspace