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

Geometry, two sides of an equilateral triangle have lengths 3x-1 and 3x-1. ...

two sides of an equilateral triangle have lengths 3x-1 and 3x-1. Which of 27-x or 2x-4 could be the length of the third side?

Show that 571 is a prime number, Show that 571 is a prime number. Ans: ...

Show that 571 is a prime number. Ans:    Let x=571⇒√x=√571 Now 571 lies between the perfect squares of  (23)2 and (24)2 Prime numbers less than 24 are 2,3,5,7,11,13,17,1

Properties of summation notation, Properties Now there are a couple of ...

Properties Now there are a couple of formulas for summation notation. 1. here c is any number. Therefore, we can factor constants out of a summation. 2. T

Angles, samuel left mauritius at 22:30 on saturday and travelled to london ...

samuel left mauritius at 22:30 on saturday and travelled to london (GMT) for 14h30min he had a stopover for 4 h in london and he continued to travel to toronto for another 6h20min

How many 6-inch tiles are required to tile the floor, Mark intends to tile ...

Mark intends to tile a kitchen floor, which is 9 by 11 ft. How many 6-inch tiles are required to tile the floor? a. 60 b. 99 c. 396 c. Since the tiles are calculated in

Determine the determinant of matrix, Example Determinant:   Determine ...

Example Determinant:   Determine the determinant of each of the following matrices. Solution : For the 2 x 2 there isn't much to perform other than to plug this in

Solve sin (a /7) =0 trig function, Solve sin (α /7) =0 . Solution B...

Solve sin (α /7) =0 . Solution By Using a unit circle it isn't too difficult to see that the solutions to this equation are, α /7 = 0 + 2 ? n     ⇒   α = 14 ? n

difference between two sample means (large sample), Testing The Difference...

Testing The Difference Between Two Sample Means (Large Samples) A large sample is defined as one which have 30 or more items as n≥30 whereas n is the sample size In a busine

Circles, assignment on theorems on circle for class 9

assignment on theorems on circle for class 9

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