How to evaluate the following equation with Ei function?

1 次查看（过去 30 天）
Daniyar
Daniyar 2013-5-12
Hey, guys
I am having trouble in evaluation of the following equation: rd=r/rw, td=k*t/(phi*mu*c*rw^2), pd=((2*pi*k*h)/(q*B*mu)*(IP-p(r,t))
pd(rd,td)=-1/2*Ei*(-rd^2/4*td)> the following equation needs to be evaluated
t=3000000 nr=nx=161 (to be introduced instead of r) q=0.02 k=5e-13 mu=0.00035 c=1.5e-8 phi=0.25 rw=0.1 h=4.5 B=1.1
My code is below: st=size(t) nr=nx;
if true
% for i=1:1:st(1)
rd=r(1:nr)/wb
td=0.0036*k*t(i)/(phi*mu*c*wb^2)
pd=((pi*h*k)/(1.842*q*1.1*mu))*(IP-p(i,1:nr))
Ei_from_mfun=mfun('Ei',1,-rd^2/(4*td))
pd=1/2*Ei_from_mfun
plot(time,pd,'bo')
end
end
The following command appears: if true % ??? Index exceeds matrix dimensions. Error in ==> radialflow1 at 103 pd=((pi*h*k)/(1.842*q*1.1*mu))*(IP-p(i,1:nr))
end
Can you please help to fix the problem for every time, t and pressure,p?
Thanks

请先登录，再回答此问题。

回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Elementary Math 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by