خطای قبلی رو هم برطرف کردم و کد جدید بصورت زیر شد:
clc;
clear;
close all;
alpha = 0.5;
a = 0;
k = 1;
t = 10;
h = 0.1;
upper_band = round((t-a)/h);
expr = (1./(h.^alpha))*symsum(((-1).^(k))*nchoosek(alpha,k)*f(t-k*h), k , 0 ,upper_band);
f(t) = limit(expr,h,0)
و الآن خطای جدید روی f هست.
Undefined function 'f' for input arguments of type 'double'.
Error in letnikov (line 12)
expr = (1./(h.^alpha))*symsum(((-1).^(k))*nchoosek(alpha,k)*f(t-k*h), k , 0
,upper_band);