plot_PM.m

function []=plot_PM(opt_data)

[m,n]=size(opt_data.parameter_sets);

%hold on

switch n
    case 2
        for k=1:(m-1)
            vector_comp = opt_data.parameter_sets(k+1,:)-opt_data.parameter_sets(k,:);
            quiver(opt_data.parameter_sets(k,1), opt_data.parameter_sets(k,2), vector_comp(1), vector_comp(2),0,'LineWidth',2)
            hold on
        end
        if isfield(opt_data,'Delta_mu')
            for k=2:m
                circle_center_x = opt_data.parameter_sets(k,1);
                circle_center_y = opt_data.parameter_sets(k,2);
                [x,y]=circle_points(circle_center_x,circle_center_y,opt_data.Delta_mu(k));
                plot(x,y)
                
            end
        end
        
    case 3
        
        for k=1:(m-1)
            vector_comp = opt_data.parameter_sets(k+1,:)-opt_data.parameter_sets(k,:);
            quiver3(opt_data.parameter_sets(k,1), opt_data.parameter_sets(k,2), opt_data.parameter_sets(k,3), vector_comp(1), vector_comp(2), vector_comp(3),0,'LineWidth',2)            
            hold on
        end  
        
    otherwise
        disp('only 2D or 3D Plots possible, size of opt_data.PM is too large or 1')
end
        
hold off
axis equal
end %plot_PM

function [x,y] = circle_points(x_m,y_m,r)

steps=100;
t=0:steps;

x= x_m + r*sin(t*2*pi/steps);
y= y_m + r*cos(t*2*pi/steps);
end