matlab file, forward euler method