1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
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
|
