微分方程数值解向前差分.docVIP

宋玉霞 座机电话号码16 统数 信科0802 实验三__1:向前差分 一、实验题目： 二、实验要求：用向前差分格式求解一维抛物型偏微分方程 三、程序运行结果： 只取前面与后面的结果：其数值解为 u Columns 1 through 8 1.0000 1.0025 1.0050 1.0075 1.0101 1.0126 1.0151 1.0177 1.1052 1.1079 1.1107 1.1135 1.1163 1.1190 1.1218 1.1244 1.2214 1.2245 1.2275 1.2306 1.2336 1.2365 1.2393 1.2417 1.3499 1.3532 1.3566 1.3599 1.3629 1.3655 1.3672 1.3680 1.4918 1.4955 1.4991 1.5020 1.5036 1.5033 1.5005 1.4950 1.6487 1.6526 1.6551 1.6540 1.6476 1.6352 1.6173 1.5949 1.8221 1.8241 1.8163 1.7924 1.7538 1.7056 1.6525 1.5977 2.0138 1.9936 1.9164 1.8007 1.6811 1.5695 1.4693 1.3806 2.2255 1.9820 1.5802 1.3199 1.1401 1.0091 0.9093 0.8306 2.4596 0 0 0 0 0 0 0 2.7183 2.7251 2.7319 2.7387 2.7456 2.7525 2.7594 2.7663 Columns 9 through 16 1.0202 1.0228 1.0253 1.0279 1.0305 1.0330 1.0356 1.0382 1.1270 1.1294 1.1316 1.1335 1.1351 1.1364 1.1374 1.1380 1.2437 1.2452 1.2461 1.2462 1.2457 1.2444 1.2424 1.2397 1.3675 1.3658 1.3627 1.3582 1.3526 1.3457 1.3379 1.3292 1.4870 1.4766 1.4641 1.4499 1.4343 1.4177 1.4002 1.3821 1.5693 1.5413 1.5119 1.4818 1.4513 1.4211 1.3912 1.3619 1.5433 1.4905 1.4400 1.3920 1.3467 1.3039 1.2636 1.2257 1.3022 1.2327 1.1709 1.1155 1.0658 1.0207 0.9798 0.9424 0.7666 0.7135 0.6684 0.6297 0.5960 0.5662 0.5398 0.5161 0 0 0 0 0 0 0 0 2.7732 2.7801 2.7871 2.7941 2.8011 2.8081 2.8151 2.8222 Columns 385 through 392 2.6117 2.6182 2.6248 2.6314 2.6379 2.6445 2.6512 2.6578 2.8863 2.8935 2.9008 2.9080 2.9153 2.9226 2.9299 2.9373 3.1898 3.1978 3.2058 3.2138 3.2219 3.2299 3.2380 3.2461 3.5253 3.5341 3.54