MathGroup Archive 2013

[Date Index] [Thread Index] [Author Index]

Search the Archive

ParametricNDSolve for a contact problem

  • To: mathgroup at
  • Subject: [mg131832] ParametricNDSolve for a contact problem
  • From: Simon Pearce <Simon.Pearce at>
  • Date: Sat, 12 Oct 2013 02:38:26 -0400 (EDT)
  • Delivered-to:
  • Delivered-to:
  • Delivered-to:
  • Delivered-to:
  • References: <>

Hi Mathgroup,

I've just installed v9.0.1.0, and I was hopeful that ParametricNDSolve would help simplify my code for solving a contact problem. I have an ODE which I don't know the starting point at which my initial conditions are to be set, which I need to vary.

The following toy model works, and shows that the range of integration may be set as a parameter:

sol1 = ParametricNDSolve[{y'[t] == y[t], y[0] == 0}, y, {t, b, 1}, {b}];
y[0] /. sol1

Whereas when the initial condition is to be set at the parameterised starting point it does not work:

sol2 = ParametricNDSolve[{y'[t] == y[t], y[b] == 0}, y, {t, b, 1}, {b}];
y[0] /. sol2

Any idea on how to get ParametricNDSolve to do this for me? My current method is to use NDSolve`ProcessEquations and NDSolve`Reinitialize inside a function which I send to FindRoot.

My actual problem is a highly nonlinear fifth order system, with one of the initial conditions at the unknown contact point also unknown.

Simon Pearce
University of Nottingham, UK
This message and any attachment are intended solely for the addressee and may contain confidential information. If you have received this message in error, please send it back to me, and immediately delete it.   Please do not use, copy or disclose the information contained in this message or in any attachment.  Any views or opinions expressed by the author of this email do not necessarily reflect the views of the University of Nottingham.

This message has been checked for viruses but the contents of an attachment
may still contain software viruses which could damage your computer system, you are advised to perform your own checks. Email communications with the University of Nottingham may be monitored as permitted by UK legislation.

  • Prev by Date: Re: Decision whether point lies within closed contour
  • Next by Date: Re: Image color ratio
  • Previous by thread: Re: bug in Print
  • Next by thread: Re: bug in Print