您好 我想请问一个一维热传导的偏微分的方程差分格式 能否帮忙？ Grank-Nicholson方法源程序：function[u，x，t]=Grank_Nicholson(A，xf，T，it0，bx0，bxf，M，N)解方程 A u_xx=u_t，0，0初值：u(x，0)=it0(x)边界条件：u(0，t)=bx0(t)，u(xf，t)=bxf(t)M：x 轴的等分段数N：t 轴的等分段数dx=xf/M；x=[0：M]*dx；dt=T/N；t=[0：N]'*dt；for i=1：M+1u(i，1)=it0(x(i))；endfor n=1：N+1u([1 M+1]，n)=[bx0(t(n))；bxf(t(n))]；endr=A*dt/dx/dx；r1=2*(1+r)；r2=2*(1-r)；for i=1：M-1P(i，i)=r1；(9.2.17)Q(i，i)=r2；if i>；1P(i-1，i)=-r；P(i，i-1)=-r；(9.2.17)等式左边矩阵Q(i-1，i)=r；Q(i，i-1)=r；(9.2.17)等式右边矩阵endendfor k=2：N+1b=Q*u(2：M，k-1)+[r*(u(1，k)+u(1，k-1))；zeros(M-2，1)]；u(2：M，k)=linsolve(P，b)；(9.2.17)endu=u'；例2.1 Grank-Nicholson方法求解一维抛物性方程应用实例。求满足以下条件的热传导数值解：自变量取值：边界：解：在MATLAB中编写脚本文件：A=0.5；方程系数it0=inline('sin(pi*x)'，'x')；初始条件bx0=inline('0')；bxf=inline('0')；边界条件xf=2；M=25；T=0.1；N=100；[u1，x，t]=Grank_Nicholson(A，xf，T，it0，bx0，bxf，M，N)；mesh(u1)xlabel('x')ylabel('t')zlabel('U')
用matlab求解抛物型方程，急啊！！用最简隐格式（向后差分格式）求解抛物型方程 用matlab求解抛物型方程，急啊！用最简隐格式（向后差分格式）求解抛物型方程 要用matlab求解，但是不能用里面的求微分方程的工具来求解，就是自己编程序，要有图示的啊，。
急求！！！ 大学数学，用matlab解决问题，题目是一维抛物型偏微分方程差分解法
请求各位大虾，帮忙指点一下，如何用MATLAB解一维抛物型方程？谢谢！ 可以用命令pdetool进入pde工具箱来用数值方法求解
旋转抛物面方程怎么写？
热传导方程为何是抛物型方程 一维热传导方程是抛物型的，因为a12^2-a11*a22=0。书上有
