Döngüler ile Taylor Serisi Hesaplama

Merhabalar, şöyle bir sorum var;
Klavyeden girilecek döngü sayısına ve  x= 3*pi/(numaranın son 4 rakamı 1030129999 ise 9999 alınacak)  değeri için

cos(x)=1-x^2/2!+ x^4/4!-x^6/6!+..
sin(x)=x-x^3/3!+ x^5/5!-x^7/7!+..

ve tan(x)=sin(x)/cos(x) değerini hesaplayarak, sin(x), cos(x) ve tan(x) değerlerini ekrana yazdıran programı yazınız.
BENİM YAZDIĞIM PROGRAM ŞU

for s=1:1:D
i1=x^ust1;
j1=factorial(ust1);
c=c+((i1/j1)*((-1)^(s+1)));
ust1=ust1+2;
end
ust2=0;
s=0;
for t=1:1:D
i2=x^ust2;
j2=factorial(ust2+1);
s=s+((i2/j2)*((-1)^(t+1)));
ust2=ust2+2;
end
s=s*x;
t=s/c;
fprintf('cosx = %f\n', j1);
fprintf('sinx = %f\n', j2);
fprintf('tanx = %f\n', j2/j1);

LAKİN SONUÇ ELDE EDEMİYORUM YARDIM EDER MİSİNİZ TEŞEKKÜRLER.