Initial value problems, math, Mathematics

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







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
Related Questions

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 ∼ Δ


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

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.

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 : 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 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

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