octave-3.2.2:1> function y=f(x)
> y(1)=4*(2*(log(sin(x(1)))+1)*(1+2*(cos(x(1)))^4-4*(cos(x(1)))^6)/(cos(x(1)))^4+sin(x(1))+2*sin(x(1))/(cos(x(1)))^2-2+4*(sin(x(1)))^2);
> y(2)=4*(2*(log(sin(x(2)))+1)*(1+2*(cos(x(2)))^4-4*(cos(x(2)))^6)/(cos(x(2)))^4+sin(x(2))+2*sin(x(2))/(cos(x(2)))^2-2+4*(sin(x(2)))^2);
> y(3)=4*(2*(log(sin(x(3)))+1)*(1+2*(cos(x(3)))^4-4*(cos(x(3)))^6)/(cos(x(3)))^4+sin(x(3))+2*sin(x(3))/(cos(x(3)))^2-2+4*(sin(x(3)))^2);
> endfunction;
octave-3.2.2:2> [x, info]=fsolve("f",[1;2;3;1])
x =

   0.41127
   2.73032
   2.98125
   1.00000

info =

  -4.2945e-07   5.3428e-06   4.4409e-16

octave-3.2.2:3> function y=f(x)
> y(1)=4*(2*(log(sin(x(1)))+1)*(1+2*(cos(x(1)))^4-4*(cos(x(1)))^6)/(cos(x(1)))^4+sin(x(1))+2*sin(x(1))/(cos(x(1)))^2-2+4*(sin(x(1)))^2);
> y(2)=4*(2*(log(sin(x(2)))+1)*(1+2*(cos(x(2)))^4-4*(cos(x(2)))^6)/(cos(x(2)))^4+sin(x(2))+2*sin(x(2))/(cos(x(2)))^2-2+4*(sin(x(2)))^2);
> y(3)=4*(2*(log(sin(x(3)))+1)*(1+2*(cos(x(3)))^4-4*(cos(x(3)))^6)/(cos(x(3)))^4+sin(x(3))+2*sin(x(3))/(cos(x(3)))^2-2+4*(sin(x(3)))^2);
> y(1)=4*(2*(log(sin(x(1)))+1)*(1+2*(cos(x(1)))^4-4*(cos(x(1)))^6)/(cos(x(1)))^4+sin(x(1))+2*sin(x(1))/(cos(x(1)))^2-2+4*(sin(x(1)))^2);
> endfunction;
octave-3.2.2:4> [x, info]=fsolve("f",[1;2;3;1])
x =

   0.41127
   2.73032
   2.98125
   1.00000

info =

  -4.2945e-07   5.3428e-06   4.4409e-16

octave-3.2.2:5>