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

Prove that its inclination theta to the horizontal, Two stations due south ...

Two stations due south of a tower, which leans towards north are at distances 'a' and 'b' from its foot. If α and β be the elevations of the top of the tower from the situation, Pr

Find the number of students in the class, Students are made to stand in row...

Students are made to stand in rows. If one student is extra in a row there would be 2 rows less. If one student is less in a row there would be 3 rows more. Find the number of stud

Relations, Suppose A and B be two non-empty sets then every subset of A Χ B...

Suppose A and B be two non-empty sets then every subset of A Χ B describes a relation from A to B and each relation from A to B is subset of AΧB. Normal 0 fals

Profit and loss, a shopkeeper buys two cameras at the same price . he sells...

a shopkeeper buys two cameras at the same price . he sells one camera at a profit of 18% and the other at a price of 10% less than the selling price of the first camera. find his p

Find the value of delta, Consider the given graph G below. Find δ( G )=__...

Consider the given graph G below. Find δ( G )=_____ , λ( G )= _____ , κ( G )= _____, number of edge-disjoint AF -paths=_____ , and number of vertex-disjoint AF -paths= ______

Vectors, A triangle has vertices A (-1, 3, 4) B (3, -1, 1) and C (5, 1, 1)....

A triangle has vertices A (-1, 3, 4) B (3, -1, 1) and C (5, 1, 1). The area of ABC is a) 30.1 b) 82.1 c) 9.1 d) 52.1

Inventory record, a) Complete the inventory record below for an FOQ of 100 ...

a) Complete the inventory record below for an FOQ of 100 units. b) Talk about weaknesses of MRP. List at least 3 and describe each in a sentence or two. Item: A

Draw a graph model with the adjacency matrix, QUESTION (a) Draw a graph...

QUESTION (a) Draw a graph model with the following adjacency matrix.                         (b) The diagram below shows different cities labelled a to g and z. Also sh

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