最速下降法matlab程序实现

function [drag,ddrag]=fxy(a0,s0)%要优化的目标函数%返回值： drag，函数在（ a0，s0）点处的值； ddrag ，函数在（ a0，s0，drag）处的梯度（行向量