diff --git a/ImpedanceFM.mod b/ImpedanceFM.mod
index 7d307c6..80842ab 100644
--- a/ImpedanceFM.mod
+++ b/ImpedanceFM.mod
@@ -58,12 +58,11 @@ ASSIGNED {
 
 
 VERBATIM
-#include <malloc.h>
+#include <stdlib.h>
 
 // calculate impedance for the whole range of frequencies
 
-void calc_impedances(Zr,Zi,fmax,df,rext,rmax,dr,R,sigma1,sigma2,lambda,epsilon,sigmaR)
-	double Zr[],Zi[],fmax,df,rext,rmax,dr,R,sigma1,sigma2,lambda,epsilon,sigmaR;
+void calc_impedances(double Zr[], double Zi[], double fmax, double df, double rext, double rmax, double dr, double R, double sigma1, double sigma2, double lambda,double epsilon, double sigmaR)
 	// vector Z is impedance, Z[0]=real part, Z[1]=imaginary part
    {
         double epsR,sigR,w,w2,sig,eps,den,ReZ,ImZ,r,f;
diff --git a/README b/README
index 9d117ee..78180a0 100644
--- a/README
+++ b/README
@@ -40,3 +40,9 @@ please cite the above paper.
 Alain Destexhe, 
 CNRS, Gif sur Yvette, France
 http://cns.iaf.cnrs-gif.fr
+
+Changelog
+---------
+2022-05: Updated MOD files to contain valid C++ and be compatible
+         with the upcoming versions 8.2 and 9.0 of NEURON.
+