function y=sline(x,a,b)
% sline.m
y=a+b*x;