Hello all, I need to solve Heat equation, Cylindrical coordinates with the crank-nicholson numerical scheme and by the ADI method to solve the system pf linear algebric equations. I dont know how to discerte the equation: T_t=1/r* _r(r*k*T_t) + 1/r^2 * _teta(k*T_teta)+ S(T) T_t , _r and T_teta arte the partial diferantial in circle. boundary condition: T_r(r=0,teta)=0 ; T(r=1,teta)=1 init cond: T(t=0, teta)=r The discretization shold be with the crank=Nich schema. Does some have a ready script for the numerical program ( in matlab)?

check out matlabcentral.com

Leo, what shoud I look for at matlab.com ? I didn't find anything usefull for my problem.

You should be looking at http://www.matlabcentral.com. You have to do a bit of searching. Click on the File exchange link. I am giving the link below. You may find a matlab code that matches your requirements. http://www.mathworks.com/matlabcentr...oadCategory.do

Also check out their newsgroup.