假定有U(r,z)在z方向上分3层,已知U1(r,z)、U2(r,z)、U3(r,z)
U1(r,z)含待定参数a1、b1
U2(r,z)含待定参数a2、b2
U3(r,z)含待定参数a3、b3
根据连续性和边界条件求待定参数
本来想用solve('U1(r,z1)=0','U1(r,z2)=U2(r,z2)','U3(r,z3)=0')解可是出现错误
其中z1、z2、z3为已知,假定为1、4、8
U1=inline(Q1,'r','z') %%Q1为字符型表达式,U2、U3类似
请问如何编程,关键是用什么函数?(U(r,z)比较复杂,根贝塞尔函数有关) |