Simple routine won't run.

[David Mebane]

Hello

Have written a (very first) octave routine:

[begin code]
function fun=ent3(d)

x31=0;
for i = 0:1
for j = 0:1
for k = 0:1
x31
+=(d^i*(1-d)^(1-i)*(d/2)^(k+j)*(1-d/2)^(2-k-j))*log(d^i*(1-d)^(1-i)*(d/2)^(k+j)*(1-d/2)^(2-k-j));
endfor
endfor
endfor

x21=0;
for i=0:1
for j=0:1
x21 +=
(d^i*(1-d)(1-i)*(d/2)^j*(1-d/2)^(1-j))log(d^i*(1-d)(1-i)*(d/2)^j*(1-d/2)^(1-j));
endfor
endfor

x22=0;
for j=0:1
for k=0:1
x22 += (d/2)^(k+j)*(1-d/2)^(2-k-j)log((d/2)^(k+j)*(1-d/2)^(2-k-j));
endfor
endfor

fun = -12*x31 + 16*x21 + 6*x22 - 5(d*log(d) + (1-d)*log(1-d));

endfunction
[end code]

Running this gets me:

[begin octave]
octave:5> ent3(0.1)
parse error near line 16 of file /home/dm215/Desktop/ent3.m

  syntax error

>>> x21 +=
(d^i*(1-d)(1-i)*(d/2)^j*(1-d/2)^(1-j))log(d^i*(1-d)(1-i)*(d/2)^j*(1-d/2)^(1-j));
                                        ^

error: `ent3' undefined near line 5 column 1
[end octave]

which I don't completely understand.  Any thoughts?


Thanks,
-David