function plot1
x=[0 1 2 3];
y=[2 3 3 5];
s=-1:0.05:4;
f=0.9*s+1.9;
plot(x,y,'*',s,f)