هفتِ بیجار

[coronaa]

سلام
تابع زیر ظاهرا بدون مشکل کار میکنه. فقط باید دامنه و ثابت ها رو درست انتخاب کرده باشید.

function y=fun_fun(t)

alpha = 101;

a = 0;

k = 1;

h = 0.1;

y=0;

f=@(x) sin(x)

for k=0:round((t-a)/h)

    y=y+(1./(h.^alpha))*((-1).^(k))*nchoosek(alpha,k)*f(t-k*h);

end

end

آلفا نمیشه اعشاری باشه؟
این معادله برای مبحث مرتبه کسری هست و آلفا معمولا یک عدد کسری و بین 0 تا 2 انتخاب میشه.