Write a Matlab function MyIVP that solves an initial-value problem (IVP) for a system of ordinary
differential equations (ODEs) of the form
x ?(t) = f (t, x(t)),
where f : R × Rn ? Rn is an arbitrary function with one one-dimensional input (for time t) and one n-dimensional input x, and n-dimensional output. The function should implement a Runge-Kutta formula (for example, the rk34 formula or the Dormand & Prince formula).
The initial value x0 is provided by the user of MyIVP. The first line of MyIVP (saved in a file MyIVP.m) should look like this
function [xend,t,xt]=MyIVP(f,x0,tspan,N) Inputs
• f: function defining the right-hand side of the ODE. f should accept two arguments: t (a number) and x (an n-dimensional vector). The function f should return an n-dimensional vector y (the time derivative). Typical calling sequence: y=f(t,x), returning the value of f at time t in position x.
• x0: initial value where integration starts from (n-dimensional vector).
• tspan: Starting time and end time for integration. Integration has to run from time t =tspan(1)
to time t =tspan(2).
• N: number of steps for integration. The integration stepsize h=(tspan(2)-tspan(1))/N should
be small.
Outputs
• xend: result of integration at t =tspan(2).
• t: vector of times at which intermediate values have been computed (this should have N + 1
entries).
• xt: intermediate values (n × (N + 1)-array). xt(:,k) should be the solution at t(k).
You can check the built-in variable nargout inside your function to see if the user wants to get three outputs or only the end value xend. If nargout==1 you don’t need to store the intermediate values.

Posted Date: 3/12/2012 10:19:33 AM | Location : United States

Ask an Expert

Related Discussions:- Initial value problems, math, Assignment Help, Ask Question on Initial value problems, math, Get Answer, Expert's Help, Initial value problems, math Discussions

Write discussion on Initial value problems, math
Your posts are moderated

Write your message here..

Related Questions

Math problem , -255=-14t-t -255=-14t-t

Prove that ar= 3/7 ac of parallelogram , ABCD is a parallelogram in the giv... ABCD is a parallelogram in the given figure, AB is divided at P and CD and Q so that AP:PB=3:2 and CQ:QD=4:1. If PQ meets AC at R, prove that AR= 3/7 AC. Ans:    ΔAPR ∼ Δ

Solid mensuration, Formulas Formulas

What is the probability that they will both come to a party, Q.  Suppose th... Q.  Suppose the probability of David coming to a party is 75% and the probability of Jason coming to a party is 85%. What is the probability that they will both come to a party, a

Math makes sense pg 261 #3 c., A seahorse layes about 200 eggs.How would yo... A seahorse layes about 200 eggs.How would you include this data on your pictograph.would you need to change anything.Explain the change.show your work.

Find how much more space than the toy it will cover, A Solid toy in the for... A Solid toy in the form of a hemisphere surmounted by the right circular cone of height  2cm  and  diameter  of  the  base  4  cm .If  a right  circular  cylinder circumscribes the

Continuity requirement, Continuity requirement : Let's discuss the continu... Continuity requirement : Let's discuss the continuity requirement a little. Nowhere in the above description did the continuity requirement clearly come into play.  We need that t

Evaluate negative infinity, Evaluate both of the following limits. ... Evaluate both of the following limits. Solution : Firstly, the only difference among these two is that one is going to +ve infinity and the other is going to negative inf

Dynamical system and differential equations, 1. Discuss lyapunov function t... 1. Discuss lyapunov function theory and how it can be used to prove global assmptotic stability of solutions.(Give an example form natural and engineering sciences.) --- Draw le

Dividing by 2 digits, 3 divided bye 24 3 divided bye 24