1 files changed, 49 insertions, 0 deletions
diff --git a/3845/CH2/EX2.16/Ex2_16.sce b/3845/CH2/EX2.16/Ex2_16.sce
new file mode 100644
index 000000000..c397146ff
--- /dev/null
+++ b/3845/CH2/EX2.16/Ex2_16.sce
@@ -0,0 +1,49 @@
+//Example 2.16
+y_0=0;//Initial position (m)
+y=-1;//Position after certain time t (m)
+t=0.45173;//Time (s)
+v_0=0;//Initial velocity (m/s)
+a=((y-y_0)-(v_0*t))/(1/2*t^2);//Acceleration (m/s^2)
+g=-a;//Acceleration due to gravity (m/s^2)
+printf('Acceleration due to gravity = %0.4f m/s^2',g)
+
+///////////////////////////////////////
+//To plot the graphs
+t=[0 0.1 0.2 0.3 0.4 0.5];//Time (s)
+y=[0 -0.049 -0.196 -0.441 -0.784 -1.225];//Position (m)
+v=[0 -0.98 -1.96 -2.94 -3.92 -4.90];//Velocity (m/s)
+a=repmat(-9.80,[1 6]);//Acceleration due to gravity (downwards) (m/s^2)
+
+
+subplot(3,1,1)
+A=gca();//Get the current axes
+A.x_location= "origin";//Set x-axis position
+A.data_bounds=[0,-1.4;0.6,0];//Set data bounds (as seen in the textbook plot)
+A.margins=[0.125,0.125,0.4,0.125];//Adjusting the margins
+plot(t,y,'-rd');//Plotting the graph with a red, solid line with diamond markers at data points
+title('Position vs. Time for Falling Sphere','position',[0.2 0.75]);//Title and its position
+xlabel('Time t (s)','position',[0.275 0.5]);//x-axis label and its position
+ylabel('Position y (m)');//y-axis label
+
+subplot(3,1,2)
+A=gca();//Get the current axes
+A.x_location= "origin";//Set x-axis position
+A.data_bounds=[0,-6;0.6,0];//Set data bounds (as seen in the textbook plot)
+A.margins=[0.125,0.125,0.4,0.125];//Adjusting the margins
+plot(t,v,'-gd');//Plotting the graph with a green, solid line and diamond markers at data points
+title('Velocity vs. Time for Falling Sphere','position',[0.2 3]);//Title and its position
+xlabel('Time t (s)','position',[0.275 2]);//x-axis label and position
+ylabel('Velocity v (m/s)');//y-axis label
+
+subplot(3,1,3)
+A=gca();//Get the current axes
+A.x_location= "origin";//Set x-axis position
+A.data_bounds=[0,-11;0.6,0];//Set data bounds 
+A.margins=[0.125,0.125,0.4,0.1];//Adjusting the margins
+plot(t,a,'-bd');//Plotting the graph with a blue, solid line and diamond markers at data points
+title('Acceleration vs. Time for Falling Sphere','position',[0.19 5]);//Title and its position
+xlabel('Time t (s)','position',[0.275 3]);//x-axis label and its position
+ylabel('Acceleration a (m/s^2)');//y-axis label
+///////////////////////////////////////////////////////
+//Openstax - College Physics
+//Download for free at http://cnx.org/content/col11406/latest