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

Which mathematical property did marty use to get similar ans, Marty used th...

Marty used the subsequent mathematical statement to show he could change an expression and still get the similar answer on both sides: 10 × (6 × 5) = (10 × 6) × 5 Which mathematica

Binomial distribution, Binomial Distribution Consider a batch of N ligh...

Binomial Distribution Consider a batch of N light bulbs. Each bulb may be defective (S) or non-defective (F). The experiment involves selecting a light bulb and checking whethe

Iceviga, Jess had a book with 100 pages to read she only read 10 how many p...

Jess had a book with 100 pages to read she only read 10 how many pages does she have to read?

Index numbers, advantages and disadvantages of index numbers

advantages and disadvantages of index numbers

Divide by 1-digit numbers, which experession can be used to check the quoti...

which experession can be used to check the quotient 646 divided by 3

Explain basic concepts of parallel lines, Explain Basic Concepts of Paralle...

Explain Basic Concepts of Parallel Lines ? Parallel lines are defined in section 1.2 and we use "//" to denote it. From the definition, we can get the following two consequenc

Tchebecheffs ineqality theorom, what are the advantages and disadvantages o...

what are the advantages and disadvantages of tchebycheffs inequality theorem

Decimals, how to multiply 8654.36*59

how to multiply 8654.36*59

the height of the tower, A Stone is dropped from the top of the tower and ...

A Stone is dropped from the top of the tower and travel 24.5 m in last second of its journey. the height of the tower is ...?

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