1 files changed, 38 insertions, 0 deletions
diff --git a/1415/CH1/EX1.4.1/ex1.sce b/1415/CH1/EX1.4.1/ex1.sce
new file mode 100644
index 000000000..660d25bcb
--- /dev/null
+++ b/1415/CH1/EX1.4.1/ex1.sce
@@ -0,0 +1,38 @@
+//Example 1 Page 96
+clc
+clear
+function y1=Y1(t)//function for y=0.5t+8
+    y1=0.5*t+8
+endfunction
+
+function y2=Y2(t)//function for y=0.25+9
+    y2=0.25*t+9
+endfunction
+
+t=[0 2 4 6 8 10 12 14]//t values as given in question
+disp(t,'Year:')//displaying t values
+y=[9 9 10 11 11 12 13 13]//y values as given in question
+disp(y,'observed:')//displaying y values
+y1=Y1(t)//function calling
+disp(y1,'Predicted:')//printing predicted values
+disp(y-y1,'Residual:')//printing residual values
+disp((y-y1)^2,'Residual square')//printing residual square values
+
+t=[0 2 4 6 8 10 12 14]//t values as given in question
+disp(t,'Year:')//displaying t values
+y=[9 9 10 11 11 12 13 13]//y values as given in question
+disp(y,'observed:')//dispalying observed values
+y2=Y2(t)//function valling
+disp(y2,'Predicted:')//printing predicted values
+disp(y-y2,'Residual:')//printing residual values
+disp((y-y2)^2,'Residual square')//printing residual square values
+
+
+plot(t,y1)//plotting graph between t and y
+xtitle('y=0.5t+8','t','y')//naming the axes
+scf;//setting the current graphic window
+plot(t,y2)//plotting the graph between t and y
+xtitle('y=0.25t+9','t','y')//naming the axes
+
+
+