B=5:0.02:12;
z=zeros(size(B));
i=0;
for b3=B
	i=i+1;
	b=[15,12,b3];
	z(i)=prodf(b);
end;
plot(B,z);
