Advanced Search

Question: substituting a procedure in place of a function

Question:substituting a procedure in place of a function

rhaynes 5
procedure
April 20 2013
1

Hi folks, I am solving a boundary value problem numerically using dsolve. I would like to pass the resulting solution (which I gather is really a procedure) back into another differential equation as a coefficient.

 

M:=x(xi)^2+1:

ode1:=diff(M*diff(x(xi),xi),xi)-0*x(xi)=0:

ics:=x(0)=0,x(1)=1:

sol1:=dsolve({ode1,ics},numeric,output=listprocedure):

y_sol1:=sol1[2]:

y_sol1:=rhs(y_sol1):

% I could then plot the solution using 

plot(y_sol1(xi),xi=0..1);

% I would then like to do something like

M:=x(xi)^2+1:

M:=subs(x(xi)=y_sol1,M);    %% but this doesn't work

ode2:=diff(M*diff(x(xi),xi),xi)-0*x(xi)=0:

etc...  

 

Any suggestions?

 

Thanks,

RH
Please Wait...