[hejiefabio]

问题已经解决程序如下：
x=0:10:920410;
y=.1989236468e-12*x.^2+.1822893318e-6*x-.5334621633e-17;
plot(x,y);
hold on
x=[0,657610,920410];
y=[0,0.2059,0.3363];
plot(x,y,'*')