function VF(p,min,max)

    x = min:0.1:max ;
    y = F_Mp(x,p) ;
    plot(x,y) ;
    axis([min max 0 1]) ;

end


function erg = F_Mp(x,p)            
	erg = [] ;                  
	for i=1:length(x)
		t = x(i) ;
		if( (t==0)||(t<0) )            % Fall t <= 0
			y1 = 1/quadgk( @(r) r.*exp(-r.^p/p) ,0,Inf) * quadgk( @(r) 1/(quadgk( @(phi) ( abs(sin(phi)).^p + abs(cos(phi)).^p ).^(1/p).^(-2) ,0,2*pi)) * quadgk( @(phi) ( abs(sin(phi)).^p + abs(cos(phi)).^p ).^(1/p).^(-2) , pi+atan(-t./(r.^p-(-t).^p).^(1/p)) , 3*pi/2-atan(-t./(r.^p-(-t).^p).^(1/p)))*r.*exp(-r.^p/p) ,-2^(1/p)*t,Inf) ;
		else                            % Fall t > 0
			y1 = 1/quadgk( @(r) r.*exp(-r.^p/p) ,0,Inf) * ( quadgk( @(r) 1/(quadgk( @(phi) ( abs(sin(phi)).^p + abs(cos(phi)).^p ).^(1/p).^(-2) ,0,2*pi)) * quadgk( @(phi) ( abs(sin(phi)).^p + abs(cos(phi)).^p ).^(1/p).^(-2) ,pi-atan(t/(r^p-t^p)^(1/p)),3/2*pi+atan(t/(r^p-t^p)^(1/p)))*r.*exp(-r.^p/p) ,2^(1/p)*t,Inf) ...
                + quadgk( @(r) ( 1/(quadgk( @(phi) ( abs(sin(phi)).^p + abs(cos(phi)).^p ).^(1/p).^(-2) ,0,2*pi)) * quadgk( @(phi) ( abs(sin(phi)).^p + abs(cos(phi)).^p ).^(1/p).^(-2) ,pi-atan(t/(r^p-t^p)^(1/p)),3/2*pi+atan(t/(r^p-t^p)^(1/p))) + 1/(quadgk( @(phi) ( abs(sin(phi)).^p + abs(cos(phi)).^p ).^(1/p).^(-2) ,0,2*pi)) * quadgk( @(phi) ( abs(sin(phi)).^p + abs(cos(phi)).^p ).^(1/p).^(-2) ,pi/2-atan(t/(r^p-t^p)^(1/p)),atan(t/(r^p-t^p)^(1/p))) )*r.*exp(-r.^p/p) ,t,2^(1/p)*t) ...
                + quadgk( @(r) r.*exp(-r.^p/p) ,0,t) ) ;
		end
		erg = [ erg , y1 ] ;
	end
end

% function erg = DAF(p,alpha,beta)		% Durchschnittsanteilsfunktion für p > 1
% 	erg = 1/(quadgk( @(phi) ( abs(sin(phi)).^p + abs(cos(phi)).^p ).^(1/p).^(-2) ,0,2*pi)) * quadgk( @(phi) ( abs(sin(phi)).^p + abs(cos(phi)).^p ).^(1/p).^(-2) ,alpha,beta) ;
% end
% 
% function erg = N(p,phi)					% Normierung für p-verallgemeinerte Sinus/Kosinus
% 	erg = ( abs(sin(phi)).^p + abs(cos(phi)).^p ).^(1/p) ;
% end
% 
% function erg = PI_p(p)					% p-verallgemeinerte Kreiszahl
%     erg = 1/2*quadgk( @(phi)N(p,phi).^(-2) ,0,2*pi) ;
% end
% 
% function erg = g(x,p)					% Dichtegenerierende Funktion
% 	erg = exp(-x/p) ;
% end
%