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

Measurement story problem, Seth has a pet goldfish. When he got his goldfis...

Seth has a pet goldfish. When he got his goldfish , it was only 5 centimeters long . Now it has grown to be 92 millimeters long. How many millimeters has the goldfish grown since

Truth criteria-nature of mathematics, Truth Criteria :  Consider the follo...

Truth Criteria :  Consider the following statements: i) Peahens (i.e., female peacocks) lay eggs around September. ii) Water boils at 100°C. iii) 5 divides 15 without lea

Variation and proportion, i am not getting what miss has taught us please w...

i am not getting what miss has taught us please will you will help me in my studies

Prove gcd value, Let a, b, c 2 Z + . (a) Prove that if a|b, then ac|bc f...

Let a, b, c 2 Z + . (a) Prove that if a|b, then ac|bc for all c. (b) If a|bc, can you conclude that either a|b or a|c? Justify your answer with a proof or a counter example.

Example of learning to count, A parent shows his child four pencils. He pla...

A parent shows his child four pencils. He places them in a row in front of her and says "one" as he points to the first pencil, "two" as he points to the second one, "three" as he

Magnitude - vector, Magnitude - Vector The magnitude, or length, of th...

Magnitude - Vector The magnitude, or length, of the vector v → = (a1, a2, a3) is given by, ||v → || = √(a 1 2 + a 2 2 + a 2 3 ) Example of Magnitude Illus

Trigonometry, explain the formular for finding trigonometry

explain the formular for finding trigonometry

Factors, Question Suppose that f(x) has (x - 2) 2 and (x + 1) as its on...

Question Suppose that f(x) has (x - 2) 2 and (x + 1) as its only factors. Sketch the graph of f. State all the zeros of f.

Partial fractions - integration techniques, Partial Fractions - Integration...

Partial Fractions - Integration techniques In this part we are going to take a look at integrals of rational expressions of polynomials and again let's start this section out w

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