This simple function will use a system call to run Xfoil. It then loads the results into a data structure which you can then use in your Matlab script.
For example to deflect the trailing edge of a NACA0009 by 5deg at 65% chord and run multiple incidence angles...
[pol foil] = xfoil('NACA0009',[-5:15],1e6,0.2,'oper iter 150','gdes flap 0.65 0 5 exec')
The coefficents will be returned in the pol structure while the foil data and pressure distribution will be returned in foil structure.
XFOIL, first developed by Mark Drela at MIT, is a program for the design and analysis of subsonic isolated airfoils and is released under the GNU General Public License.
You can find Xfoil here: