0
|
10968 Views
|
5 Replies
|
5 Total Likes
View groups...
Share
GROUPS:

# Differential Algebraic Equation

Posted 10 years ago
 There is an error in solving the attached differential algebraic equation. Someone kindly suggest me how to remove the error. Attachments:
5 Replies
Sort By:
Posted 10 years ago
 Hi,what is f[zeta]? This function is not defined.Suvadip.
Posted 10 years ago
 Sir, There are two coupled equations of f[zeta] and eta[zeta]. From these the equation of f[zeta] is second order differential equation and the equation of eta[zeta] is algebraic equation. I want to plot the graph of eta[zeta] with zeta. So f[zeta] is a function which is to be used in eta[zeta].
Posted 10 years ago
 Hi,I have looked at your attachment. One of the equations you wish to solve has the form {-(1/f[z]^3) + 2.03742*10^7/((1/3 - 0.166597 (1 + 8.14968/f[z]^2)^3.)^2 (1 + 8.14968/f[z]^2) f[z]^3) + (f^\[Prime]\[Prime])[z]} == 0 The left hand side is a list and the right hand side is a number. Therefore, it is impossible for a solution to exist. There is a coding mistake somewhere.
Posted 10 years ago
 Perhaps I have misunderstood some part of your attachement.If I replace your NDSolve and the plots that follow it with this f = First[f/.NDSolve[{f''[\[Xi]]-1/f[\[Xi]]^3+1/eo f[\[Xi]] phir[\[Xi]] 10^7==0, f[0]==1, f'[0]==0}, f,{\[Xi],0,1}]]; Plot[f[\[Xi]], {\[Xi], 0, 1}, AxesOrigin->{0,0}, PlotStyle->{Thickness[0.006]}, AxesLabel->{f}, PlotRange->All] Plot[eta[\[Xi]]-Abs[(\[Xi] Exp[b1[\[Xi]]] a1[\[Xi]] Sqrt[pi/2])/2*Exp[b1[\[Xi]]^2/8] (Erf[Sqrt[2]*\[Xi]+b1[\[Xi]]/ (2 Sqrt[2])]-Erf[Sqrt[2]*(\[Xi]*(1-\[Beta][\[Xi]])-\[Beta][\[Xi]]*tdash)+b1[\[Xi]]/(2 Sqrt[2])])], {\[Xi], 0, 1}, AxesOrigin->{0,0}, PlotStyle->{Thickness[0.006]}, AxesLabel->{\[Xi]}, PlotRange->All] then it seems to find the desired function f and it then uses that f to evaluate your next expression.I am concerned that in your attachment you seem to want that next expression to equal zero, but from the second plot that does not seem possible.If I have made some mistake and you can explain what that was then I will see if I can correct it.
Posted 10 years ago
 Thanks Sir for your reply. Sir, actually there are two coupled equations of f[zeta] and eta[zeta]. From these the equation of f[zeta] is second order differential equation and the equation of eta[zeta] is algebraic equation. I want to plot the graph of eta[zeta] with zeta and in the equation of eta[zeta], f[zeta] exists. So, my purpose is to solve two coupled equations in f[zeta] and eta[zeta]. Although i can solve these two equations individually, but if i try to solve them collectively, error occurs. Kindly try to help me.