%function mcdist

load mcsim001

L=length(J);
mF=mF(J);
mD=mD(J,:)-ones(L,1)*[10 5 40 25];
mX=mX(J,:);
%inter=60;
inter=300;
%inter=120;
dynsim=10;
dynsim2=100;

t=0:inter:(L-1)*inter;
t1=zeros(2*L-1,1);
d1=zeros(2*L-1,4);
f1=zeros(2*L-1,1);
x1=zeros(2*L-1,5);
t1(2:2:end-1)=t(2:end);
t1(3:2:end)=t(2:end);
d1(1:2:end,:)=mD;
d1(2:2:end-1,:)=d1(1:2:end-2,:);
f1(1:2:end)=mF;
f1(2:2:end-1)=mF(1:end-1);
% f1(1:end)=mF(4);
% d1=mD(4*ones(2*L-1,1),:);
DistData=[t1 d1];

%mX: P2 F100 T201 F200 F3
x1(1:2:end,:)=mX(:,1:5);
x1(2:2:end-1,:)=x1(1:2:end-2,:);
P2Set=[t1 x1(:,1)-56.4245];