updated README

Author: matteoicardi

README.md

@@ -1 +1,19 @@
-# fractionalPorousMedia
+# fractionalPorousMedia
+
+Python code for solving fractional porous media equation using finite difference method and product-integration rule.
+
+## Connected paper
+Berardi, Garrappa, Icardi, Nuca
+Basset-type fractional models for heterogeneous porous media: model derivation, analysis and numerical simulations
+
+## Instructions
+The main code to reproduce the results of the paper is `main_BassetEqn.py`. The code is written in Python 3.7 and requires the following packages:
+- numpy
+- matplotlib
+
+Other extra files in python and MATLAB are only for testing and comparison purposes.
+
+
+## Authors
+- Matteo Icardi
+- Roberto Nuca

matlab/eigenvalue.m

@@ -1,5 +1,5 @@
 %% Setup
-b = 1;
+b = -1;
 x=chebfun('x',[0,1]);
 y=chebfun('y',[-b.^2/4+1e-5,100]);
 
@@ -27,7 +27,7 @@
 
 %% Eigenvalues with chebfun
 
-L=chebop(@(u) -diff(u,2) + b*diff(u),[0,1]);
+L=chebop(@(u) diff(u,2) + b*diff(u),[0,1]);
 L.lbc='dirichlet';
 L.rbc='neumann';