Re: Optimise function

[Nir Krakauer]

How about 

function result = obj(Angletweak)

  .... [obtain D]

  result = D(4, 9)

endfunction

[xmin, fval] = fminsearch (@obj, Angletweak_initial_guess);