diff --git a/src/octave/curve.m b/src/octave/curve.m index 62f1c9f81e19763a511003bfc7b51a47b67f0ae0..9e2583400d6d549ff52952a4641c335715025224 100644 --- a/src/octave/curve.m +++ b/src/octave/curve.m @@ -40,10 +40,10 @@ vecs(1,:) = [0; 0]; % So that we can always calculate a difference vecs(2,:) = [0; 1]; % Start for j = 2:steps+1 - vecs(j+1,:) = duneminimise(A,b,func,vecs(j,:)); - diffs(j,:) = vecs(j+1,:) - vecs(j,:); - line([vecs(j,1) vecs(j+1,1)], ... - [vecs(j,2) vecs(j+1,2)], ... + vecs(j+1,:) = duneminimise(A,b,func,vecs(j,:)); + diffs(j,:) = vecs(j+1,:) - vecs(j,:); + line([vecs(j,1) vecs(j+1,1)], ... + [vecs(j,2) vecs(j+1,2)], ... [duneevaluate(A,b,func,vecs(j,:)') duneevaluate(A,b,func,vecs(j+1,:)')], ... 'color', 'y'); end @@ -52,4 +52,3 @@ axis tight; xlabel x; ylabel y; hold off; - diff --git a/src/octave/zigzag.m b/src/octave/zigzag.m index 28e75be116ffaf02d819fac5b542ea1615bbca67..a6667c16ae6a07ad2a40315d06f8b1766bb6e9f2 100644 --- a/src/octave/zigzag.m +++ b/src/octave/zigzag.m @@ -40,10 +40,10 @@ vecs(1,:) = [ 0; 0]; % So that we can always calculate a difference vecs(2,:) = [279; 0]; % Start for j = 2:steps+1 - vecs(j+1,:) = duneminimise(A,b,func,vecs(j,:)); - diffs(j,:) = vecs(j+1,:) - vecs(j,:); - line([vecs(j,1) vecs(j+1,1)], ... - [vecs(j,2) vecs(j+1,2)], ... + vecs(j+1,:) = duneminimise(A,b,func,vecs(j,:)); + diffs(j,:) = vecs(j+1,:) - vecs(j,:); + line([vecs(j,1) vecs(j+1,1)], ... + [vecs(j,2) vecs(j+1,2)], ... [duneevaluate(A,b,func,vecs(j,:)') duneevaluate(A,b,func,vecs(j+1,:)')], ... 'color', 'y'); end @@ -52,4 +52,3 @@ axis tight; xlabel x; ylabel y; hold off; -