Group Abstract Group Abstract

Message Boards Message Boards

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:
POSTED BY: Shivi Sv
5 Replies
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 BY: Bill Simpson
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.

POSTED BY: Shivi Sv

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 BY: David Vasholz

Hi,

what is f[zeta]? This function is not defined.

Suvadip.

POSTED BY: suvadip mandal
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 BY: Shivi Sv
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
Attachments
Remove
or Discard