function?[Uxt]=rechuandao（cfg1g2xspantspanngrid）?%?
%?ngrid=[nm]
n=ngrid（1）;?
m=ngrid（2）;?
h=range（xspan）/（m-1）;?
x=linspace（xspan（1）xspan（2）m）;?
k=range（tspan）/（n-1）;
t=linspace（tspan（1）tspan（2）n）;??
r=c*k/h^2