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

SURFACE AREA AND VOLUMES, Metallic spheres of radii 6 centimetre, 8 centime...

Metallic spheres of radii 6 centimetre, 8 centimetre and 10 centimetres respectively are melted to form a single solid sphere. Find the radius of the resulting sphere.

What are the three sides of a right triangle, What are the Three Sides of a...

What are the Three Sides of a Right Triangle? Each side of a right triangle can be labeled opposite, adjacent, or hypotenuse, based on its relationship to the right angle and o

Tangents, two circle of radius of 2cm &3cm &diameter of 8cm dram common tan...

two circle of radius of 2cm &3cm &diameter of 8cm dram common tangent

ALgebra, Please quote me a price

Please quote me a price

Trigonometric approximation grid, With a compass draw the arc associated wi...

With a compass draw the arc associated with a 720° angle, it looks like a circle. With a protractor, label the angle in multiples of 45° and 30° up to 720°.  Notice 30° and 390° ar

Hasse diagram, The digraph D for a relation R on V = {1, 2, 3, 4} is shown ...

The digraph D for a relation R on V = {1, 2, 3, 4} is shown below (a) show that (V,R) is a poset. (b) Draw its Hasse diagram. (c) Give a total order that have R.

Factors or multiples, long ago, people decided to divide the day into units...

long ago, people decided to divide the day into units called hours. they choose 24 as the number of hours in one day. why is 24 a more convenient choice than 23 or 25?

Linear programming problem., Ask question #Minimum 100 words acca paper mil...

Ask question #Minimum 100 words acca paper mill produces two grades of paper viz.,xand y.Bacause of raw material restrictions, it cannot produce more than 400 tones of grade x pape

Calculate the height of the tunnel and the perimeter, The adjoining figure...

The adjoining figure shows the cross-section of a railway tunnel. The radius of the tunnel is 3.5m (i.e., OA=3.5m) and ∠AOB=90 o . Calculate : i.       the height of the

Area and perimeter, a garden is constructed with a 3ft patio all around how...

a garden is constructed with a 3ft patio all around how would you give the expression for the area of the garden, excluding the patio

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