a1=-7;a2=2.22;a3=5.15;a4=1.26;
x=1:.1:8.6;
for i=1:size(x)
g(i)=a1+a4*x(i)-(a4/log(a2))*a2^(x(i)-a3);
end
plot(x,g)