Several issues:
- new_vec - old_vec is a vector. How do you want to compare it to your (scalar) error margin? (one option: sum the absolute differences). Or, do you want to compare elements of a single vector?
- you don't actually update p anywhere in this code
- use code formatting blocks:
p = [ 1 1 0.5]';
marg= 1e-3;
[d,g]= func(x,u,p)
for i = 1:100
err(i+1) = p(i+1)-p(i).
end
if err <=marg
d(i)= d(i-1) %d and g remains the same as the previous
else
[d,g]= func(x,u,p) % compute a new d and g with respect to the new p values.
end
