Bu çalışmada, iç basınca maruz bırakılmış kalın cidarlı fonksiyonel kademeli silindirin statik gerilme analizi sayısal olarak ele alınmıştır. İzotropik malzeme için elastik modülün ve Poisson oranının kalınlık boyunca eksponansiyel olarak değiştiği varsayılmıştır. Bu koşullar altında, geleneksel analitik yöntemlerle çözülemeyen değişken katsayılı sınır değer problemi elde edilir. Bu diferansiyel denklemin çözümü hem Tamamlayıcı Fonksiyonlar Yöntemi (TFY) hem de Chebyshev Pseudospektral Yöntemi ile yapılmıştır. Tamamlayıcı Fonksiyonlar Yöntemi, sınır değer problemini, kolayca çözülebilen başlangıç değer problemi sistemine dönüştürerek çözülmesine imkan veren bir metot iken Chebyshev Pseudospektral Yöntemi diferansiyel matris yaklaşımına dayanan ve bu yaklaşım ile diferansiyel denklemi lineer denklem sistemine dönüştüren bir yöntemdir. Literatürde mevcut olan çözümler, bu çalışmada elde edilen sonuçları doğrulamak için kullanılmıştır. Rasgele seçilmiş homojensizlik parametrelerinin stres ve deplasman dağılımları üzerindeki etkieri gösterilmiştir.
