home top contents previous up next

        
        q=-4.0:0.1:4.0;
        f=[-4.0:0.4:4.0];
        %global erfScale
        erfScale=0.5;
        for i=1:size(q,2)
           f(i)=erfh(q(i),erfScale);
        end
        f
        plot(q,f,'b');
        hold on
        
        erfScale=-1.0;
        for i=1:size(q,2)
           f(i)=0.95*erfh(q(i),erfScale)+0.025;
        end
        plot(q,f,'m');
                

Copyright (C) 2009 Konstantin Kirillov