mainTPflow.m

global zer n A tau

zerobes

for n=1:30
    A(n)=-quadl(@un,0,1)./quadl(@ud,0,1);
end

A

x=-1:0.01:1;

tau=0.5;

y=feval(@trans,x);

plot(x,y,'b')

hold on

tau=0.2;

y=feval(@trans,x);

plot(x,y,'r')

tau=0.1;

y=feval(@trans,x);

plot(x,y,'g')

tau=0.02;

y=feval(@trans,x);

plot(x,y,'m')

axis([-1 1 0 2])