Initial value problems, math, Mathematics

Assignment Help:
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.

Related Discussions:- Initial value problems, math

Find the instantaneous rate, The time t required to test a computer memor...

The time t required to test a computer memory unit is directly proportional to the square of the number n of memory cells in the unit. For a particular type of unit, n = 6400

Probability of chosen number from 1st 500 divisble by 3or5 , IN THIS WE HAV...

IN THIS WE HAVE TO ADD THE PROBABILITY of 3 and 5  occuring separtely and subtract prob. of 3 and 5 occuring together therefore p=(166+100-33)/500=233/500=0.466

Solid mensuration., assuming that the earth''s sphere with a radius of 6400...

assuming that the earth''s sphere with a radius of 6400 km.. find the distance along a 3 degree arc at the equator of the earth''s surface?

Derivative for the trig function, Derivative for the trig function: We'll ...

Derivative for the trig function: We'll begin with finding the derivative of the sine function. To do this we will have to utilize the definition of the derivative. It's been wher

Travel time, you are driving on a freeway to a tour that is 500 kilometers ...

you are driving on a freeway to a tour that is 500 kilometers from your home. after 30 minutes , you pass a freeway exit that you know is 50 kilometer from your home. assuming that

What is pythagorean triples, What is Pythagorean Triples? A set of thre...

What is Pythagorean Triples? A set of three numbers a, b, and c that can satisfy the equation A 2 +b 2 = c 2 , is called a Pythagorean triple. The following is a list of

Proof of root test - sequences and series, Proof of Root Test  Firstly...

Proof of Root Test  Firstly note that we can suppose without loss of generality that the series will initiate at n = 1 as we've done for all our series test proofs.  As well n

Solving decimal equations, The distance around a square photograph is 12.8 ...

The distance around a square photograph is 12.8 centimeters. What is the langth of each side of the fotograph?

Write Your Message!

Captcha
Free Assignment Quote

Assured A++ Grade

Get guaranteed satisfaction & time on delivery in every assignment order you paid with us! We ensure premium quality solution document along with free turntin report!

All rights reserved! Copyrights ©2019-2020 ExpertsMind IT Educational Pvt Ltd