I've run the program with "Matlab R2008a" but I've received the next error message:
?? Input argument "poly" is undefined.
Error in ==> odegalerkin at 54
I think the problem is with the poly function. Do you have this poly function or this is the same function from the matlab library ?
Also please check the email because I sent you some questions/proposal regarding this method.