More like this
vin=linspace(0.2,2.8);
current = zeros(1,numel(vin));
for i = 1:numel(vin)
v = vin(i);
if v < 0.5
current(i)=0;
elseif 1.5 > v && v >= 0.5
current(i)=(860*10^-6)*(v-0.5).^2;
elseif 2.5 > v && v>= 1.5
current(i)=(860*10^-6)*(2*(v-0.5)*0.049-(0.049)^2) ;
elseif v >= 2.5
current(i)=(860*10^-6)*(2*(v-0.5)*0-(0)^2) ;
end
end
plot(vin,current)