Update c_hpc.m

HCS_Toolbox_ver01/c_hpc.m

@@ -34,7 +34,7 @@ hn=hnorm(x,Gd,P);
 if hn^(-mu)>-mu*rho*h*n-tol
 F0=A+B*(K0+K)+rho*Gd;
 s=log(1+mu*rho*h*n*hn^mu)/(rho*mu);
-Q=expm(rho*Gd*log(hn))*expm(-rho*Gd*s)*expm(F0*s)*expm(-rho*Gd*log(hn));
+Q=expm(Gd*log(hn))*expm(-rho*Gd*s)*expm(F0*s)*expm(-Gd*log(hn));
 else Q=zeros(size(A));
 end;