function [mx,my,dx,dy,ax,ay,bx,by,iprob,fname] = read_dat();

fid = fopen('fort.matlab','r');

[data,info] =fscanf(fid,'%d %d',2);
mx = data(1);
my = data(2);

[data,info] =fscanf(fid,'%f %f',2);
dx = data(1);
dy = data(2);

[data,info] = fscanf(fid,'%f %f',4);
ax = data(1);
ay = data(2);
bx = data(3);
by = data(4);

[data,info] = fscanf(fid,'%d',1);
iprob = data(1);

[data,info] = fscanf(fid,'%s',1);
fname = data;
