function err = sol6(x,param) % To solve ODE y(t)= (D*A+beta1*J1+beta2*J2)*y(t)
z=param.z;A=param.A;J1=param.J1;J2=param.J2;l10=param.l10;l35=param.l35;l50=param.l50;lpel=param.lpel;lbrk2=param.lbrk2;T=param.T;
C11=param.C11;C12=param.C12;C13=param.C13;SD11=param.SD11;SD12=param.SD12;SD13=param.SD13;
y0=x(2)*100*[(exp(-((z(1:lpel-1)-17.5).^2)/x(5)));(exp(-((z(lpel:lbrk2-1)-17.5).^2)/x(6)));exp(-((z(lbrk2)-17.5).^2)/x(6))];
y0=[y0;.015*(x(3)/x(4))+0*y0];
[t,y] = ode45(@(t,y)(x(1)*A+x(3)*J1+x(4)*J2)*y,[0,323],y0);
kdp=zeros(size(T)); % Indices corresponding to the data points
for j=1:size(kdp,1)
    kdp(j)=find(abs(t-T(j))==min(abs(t-T(j))),1);
end
err = (norm((y(kdp,l10)-C11)./SD11))^2+(norm((y(kdp,l35)-C12)./SD12))^2+(norm((y(kdp,l50)-C13)./SD13))^2; 
end