Something has changed in octave 3.3.51?

[Roessli Bertrand]

Hello,

I have an oct-function that used to work until octave 3.3.51, so I guess 
something changed that I cannot figure out. If somebody can help 
I would be thankful. Note that mkoctfile can compile the function. 

In the cc-code I have :

Matrix Ipp_ij=Matrix(3,3,0);
Matrix Ipm_ij=Matrix(3,3,0);
Matrix polar=Matrix(3,3,0);

that I fill somehow and then return the result as a list :

octave_value_list f_return;

through :

f_return(0)=Ipp_list;
f_return(1)=Ipm_list;
f_return(2)=polar_list;

return f_return;


In the m-file I have :

[I_pp_1d_list,I_pm_1d_list,polar_1d_list]=__polmat(poln_cs,P_1,P_2,Pi_all);
I_pp_1d_list

which gives 

I_pp_1d_list =

error: octave_base_value::print (): wrong type argument `cs-list'


thank you in advance

Bertrand Roessli