110 lines
32 KiB
Plaintext
110 lines
32 KiB
Plaintext
clc;clear;clf;rng(0);tic;set(0,'DefaultFigureWindowStyle','docked');
|
|
O_EGAMI_YARA_O_ARAY_IMAGE_O=...
|
|
[...
|
|
1,1,1,1,1,1,1,1,1,1,1,1,1,1,1;...
|
|
1,1,1,1,1,1,1,1,1,1,1,1,1,1,1;...
|
|
1,1,1,1,0,0,0,0,0,0,0,1,1,1,1;...
|
|
1,1,1,0,0,0,0,0,0,0,0,0,1,1,1;...
|
|
1,1,0,0,0,1,1,1,1,1,0,0,0,1,1;...
|
|
1,1,0,0,1,1,1,1,1,1,1,0,0,1,1;...
|
|
1,1,0,0,1,1,0,0,0,1,1,0,0,1,1;...
|
|
1,1,0,0,1,1,0,1,0,1,1,0,0,1,1;...
|
|
1,1,0,0,1,1,0,0,0,1,1,0,0,1,1;...
|
|
1,1,0,0,1,1,1,1,1,1,1,0,0,1,1;...
|
|
1,1,0,0,0,1,1,1,1,1,0,0,0,1,1;...
|
|
1,1,1,0,0,0,0,0,0,0,0,0,1,1,1;...
|
|
1,1,1,1,0,0,0,0,0,0,0,1,1,1,1;...
|
|
1,1,1,1,1,1,1,1,1,1,1,1,1,1,1;...
|
|
1,1,1,1,1,1,1,1,1,1,1,1,1,1,1;...
|
|
]...
|
|
%{
|
|
[...
|
|
0,0,0;...
|
|
0,1,0;...
|
|
0,0,0;...
|
|
]...
|
|
%}
|
|
;
|
|
|
|
O_EGAMI_TUPNI_O_INPUT_IMAGE_O=O_EGAMI_YARA_O_ARAY_IMAGE_O;
|
|
|
|
O_OITAR_ELACS_O_SCALE_RATIO_O=3^(4);
|
|
|
|
O_EZISER_SELNOITALOPRETNI_O_INTERPOLATIONLES_RESIZE_O=imresize(O_EGAMI_TUPNI_O_INPUT_IMAGE_O,O_OITAR_ELACS_O_SCALE_RATIO_O,'NEAREST');
|
|
O_EZISER_NOITALOPRETNI_ERAUQS_O_SQUARE_INTERPOLATION_RESIZE_O=imresize(O_EGAMI_TUPNI_O_INPUT_IMAGE_O,O_OITAR_ELACS_O_SCALE_RATIO_O,{@O_NOITALOPRETNI_ERAUQS_O_SQUARE_INTERPOLATION_O,1});
|
|
O_EZISER_NOITALOPRETNI_RAENIL_O_LINEAR_INTERPOLATION_RESIZE_O=imresize(O_EGAMI_TUPNI_O_INPUT_IMAGE_O,O_OITAR_ELACS_O_SCALE_RATIO_O,{@O_NOITALOPRETNI_RAENIL_O_LINEAR_INTERPOLATION_O,1});
|
|
O_EZISER_NOITALOPRETNI_CISOC_O_COSINE_INTERPOLATION_RESIZE_O=imresize(O_EGAMI_TUPNI_O_INPUT_IMAGE_O,O_OITAR_ELACS_O_SCALE_RATIO_O,{@O_NOITALOPRETNI_CISOC_O_COSINE_INTERPOLATION_O,1});
|
|
O_EZISER_NOITALOPRETNI_SUIBAF_O_FABIUS_INTERPOLATION_RESIZE_O=imresize(O_EGAMI_TUPNI_O_INPUT_IMAGE_O,O_OITAR_ELACS_O_SCALE_RATIO_O,{@O_NOITALOPRETNI_8SUIBAF_O_FABIUS8_INTERPOLATION_O,1});
|
|
O_EZISER_NOITALOPRETNI_CRA_O_ARC_INTERPOLATION_RESIZE_O=imresize(O_EGAMI_TUPNI_O_INPUT_IMAGE_O,O_OITAR_ELACS_O_SCALE_RATIO_O,{@O_NOITALOPRETNI_CRA_O_ARC_INTERPOLATION_O,1});
|
|
O_EZISER_NOITALOPRETNI_CIBUCIB_O_BICUBIC_INTERPOLATION_RESIZE_O=imresize(O_EGAMI_TUPNI_O_INPUT_IMAGE_O,O_OITAR_ELACS_O_SCALE_RATIO_O,'BICUBIC');
|
|
|
|
fprintf('%d ',size(O_EZISER_SELNOITALOPRETNI_O_INTERPOLATIONLES_RESIZE_O));fprintf('%s ',upper(class(O_EZISER_SELNOITALOPRETNI_O_INTERPOLATIONLES_RESIZE_O)));fprintf('%d ',whos('O_EZISER_SELNOITALOPRETNI_O_INTERPOLATIONLES_RESIZE_O').bytes);
|
|
|
|
O_SGNIZISER_O_RESIZINGS_O={ ...
|
|
O_EZISER_NOITALOPRETNI_CIBUCIB_O_BICUBIC_INTERPOLATION_RESIZE_O,...
|
|
O_EZISER_NOITALOPRETNI_CRA_O_ARC_INTERPOLATION_RESIZE_O,...
|
|
O_EZISER_NOITALOPRETNI_SUIBAF_O_FABIUS_INTERPOLATION_RESIZE_O,...
|
|
O_EZISER_NOITALOPRETNI_CISOC_O_COSINE_INTERPOLATION_RESIZE_O,...
|
|
O_EZISER_NOITALOPRETNI_RAENIL_O_LINEAR_INTERPOLATION_RESIZE_O,...
|
|
O_EZISER_NOITALOPRETNI_ERAUQS_O_SQUARE_INTERPOLATION_RESIZE_O,...
|
|
O_EZISER_SELNOITALOPRETNI_O_INTERPOLATIONLES_RESIZE_O,...
|
|
O_EZISER_NOITALOPRETNI_ERAUQS_O_SQUARE_INTERPOLATION_RESIZE_O,...
|
|
O_EZISER_NOITALOPRETNI_RAENIL_O_LINEAR_INTERPOLATION_RESIZE_O,...
|
|
O_EZISER_NOITALOPRETNI_CISOC_O_COSINE_INTERPOLATION_RESIZE_O,...
|
|
O_EZISER_NOITALOPRETNI_SUIBAF_O_FABIUS_INTERPOLATION_RESIZE_O,...
|
|
O_EZISER_NOITALOPRETNI_CRA_O_ARC_INTERPOLATION_RESIZE_O,...
|
|
O_EZISER_NOITALOPRETNI_CIBUCIB_O_BICUBIC_INTERPOLATION_RESIZE_O ...
|
|
};
|
|
|
|
O_SNOITALOPRETNI_O_INTERPOLATIONS_O={ ...
|
|
@(X) O_NOITALOPRETNI_CIBUCIB_O_BICUBIC_INTERPOLATION_O(X),...
|
|
@(X) O_NOITALOPRETNI_CRA_O_ARC_INTERPOLATION_O(X),...
|
|
@(X) O_NOITALOPRETNI_8SUIBAF_O_FABIUS8_INTERPOLATION_O(X),...
|
|
@(X) O_NOITALOPRETNI_CISOC_O_COSINE_INTERPOLATION_O(X),...
|
|
@(X) O_NOITALOPRETNI_RAENIL_O_LINEAR_INTERPOLATION_O(X),...
|
|
@(X) O_NOITALOPRETNI_ERAUQS_O_SQUARE_INTERPOLATION_O(X),...
|
|
@(X) O_SELNOITALOPRETNI_O_INTERPOLATIONLES_O(X),...
|
|
@(X) O_NOITALOPRETNI_ERAUQS_O_SQUARE_INTERPOLATION_O(X),...
|
|
@(X) O_NOITALOPRETNI_RAENIL_O_LINEAR_INTERPOLATION_O(X),...
|
|
@(X) O_NOITALOPRETNI_CISOC_O_COSINE_INTERPOLATION_O(X),...
|
|
@(X) O_NOITALOPRETNI_8SUIBAF_O_FABIUS8_INTERPOLATION_O(X),...
|
|
@(X) O_NOITALOPRETNI_CRA_O_ARC_INTERPOLATION_O(X),...
|
|
@(X) O_NOITALOPRETNI_CIBUCIB_O_BICUBIC_INTERPOLATION_O(X) ...
|
|
};
|
|
|
|
O_YTITNAUQ_EGAMI_O_IMAGE_QUANTITY_O=length(O_SGNIZISER_O_RESIZINGS_O);
|
|
set (gcf,'Name', 'O____STOLP_NOITALOPRETNI_HTIW_SGNIZISER____O____RESIZINGS_WITH_INTERPOLATION_PLOTS____O')
|
|
tiledlayout(2,O_YTITNAUQ_EGAMI_O_IMAGE_QUANTITY_O,'TILESPACING','NONE','PADDING','NONE');
|
|
for i=1:O_YTITNAUQ_EGAMI_O_IMAGE_QUANTITY_O nexttile;imshow(O_SGNIZISER_O_RESIZINGS_O{i});end
|
|
for i=1:O_YTITNAUQ_EGAMI_O_IMAGE_QUANTITY_O nexttile;X=-2.5:1/2^8:2.5;Y=O_SNOITALOPRETNI_O_INTERPOLATIONS_O{i}(X);plot(X,Y,'COLOR','#9C9C9C');pbaspect([1 1 1]);O_SEXA_O_AXES_O=gca;O_SEXA_O_AXES_O.XColor='NONE';O_SEXA_O_AXES_O.YColor='NONE';
|
|
O_NIGRAM_X_O_X_MARGIN_O=8/256;O_NIGRAM_Y_O_Y_MARGIN_O=8/256;xlim([min(X)-O_NIGRAM_X_O_X_MARGIN_O,max(X)+O_NIGRAM_X_O_X_MARGIN_O]);ylim([min(Y)-O_NIGRAM_Y_O_Y_MARGIN_O,max(Y)+O_NIGRAM_Y_O_Y_MARGIN_O]);
|
|
O_SKCIT_X_O_X_TICKS_O=get(gca,'XTICK');O_SKCIT_Y_O_Y_TICKS_O=get(gca,'YTICK');
|
|
|
|
for i=1:length(O_SKCIT_X_O_X_TICKS_O) text(O_SKCIT_X_O_X_TICKS_O(i),max(Y)+((8/256)*(max(Y)-min(Y))),num2str(O_SKCIT_X_O_X_TICKS_O(i)),'HORIZONTALALIGNMENT','CENTER','VERTICALALIGNMENT','TOP','FONTSIZE',7,'COLOR','#9C9C9C');end
|
|
for i=1:length(O_SKCIT_Y_O_Y_TICKS_O) text(max(X)-((1.5/256)*(max(X)-min(X))),O_SKCIT_Y_O_Y_TICKS_O(i),num2str(O_SKCIT_Y_O_Y_TICKS_O(i)),'HORIZONTALALIGNMENT','RIGHT','VERTICALALIGNMENT','MIDDLE','FONTSIZE',7,'COLOR','#9C9C9C');end
|
|
for i=1:length(O_SKCIT_X_O_X_TICKS_O) text(O_SKCIT_X_O_X_TICKS_O(i),min(Y)-((8/256)*(max(Y)-min(Y))),num2str(O_SKCIT_X_O_X_TICKS_O(i)),'HORIZONTALALIGNMENT','CENTER','VERTICALALIGNMENT','BOTTOM','FONTSIZE',7,'COLOR','#9C9C9C');end
|
|
for i=1:length(O_SKCIT_Y_O_Y_TICKS_O) text(min(X)+((1.5/256) *(max(X)-min(X))),O_SKCIT_Y_O_Y_TICKS_O(i),num2str(O_SKCIT_Y_O_Y_TICKS_O(i)),'HorizontalAlignment','left','VerticalAlignment','middle','FontSize',7,'Color','#9C9C9C');end
|
|
end
|
|
|
|
imwrite(O_EZISER_NOITALOPRETNI_SUIBAF_O_FABIUS_INTERPOLATION_RESIZE_O,'C:\GNP.O____TUPTUO_EZISER_NOITALOPRETNI_MOTSUC____O____CUSTOM_INTERPOLATION_RESIZE_OUTPUT____O.PNG','png','BitDepth',16,'Alpha',1-O_EZISER_NOITALOPRETNI_SUIBAF_O_FABIUS_INTERPOLATION_RESIZE_O)
|
|
exrwrite(O_EZISER_NOITALOPRETNI_SUIBAF_O_FABIUS_INTERPOLATION_RESIZE_O*2^32-1,'C:\rxe.O____TUPTUO_EZISER_NOITALOPRETNI_MOTSUC____O____CUSTOM_INTERPOLATION_RESIZE_OUTPUT____O.exr','OutputType','uint32','Alpha',(1-O_EZISER_NOITALOPRETNI_SUIBAF_O_FABIUS_INTERPOLATION_RESIZE_O)*2^32-1)
|
|
|
|
fprintf(1,['%.' num2str(2^6) 'f.\n'],toc);
|
|
|
|
function f=O_SELNOITALOPRETNI_O_INTERPOLATIONLES_O(X),f=(abs(X/.5)<1);end
|
|
function f=O_NOITALOPRETNI_ERAUQS_O_SQUARE_INTERPOLATION_O(X),f=(abs(X/.75)<1);end
|
|
function f=O_NOITALOPRETNI_RAENIL_O_LINEAR_INTERPOLATION_O(X),f=(1-abs(X/1)).*(abs(X)<=1);end
|
|
function f=O_NOITALOPRETNI_CISOC_O_COSINE_INTERPOLATION_O(X),f=(.5+.5*cos(X*1*(4*atan(1)))).*(abs(X)<=1);end
|
|
function f=O_NOITALOPRETNI_CRA_O_ARC_INTERPOLATION_O(X),M=1.5;f=(.5-.5.*((-1).^floor(((X.*M)-.5).*2/2).*(1-abs(mod(((X.*M)-.5).*2,2)-1).^2).^(1/2))).*(abs((X.*M))<1);end
|
|
function f=O_NOITALOPRETNI_CIBUCIB_O_BICUBIC_INTERPOLATION_O(X),M=256/256;f=(1.5*abs(X).^3-2.5*abs(X).^2+1).*(abs(X)<=1)+(-0.5*abs(X).^3+2.5*abs(X).^2-4*abs(X)+2).*((1<abs(X))&(abs(X)<=2));end
|
|
function f=O_NOITALOPRETNI_TNENOPXE_O_EXPONENT_INTERPOLATION_O(X),M=256/256;f = (1-(0.5 - 0.5.*(-1).^floor(M.*X) + (-1).^floor(M.*X)./(1+exp((-1+mod(M.*X,1)).^(-1)+mod(M.*X,1).^(-1))))).*(abs(X)<=1);end
|
|
function f=O_NOITALOPRETNI_1SUIBAF_O_FABIUS1_INTERPOLATION_O(X),M=256/256;f=(0.5-(-1).^floor(M.*X).*(-0.5+2.*(-max(0,-3/4+mod(M.*X,1))+max(0,-1/4+mod(M.*X,1))))).*(abs(X)<=1);end
|
|
function f=O_NOITALOPRETNI_2SUIBAF_O_FABIUS2_INTERPOLATION_O(X),M=256/256;f=(0.5-(-1).^floor(M.*X).*(-0.5+4.*(-2.*(-1/2.*max(0,-7/8+mod(M.*X,1)).^2+max(0,-3/8+mod(M.*X,1)).^2/2)+2.*(-1/2.*max(0,-5/8+mod(M.*X,1)).^2+max(0,-1/8+mod(M.*X,1)).^2/2)))).*(abs(X)<=1);end
|
|
function f=O_NOITALOPRETNI_3SUIBAF_O_FABIUS3_INTERPOLATION_O(X),M=256/256;f=(0.5-(-1).^floor(M.*X).*(-0.5+8.*(-4.*(-2.*(-1/6.*max(0,-15/16+mod(M.*X,1)).^3+max(0,-7/16+mod(M.*X,1)).^3/6)+2.*(-1/6.*max(0,-11/16+mod(M.*X,1)).^3+max(0,-3/16+mod(M.*X,1)).^3/6))+4.*(-2.*(-1/6.*max(0,-13/16+mod(M.*X,1)).^3+max(0,-5/16+mod(M.*X,1)).^3/6)+2.*(-1/6.*max(0,-9/16+mod(M.*X,1)).^3+max(0,-1/16+mod(M.*X,1)).^3/6))))).*(abs(X)<=1);end
|
|
function f=O_NOITALOPRETNI_4SUIBAF_O_FABIUS4_INTERPOLATION_O(X),M=256/256;f=(0.5-(-1).^floor(M.*X).*(-0.5+16.*(-8.*(-4.*(-2.*(-1/24.*max(0,-31/32+mod(M.*X,1)).^4+max(0,-15/32+mod(M.*X,1)).^4/24)+2.*(-1/24.*max(0,-23/32+mod(M.*X,1)).^4+max(0,-7/32+mod(M.*X,1)).^4/24))+4.*(-2.*(-1/24.*max(0,-27/32+mod(M.*X,1)).^4+max(0,-11/32+mod(M.*X,1)).^4/24)+2.*(-1/24.*max(0,-19/32+mod(M.*X,1)).^4+max(0,-3/32+mod(M.*X,1)).^4/24)))+8.*(-4.*(-2.*(-1/24.*max(0,-29/32+mod(M.*X,1)).^4+max(0,-13/32+mod(M.*X,1)).^4/24)+2.*(-1/24.*max(0,-21/32+mod(M.*X,1)).^4+max(0,-5/32+mod(M.*X,1)).^4/24))+4.*(-2.*(-1/24.*max(0,-25/32+mod(M.*X,1)).^4+max(0,-9/32+mod(M.*X,1)).^4/24)+2.*(-1/24.*max(0,-17/32+mod(M.*X,1)).^4+max(0,-1/32+mod(M.*X,1)).^4/24)))))).*(abs(X)<=1);end
|
|
function f=O_NOITALOPRETNI_5SUIBAF_O_FABIUS5_INTERPOLATION_O(X),M=256/256;f=(0.5-(-1).^floor(M.*X).*(-0.5+32.*(-16.*(-8.*(-4.*(-2.*(-1/120.*max(0,-63/64+mod(M.*X,1)).^5+max(0,-31/64+mod(M.*X,1)).^5/120)+2.*(-1/120.*max(0,-47/64+mod(M.*X,1)).^5+max(0,-15/64+mod(M.*X,1)).^5/120))+4.*(-2.*(-1/120.*max(0,-55/64+mod(M.*X,1)).^5+max(0,-23/64+mod(M.*X,1)).^5/120)+2.*(-1/120.*max(0,-39/64+mod(M.*X,1)).^5+max(0,-7/64+mod(M.*X,1)).^5/120)))+8.*(-4.*(-2.*(-1/120.*max(0,-59/64+mod(M.*X,1)).^5+max(0,-27/64+mod(M.*X,1)).^5/120)+2.*(-1/120.*max(0,-43/64+mod(M.*X,1)).^5+max(0,-11/64+mod(M.*X,1)).^5/120))+4.*(-2.*(-1/120.*max(0,-51/64+mod(M.*X,1)).^5+max(0,-19/64+mod(M.*X,1)).^5/120)+2.*(-1/120.*max(0,-35/64+mod(M.*X,1)).^5+max(0,-3/64+mod(M.*X,1)).^5/120))))+16.*(-8.*(-4.*(-2.*(-1/120.*max(0,-61/64+mod(M.*X,1)).^5+max(0,-29/64+mod(M.*X,1)).^5/120)+2.*(-1/120.*max(0,-45/64+mod(M.*X,1)).^5+max(0,-13/64+mod(M.*X,1)).^5/120))+4.*(-2.*(-1/120.*max(0,-53/64+mod(M.*X,1)).^5+max(0,-21/64+mod(M.*X,1)).^5/120)+2.*(-1/120.*max(0,-37/64+mod(M.*X,1)).^5+max(0,-5/64+mod(M.*X,1)).^5/120)))+8.*(-4.*(-2.*(-1/120.*max(0,-57/64+mod(M.*X,1)).^5+max(0,-25/64+mod(M.*X,1)).^5/120)+2.*(-1/120.*max(0,-41/64+mod(M.*X,1)).^5+max(0,-9/64+mod(M.*X,1)).^5/120))+4.*(-2.*(-1/120.*max(0,-49/64+mod(M.*X,1)).^5+max(0,-17/64+mod(M.*X,1)).^5/120)+2.*(-1/120.*max(0,-33/64+mod(M.*X,1)).^5+max(0,-1/64+mod(M.*X,1)).^5/120))))))).*(abs(X)<=1);end
|
|
function f=O_NOITALOPRETNI_6SUIBAF_O_FABIUS6_INTERPOLATION_O(X),M=256/256;f=(0.5-(-1).^floor(M.*X).*(-0.5+64.*(-32.*(-16.*(-8.*(-4.*(-2.*(-1/720.*max(0,-127/128+mod(M.*X,1)).^6+max(0,-63/128+mod(M.*X,1)).^6/720)+2.*(-1/720.*max(0,-95/128+mod(M.*X,1)).^6+max(0,-31/128+mod(M.*X,1)).^6/720))+4.*(-2.*(-1/720.*max(0,-111/128+mod(M.*X,1)).^6+max(0,-47/128+mod(M.*X,1)).^6/720)+2.*(-1/720.*max(0,-79/128+mod(M.*X,1)).^6+max(0,-15/128+mod(M.*X,1)).^6/720)))+8.*(-4.*(-2.*(-1/720.*max(0,-119/128+mod(M.*X,1)).^6+max(0,-55/128+mod(M.*X,1)).^6/720)+2.*(-1/720.*max(0,-87/128+mod(M.*X,1)).^6+max(0,-23/128+mod(M.*X,1)).^6/720))+4.*(-2.*(-1/720.*max(0,-103/128+mod(M.*X,1)).^6+max(0,-39/128+mod(M.*X,1)).^6/720)+2.*(-1/720.*max(0,-71/128+mod(M.*X,1)).^6+max(0,-7/128+mod(M.*X,1)).^6/720))))+16.*(-8.*(-4.*(-2.*(-1/720.*max(0,-123/128+mod(M.*X,1)).^6+max(0,-59/128+mod(M.*X,1)).^6/720)+2.*(-1/720.*max(0,-91/128+mod(M.*X,1)).^6+max(0,-27/128+mod(M.*X,1)).^6/720))+4.*(-2.*(-1/720.*max(0,-107/128+mod(M.*X,1)).^6+max(0,-43/128+mod(M.*X,1)).^6/720)+2.*(-1/720.*max(0,-75/128+mod(M.*X,1)).^6+max(0,-11/128+mod(M.*X,1)).^6/720)))+8.*(-4.*(-2.*(-1/720.*max(0,-115/128+mod(M.*X,1)).^6+max(0,-51/128+mod(M.*X,1)).^6/720)+2.*(-1/720.*max(0,-83/128+mod(M.*X,1)).^6+max(0,-19/128+mod(M.*X,1)).^6/720))+4.*(-2.*(-1/720.*max(0,-99/128+mod(M.*X,1)).^6+max(0,-35/128+mod(M.*X,1)).^6/720)+2.*(-1/720.*max(0,-67/128+mod(M.*X,1)).^6+max(0,-3/128+mod(M.*X,1)).^6/720)))))+32.*(-16.*(-8.*(-4.*(-2.*(-1/720.*max(0,-125/128+mod(M.*X,1)).^6+max(0,-61/128+mod(M.*X,1)).^6/720)+2.*(-1/720.*max(0,-93/128+mod(M.*X,1)).^6+max(0,-29/128+mod(M.*X,1)).^6/720))+4.*(-2.*(-1/720.*max(0,-109/128+mod(M.*X,1)).^6+max(0,-45/128+mod(M.*X,1)).^6/720)+2.*(-1/720.*max(0,-77/128+mod(M.*X,1)).^6+max(0,-13/128+mod(M.*X,1)).^6/720)))+8.*(-4.*(-2.*(-1/720.*max(0,-117/128+mod(M.*X,1)).^6+max(0,-53/128+mod(M.*X,1)).^6/720)+2.*(-1/720.*max(0,-85/128+mod(M.*X,1)).^6+max(0,-21/128+mod(M.*X,1)).^6/720))+4.*(-2.*(-1/720.*max(0,-101/128+mod(M.*X,1)).^6+max(0,-37/128+mod(M.*X,1)).^6/720)+2.*(-1/720.*max(0,-69/128+mod(M.*X,1)).^6+max(0,-5/128+mod(M.*X,1)).^6/720))))+16.*(-8.*(-4.*(-2.*(-1/720.*max(0,-121/128+mod(M.*X,1)).^6+max(0,-57/128+mod(M.*X,1)).^6/720)+2.*(-1/720.*max(0,-89/128+mod(M.*X,1)).^6+max(0,-25/128+mod(M.*X,1)).^6/720))+4.*(-2.*(-1/720.*max(0,-105/128+mod(M.*X,1)).^6+max(0,-41/128+mod(M.*X,1)).^6/720)+2.*(-1/720.*max(0,-73/128+mod(M.*X,1)).^6+max(0,-9/128+mod(M.*X,1)).^6/720)))+8.*(-4.*(-2.*(-1/720.*max(0,-113/128+mod(M.*X,1)).^6+max(0,-49/128+mod(M.*X,1)).^6/720)+2.*(-1/720.*max(0,-81/128+mod(M.*X,1)).^6+max(0,-17/128+mod(M.*X,1)).^6/720))+4.*(-2.*(-1/720.*max(0,-97/128+mod(M.*X,1)).^6+max(0,-33/128+mod(M.*X,1)).^6/720)+2.*(-1/720.*max(0,-65/128+mod(M.*X,1)).^6+max(0,-1/128+mod(M.*X,1)).^6/720)))))))).*(abs(X)<=1);end
|
|
function f=O_NOITALOPRETNI_7SUIBAF_O_FABIUS7_INTERPOLATION_O(X),M=256/256;f=(0.5-(-1).^floor(M.*X).*(-0.5+128.*(-64.*(-32.*(-16.*(-8.*(-4.*(-2.*(-1/5040.*max(0,-255/256+mod(M.*X,1)).^7+max(0,-127/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-191/256+mod(M.*X,1)).^7+max(0,-63/256+mod(M.*X,1)).^7/5040))+4.*(-2.*(-1/5040.*max(0,-223/256+mod(M.*X,1)).^7+max(0,-95/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-159/256+mod(M.*X,1)).^7+max(0,-31/256+mod(M.*X,1)).^7/5040)))+8.*(-4.*(-2.*(-1/5040.*max(0,-239/256+mod(M.*X,1)).^7+max(0,-111/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-175/256+mod(M.*X,1)).^7+max(0,-47/256+mod(M.*X,1)).^7/5040))+4.*(-2.*(-1/5040.*max(0,-207/256+mod(M.*X,1)).^7+max(0,-79/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-143/256+mod(M.*X,1)).^7+max(0,-15/256+mod(M.*X,1)).^7/5040))))+16.*(-8.*(-4.*(-2.*(-1/5040.*max(0,-247/256+mod(M.*X,1)).^7+max(0,-119/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-183/256+mod(M.*X,1)).^7+max(0,-55/256+mod(M.*X,1)).^7/5040))+4.*(-2.*(-1/5040.*max(0,-215/256+mod(M.*X,1)).^7+max(0,-87/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-151/256+mod(M.*X,1)).^7+max(0,-23/256+mod(M.*X,1)).^7/5040)))+8.*(-4.*(-2.*(-1/5040.*max(0,-231/256+mod(M.*X,1)).^7+max(0,-103/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-167/256+mod(M.*X,1)).^7+max(0,-39/256+mod(M.*X,1)).^7/5040))+4.*(-2.*(-1/5040.*max(0,-199/256+mod(M.*X,1)).^7+max(0,-71/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-135/256+mod(M.*X,1)).^7+max(0,-7/256+mod(M.*X,1)).^7/5040)))))+32.*(-16.*(-8.*(-4.*(-2.*(-1/5040.*max(0,-251/256+mod(M.*X,1)).^7+max(0,-123/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-187/256+mod(M.*X,1)).^7+max(0,-59/256+mod(M.*X,1)).^7/5040))+4.*(-2.*(-1/5040.*max(0,-219/256+mod(M.*X,1)).^7+max(0,-91/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-155/256+mod(M.*X,1)).^7+max(0,-27/256+mod(M.*X,1)).^7/5040)))+8.*(-4.*(-2.*(-1/5040.*max(0,-235/256+mod(M.*X,1)).^7+max(0,-107/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-171/256+mod(M.*X,1)).^7+max(0,-43/256+mod(M.*X,1)).^7/5040))+4.*(-2.*(-1/5040.*max(0,-203/256+mod(M.*X,1)).^7+max(0,-75/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-139/256+mod(M.*X,1)).^7+max(0,-11/256+mod(M.*X,1)).^7/5040))))+16.*(-8.*(-4.*(-2.*(-1/5040.*max(0,-243/256+mod(M.*X,1)).^7+max(0,-115/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-179/256+mod(M.*X,1)).^7+max(0,-51/256+mod(M.*X,1)).^7/5040))+4.*(-2.*(-1/5040.*max(0,-211/256+mod(M.*X,1)).^7+max(0,-83/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-147/256+mod(M.*X,1)).^7+max(0,-19/256+mod(M.*X,1)).^7/5040)))+8.*(-4.*(-2.*(-1/5040.*max(0,-227/256+mod(M.*X,1)).^7+max(0,-99/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-163/256+mod(M.*X,1)).^7+max(0,-35/256+mod(M.*X,1)).^7/5040))+4.*(-2.*(-1/5040.*max(0,-195/256+mod(M.*X,1)).^7+max(0,-67/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-131/256+mod(M.*X,1)).^7+max(0,-3/256+mod(M.*X,1)).^7/5040))))))+64.*(-32.*(-16.*(-8.*(-4.*(-2.*(-1/5040.*max(0,-253/256+mod(M.*X,1)).^7+max(0,-125/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-189/256+mod(M.*X,1)).^7+max(0,-61/256+mod(M.*X,1)).^7/5040))+4.*(-2.*(-1/5040.*max(0,-221/256+mod(M.*X,1)).^7+max(0,-93/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-157/256+mod(M.*X,1)).^7+max(0,-29/256+mod(M.*X,1)).^7/5040)))+8.*(-4.*(-2.*(-1/5040.*max(0,-237/256+mod(M.*X,1)).^7+max(0,-109/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-173/256+mod(M.*X,1)).^7+max(0,-45/256+mod(M.*X,1)).^7/5040))+4.*(-2.*(-1/5040.*max(0,-205/256+mod(M.*X,1)).^7+max(0,-77/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-141/256+mod(M.*X,1)).^7+max(0,-13/256+mod(M.*X,1)).^7/5040))))+16.*(-8.*(-4.*(-2.*(-1/5040.*max(0,-245/256+mod(M.*X,1)).^7+max(0,-117/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-181/256+mod(M.*X,1)).^7+max(0,-53/256+mod(M.*X,1)).^7/5040))+4.*(-2.*(-1/5040.*max(0,-213/256+mod(M.*X,1)).^7+max(0,-85/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-149/256+mod(M.*X,1)).^7+max(0,-21/256+mod(M.*X,1)).^7/5040)))+8.*(-4.*(-2.*(-1/5040.*max(0,-229/256+mod(M.*X,1)).^7+max(0,-101/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-165/256+mod(M.*X,1)).^7+max(0,-37/256+mod(M.*X,1)).^7/5040))+4.*(-2.*(-1/5040.*max(0,-197/256+mod(M.*X,1)).^7+max(0,-69/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-133/256+mod(M.*X,1)).^7+max(0,-5/256+mod(M.*X,1)).^7/5040)))))+32.*(-16.*(-8.*(-4.*(-2.*(-1/5040.*max(0,-249/256+mod(M.*X,1)).^7+max(0,-121/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-185/256+mod(M.*X,1)).^7+max(0,-57/256+mod(M.*X,1)).^7/5040))+4.*(-2.*(-1/5040.*max(0,-217/256+mod(M.*X,1)).^7+max(0,-89/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-153/256+mod(M.*X,1)).^7+max(0,-25/256+mod(M.*X,1)).^7/5040)))+8.*(-4.*(-2.*(-1/5040.*max(0,-233/256+mod(M.*X,1)).^7+max(0,-105/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-169/256+mod(M.*X,1)).^7+max(0,-41/256+mod(M.*X,1)).^7/5040))+4.*(-2.*(-1/5040.*max(0,-201/256+mod(M.*X,1)).^7+max(0,-73/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-137/256+mod(M.*X,1)).^7+max(0,-9/256+mod(M.*X,1)).^7/5040))))+16.*(-8.*(-4.*(-2.*(-1/5040.*max(0,-241/256+mod(M.*X,1)).^7+max(0,-113/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-177/256+mod(M.*X,1)).^7+max(0,-49/256+mod(M.*X,1)).^7/5040))+4.*(-2.*(-1/5040.*max(0,-209/256+mod(M.*X,1)).^7+max(0,-81/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-145/256+mod(M.*X,1)).^7+max(0,-17/256+mod(M.*X,1)).^7/5040)))+8.*(-4.*(-2.*(-1/5040.*max(0,-225/256+mod(M.*X,1)).^7+max(0,-97/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-161/256+mod(M.*X,1)).^7+max(0,-33/256+mod(M.*X,1)).^7/5040))+4.*(-2.*(-1/5040.*max(0,-193/256+mod(M.*X,1)).^7+max(0,-65/256+mod(M.*X,1)).^7/5040)+2.*(-1/5040.*max(0,-129/256+mod(M.*X,1)).^7+max(0,-1/256+mod(M.*X,1)).^7/5040))))))))).*(abs(X)<=1);end
|
|
function f=O_NOITALOPRETNI_8SUIBAF_O_FABIUS8_INTERPOLATION_O(X),M=256/256;f=(0.5-(-1).^floor(M.*X).*(-0.5+256.*(-128.*(-64.*(-32.*(-16.*(-8.*(-4.*(-2.*(-1/40320.*max(0,-511/512+mod(M.*X,1)).^8+max(0,-255/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-383/512+mod(M.*X,1)).^8+max(0,-127/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-447/512+mod(M.*X,1)).^8+max(0,-191/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-319/512+mod(M.*X,1)).^8+max(0,-63/512+mod(M.*X,1)).^8/40320)))+8.*(-4.*(-2.*(-1/40320.*max(0,-479/512+mod(M.*X,1)).^8+max(0,-223/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-351/512+mod(M.*X,1)).^8+max(0,-95/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-415/512+mod(M.*X,1)).^8+max(0,-159/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-287/512+mod(M.*X,1)).^8+max(0,-31/512+mod(M.*X,1)).^8/40320))))+16.*(-8.*(-4.*(-2.*(-1/40320.*max(0,-495/512+mod(M.*X,1)).^8+max(0,-239/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-367/512+mod(M.*X,1)).^8+max(0,-111/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-431/512+mod(M.*X,1)).^8+max(0,-175/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-303/512+mod(M.*X,1)).^8+max(0,-47/512+mod(M.*X,1)).^8/40320)))+8.*(-4.*(-2.*(-1/40320.*max(0,-463/512+mod(M.*X,1)).^8+max(0,-207/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-335/512+mod(M.*X,1)).^8+max(0,-79/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-399/512+mod(M.*X,1)).^8+max(0,-143/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-271/512+mod(M.*X,1)).^8+max(0,-15/512+mod(M.*X,1)).^8/40320)))))+32.*(-16.*(-8.*(-4.*(-2.*(-1/40320.*max(0,-503/512+mod(M.*X,1)).^8+max(0,-247/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-375/512+mod(M.*X,1)).^8+max(0,-119/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-439/512+mod(M.*X,1)).^8+max(0,-183/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-311/512+mod(M.*X,1)).^8+max(0,-55/512+mod(M.*X,1)).^8/40320)))+8.*(-4.*(-2.*(-1/40320.*max(0,-471/512+mod(M.*X,1)).^8+max(0,-215/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-343/512+mod(M.*X,1)).^8+max(0,-87/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-407/512+mod(M.*X,1)).^8+max(0,-151/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-279/512+mod(M.*X,1)).^8+max(0,-23/512+mod(M.*X,1)).^8/40320))))+16.*(-8.*(-4.*(-2.*(-1/40320.*max(0,-487/512+mod(M.*X,1)).^8+max(0,-231/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-359/512+mod(M.*X,1)).^8+max(0,-103/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-423/512+mod(M.*X,1)).^8+max(0,-167/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-295/512+mod(M.*X,1)).^8+max(0,-39/512+mod(M.*X,1)).^8/40320)))+8.*(-4.*(-2.*(-1/40320.*max(0,-455/512+mod(M.*X,1)).^8+max(0,-199/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-327/512+mod(M.*X,1)).^8+max(0,-71/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-391/512+mod(M.*X,1)).^8+max(0,-135/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-263/512+mod(M.*X,1)).^8+max(0,-7/512+mod(M.*X,1)).^8/40320))))))+64.*(-32.*(-16.*(-8.*(-4.*(-2.*(-1/40320.*max(0,-507/512+mod(M.*X,1)).^8+max(0,-251/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-379/512+mod(M.*X,1)).^8+max(0,-123/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-443/512+mod(M.*X,1)).^8+max(0,-187/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-315/512+mod(M.*X,1)).^8+max(0,-59/512+mod(M.*X,1)).^8/40320)))+8.*(-4.*(-2.*(-1/40320.*max(0,-475/512+mod(M.*X,1)).^8+max(0,-219/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-347/512+mod(M.*X,1)).^8+max(0,-91/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-411/512+mod(M.*X,1)).^8+max(0,-155/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-283/512+mod(M.*X,1)).^8+max(0,-27/512+mod(M.*X,1)).^8/40320))))+16.*(-8.*(-4.*(-2.*(-1/40320.*max(0,-491/512+mod(M.*X,1)).^8+max(0,-235/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-363/512+mod(M.*X,1)).^8+max(0,-107/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-427/512+mod(M.*X,1)).^8+max(0,-171/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-299/512+mod(M.*X,1)).^8+max(0,-43/512+mod(M.*X,1)).^8/40320)))+8.*(-4.*(-2.*(-1/40320.*max(0,-459/512+mod(M.*X,1)).^8+max(0,-203/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-331/512+mod(M.*X,1)).^8+max(0,-75/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-395/512+mod(M.*X,1)).^8+max(0,-139/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-267/512+mod(M.*X,1)).^8+max(0,-11/512+mod(M.*X,1)).^8/40320)))))+32.*(-16.*(-8.*(-4.*(-2.*(-1/40320.*max(0,-499/512+mod(M.*X,1)).^8+max(0,-243/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-371/512+mod(M.*X,1)).^8+max(0,-115/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-435/512+mod(M.*X,1)).^8+max(0,-179/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-307/512+mod(M.*X,1)).^8+max(0,-51/512+mod(M.*X,1)).^8/40320)))+8.*(-4.*(-2.*(-1/40320.*max(0,-467/512+mod(M.*X,1)).^8+max(0,-211/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-339/512+mod(M.*X,1)).^8+max(0,-83/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-403/512+mod(M.*X,1)).^8+max(0,-147/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-275/512+mod(M.*X,1)).^8+max(0,-19/512+mod(M.*X,1)).^8/40320))))+16.*(-8.*(-4.*(-2.*(-1/40320.*max(0,-483/512+mod(M.*X,1)).^8+max(0,-227/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-355/512+mod(M.*X,1)).^8+max(0,-99/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-419/512+mod(M.*X,1)).^8+max(0,-163/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-291/512+mod(M.*X,1)).^8+max(0,-35/512+mod(M.*X,1)).^8/40320)))+8.*(-4.*(-2.*(-1/40320.*max(0,-451/512+mod(M.*X,1)).^8+max(0,-195/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-323/512+mod(M.*X,1)).^8+max(0,-67/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-387/512+mod(M.*X,1)).^8+max(0,-131/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-259/512+mod(M.*X,1)).^8+max(0,-3/512+mod(M.*X,1)).^8/40320)))))))+128.*(-64.*(-32.*(-16.*(-8.*(-4.*(-2.*(-1/40320.*max(0,-509/512+mod(M.*X,1)).^8+max(0,-253/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-381/512+mod(M.*X,1)).^8+max(0,-125/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-445/512+mod(M.*X,1)).^8+max(0,-189/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-317/512+mod(M.*X,1)).^8+max(0,-61/512+mod(M.*X,1)).^8/40320)))+8.*(-4.*(-2.*(-1/40320.*max(0,-477/512+mod(M.*X,1)).^8+max(0,-221/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-349/512+mod(M.*X,1)).^8+max(0,-93/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-413/512+mod(M.*X,1)).^8+max(0,-157/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-285/512+mod(M.*X,1)).^8+max(0,-29/512+mod(M.*X,1)).^8/40320))))+16.*(-8.*(-4.*(-2.*(-1/40320.*max(0,-493/512+mod(M.*X,1)).^8+max(0,-237/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-365/512+mod(M.*X,1)).^8+max(0,-109/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-429/512+mod(M.*X,1)).^8+max(0,-173/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-301/512+mod(M.*X,1)).^8+max(0,-45/512+mod(M.*X,1)).^8/40320)))+8.*(-4.*(-2.*(-1/40320.*max(0,-461/512+mod(M.*X,1)).^8+max(0,-205/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-333/512+mod(M.*X,1)).^8+max(0,-77/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-397/512+mod(M.*X,1)).^8+max(0,-141/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-269/512+mod(M.*X,1)).^8+max(0,-13/512+mod(M.*X,1)).^8/40320)))))+32.*(-16.*(-8.*(-4.*(-2.*(-1/40320.*max(0,-501/512+mod(M.*X,1)).^8+max(0,-245/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-373/512+mod(M.*X,1)).^8+max(0,-117/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-437/512+mod(M.*X,1)).^8+max(0,-181/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-309/512+mod(M.*X,1)).^8+max(0,-53/512+mod(M.*X,1)).^8/40320)))+8.*(-4.*(-2.*(-1/40320.*max(0,-469/512+mod(M.*X,1)).^8+max(0,-213/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-341/512+mod(M.*X,1)).^8+max(0,-85/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-405/512+mod(M.*X,1)).^8+max(0,-149/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-277/512+mod(M.*X,1)).^8+max(0,-21/512+mod(M.*X,1)).^8/40320))))+16.*(-8.*(-4.*(-2.*(-1/40320.*max(0,-485/512+mod(M.*X,1)).^8+max(0,-229/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-357/512+mod(M.*X,1)).^8+max(0,-101/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-421/512+mod(M.*X,1)).^8+max(0,-165/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-293/512+mod(M.*X,1)).^8+max(0,-37/512+mod(M.*X,1)).^8/40320)))+8.*(-4.*(-2.*(-1/40320.*max(0,-453/512+mod(M.*X,1)).^8+max(0,-197/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-325/512+mod(M.*X,1)).^8+max(0,-69/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-389/512+mod(M.*X,1)).^8+max(0,-133/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-261/512+mod(M.*X,1)).^8+max(0,-5/512+mod(M.*X,1)).^8/40320))))))+64.*(-32.*(-16.*(-8.*(-4.*(-2.*(-1/40320.*max(0,-505/512+mod(M.*X,1)).^8+max(0,-249/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-377/512+mod(M.*X,1)).^8+max(0,-121/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-441/512+mod(M.*X,1)).^8+max(0,-185/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-313/512+mod(M.*X,1)).^8+max(0,-57/512+mod(M.*X,1)).^8/40320)))+8.*(-4.*(-2.*(-1/40320.*max(0,-473/512+mod(M.*X,1)).^8+max(0,-217/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-345/512+mod(M.*X,1)).^8+max(0,-89/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-409/512+mod(M.*X,1)).^8+max(0,-153/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-281/512+mod(M.*X,1)).^8+max(0,-25/512+mod(M.*X,1)).^8/40320))))+16.*(-8.*(-4.*(-2.*(-1/40320.*max(0,-489/512+mod(M.*X,1)).^8+max(0,-233/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-361/512+mod(M.*X,1)).^8+max(0,-105/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-425/512+mod(M.*X,1)).^8+max(0,-169/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-297/512+mod(M.*X,1)).^8+max(0,-41/512+mod(M.*X,1)).^8/40320)))+8.*(-4.*(-2.*(-1/40320.*max(0,-457/512+mod(M.*X,1)).^8+max(0,-201/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-329/512+mod(M.*X,1)).^8+max(0,-73/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-393/512+mod(M.*X,1)).^8+max(0,-137/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-265/512+mod(M.*X,1)).^8+max(0,-9/512+mod(M.*X,1)).^8/40320)))))+32.*(-16.*(-8.*(-4.*(-2.*(-1/40320.*max(0,-497/512+mod(M.*X,1)).^8+max(0,-241/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-369/512+mod(M.*X,1)).^8+max(0,-113/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-433/512+mod(M.*X,1)).^8+max(0,-177/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-305/512+mod(M.*X,1)).^8+max(0,-49/512+mod(M.*X,1)).^8/40320)))+8.*(-4.*(-2.*(-1/40320.*max(0,-465/512+mod(M.*X,1)).^8+max(0,-209/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-337/512+mod(M.*X,1)).^8+max(0,-81/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-401/512+mod(M.*X,1)).^8+max(0,-145/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-273/512+mod(M.*X,1)).^8+max(0,-17/512+mod(M.*X,1)).^8/40320))))+16.*(-8.*(-4.*(-2.*(-1/40320.*max(0,-481/512+mod(M.*X,1)).^8+max(0,-225/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-353/512+mod(M.*X,1)).^8+max(0,-97/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-417/512+mod(M.*X,1)).^8+max(0,-161/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-289/512+mod(M.*X,1)).^8+max(0,-33/512+mod(M.*X,1)).^8/40320)))+8.*(-4.*(-2.*(-1/40320.*max(0,-449/512+mod(M.*X,1)).^8+max(0,-193/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-321/512+mod(M.*X,1)).^8+max(0,-65/512+mod(M.*X,1)).^8/40320))+4.*(-2.*(-1/40320.*max(0,-385/512+mod(M.*X,1)).^8+max(0,-129/512+mod(M.*X,1)).^8/40320)+2.*(-1/40320.*max(0,-257/512+mod(M.*X,1)).^8+max(0,-1/512+mod(M.*X,1)).^8/40320)))))))))).*(abs(X)<=1);end
|
|
|