function y=yrsm(x,b0,b,BF)
% MAKE SURE YOU USE FULL B (not B from leasqr1 =lower B)

%Arthur Jutan 
%Dept of Chemical and Biochemical Enginering
%Univ of Western Ontario
%London Ontario Canada, N6A 5B9

xc=x';

if nargin > 1 , b0_=b0;b_=b;B_=BF; end
y=b0_ + xc*b_' +xc*B_*xc'; % calc y for rsm model
