% multA.m

function y=multA(x)
    n=length(x);
    c=[-1.0./[n-1:-1:1].^2,5,-1.0./[1:n-1].^2];
    for i=n:-1:1
        y(i,1)=c(n-i+1:2*n-i)*x;
    end
