%首先編寫函數f(x)
function f = myfun (a, x)
f = a(1)*x.^2+a(2)*sin(x)+a(3)*x.^3;

x= [3.6, 7.7, 9.3, 4.1, 8.6, 2.8, 1.3, 7.9, 10.0, 5.4 ];
y= [16.5, 150.6, 263.1, 24.7, 208.5, 9.9, 2.7, 163.9, 325.0, 54.3 ];
a0 = [10, 10, 10];
[a, residual] = lsqcurvefit(@myfun, a0, x, y);