matlab求解方程系数并代入一个值x求 y

x=[1;2;3];
y=[1;4;7];
p=fittype('ax^2+bx+c', 'independent','x');%定义函数定义自变量
opt=fitoptions(p);
opt.StartPoint=[0.1 0.1 0.1];%设置初始参数
f1=fit(x,y,p,opt)

这个代码可以求出系数但是a的值应该是0 结果却不是应该怎么做啊？还有该怎么实现代入一个x求y呢？