您好，可以尝试咨询智能问答助手或者MWORKS.Copilot。

参考这个例题的代码，进行修改即可。
function [u,x,t] = ECD_Wave(A,xf,T,it0,i1t0,bx0,bxf,M,N)
%解方程a u_xx = u_tt for 0<=x<=xf, 0<=t<=T
% 初始条件: u(x,0) = it0(x), u_t(x,0) = i1t0(x)
% 边界条件: u(0,t)= bx0(t), u(xf,t) = bxf(t)
% M ：沿x轴的等分段数
% N ：沿y轴的等分段数
dx = xf/M; x = [0:M]'dx;
dt = T/N; t = [0:N]dt;
for i = 1:M + 1
u(i,1) = it0(x(i));
end
for k = 1:N + 1
u([1 M + 1],k) = [bx0(t(k)); bxf(t(k))];
end
r = A(dt/dx)^ 2; r1 = r/2; r2 = 2(1 - r);
u(2:M,2) = r1u(1:M - 1,1) + (1 - r)u(2:M,1) + r1u(3:M + 1,1) + dti1t0(x(2:M)); %(11.3.4)
for k = 3:N + 1
u(2:M,k) = ru(1:M - 1,k - 1) + r2u(2:M,k-1) + r*u(3:M + 1,k - 1)- u(2:M,k - 2); %(11.3.3)
end
u=u';
end
A= 1;
it0 = inline('x-x^2','x'); i1t0 = inline('0');
bx0t = inline('0'); bxft = inline('0');
xf =1; M = 100; T = 2; N = 500;
[u,x,t] = ECD_Wave(A,xf,T,it0,i1t0,bx0t,bxft,M,N);
mesh(u)
xlabel('x')
ylabel('t')
zlabel('U')