function f input Real t "Time"; input Real x "a point between 0 and l"; input Real a "a=lambda/cp/rho"; input Real l "length"; input Real T_0 "fixed temperature at x=0"; input Real T_l "fixed temperature at x=l"; input Real TK "initial temperature along rod"; input Integer Nu = 200 "number of iterations"; output Real y "output temperature"; end f;