Luận văn thạc sĩ giải số hệ phương trình khuếch tán đồng thời bốn thành phần (b, as, i và v) tìm phân bố nồng độ trong vật liệu silic dựa trên lý thuyết nhiệt động học không thuận nghịch