In this paper analytical and numerical methods fur analysis of solidification of castings are described, and the matrix method, one of numerical method, where the nodal point is designated on the element boundary was adapted. The cooling curve obtained by experimental values, when cast steel (0.29%C. 0.62%Si) was poured into CO₂mold, is compared with that of computed values by exploiting computer (V77-600 Data Proceeding System, UNIVAC). The computed value is nearly approximation to the experimental. But the computed value shows a tendency that is a little higher than the experimental in solid-liquid coexisting temperature ranges and much lower than the experimental after solidification. It is considered to result from the lacks and difficulties of ultimately appropriate adaptation of various physical properties and also air cap between castings and mold.