saulyev.zip_matlab_saulyev_saulyev差分_solidkbd_抛物型方程

function [] = saulyev_main() clear; data; fid = fopen('saulyev_result.txt','a'); fprintf(fid,"1/tao, 1/h, tol\n"); %% for i=1:4 tao = 1/3000; h = 1/(10*2^(i-1)); fprintf(fid,"%d, %d, ",1/tao,1/h); % 精确解 n = (xmax-xmin)/h; u = zeros(1,n+1); for j=1:n+1 u(j) = f_u(xmin+(j-1)*h,tmax); end % 求解 [x] = saulyev_difference(tao,h); fprintf(fid,"%d\n",max(abs(x-u)./abs(x))); end %% for i=3:6 h = 1/200; tao = 1/(2^i); fprintf(fid,"%d, %d, ",1/tao,1/h); % 精确解 n = (xmax-xmin)/h; u = zeros(1,n+1); for j=1:n+1 u(j) = f_u(xmin+(j-1)*h,tmax); end % 求解 [x] = saulyev_difference(tao,h); fprintf(fid,"%d\n",max(abs(x-u)./abs(x))); end %% fclose(fid); clear; end