No BSD License
xfunction D=deriv3(M); D=zeros(M); for i=0:M-4 for j = i+3:2:M-1 D(i+1,j+1)=(1/myC(i))*(j)*(j^2-i^2)*((j-2)^2-(i-2)^2); end end
Contact us at files@mathworks.com