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

Minima, Minima, Maxima and points of inflexion a)      Test for rela...

Minima, Maxima and points of inflexion a)      Test for relative maximum Consider the given function of x whose graph is presented by the figure given below

Tent originally sold for $2 what is the percent of discount, A tent origina...

A tent originally sold for $260 and has been marked down to $208. What is the percent of discount? Find out the number of dollars off. $260 - $208 = $52. Further, determine wha

Prove - digraph of a partial order has no cycle more than 1, Prove that the...

Prove that the Digraph of a partial order has no cycle of length greater than 1. Assume that there exists a cycle of length n ≥ 2 in the digraph of a partial order ≤ on a set A

Show trigonometric functions on a graph, Q. Show Trigonometric Functions on...

Q. Show Trigonometric Functions on a Graph? Ans. By discussing the trig functions with respect to an angle in a right-angle triangle, we have only considered angles betwee

What is the meaning of the f-statistic and f test, You are given the follow...

You are given the following regression results estimating the demand for widgets based on time series data for the past 40 months. Q t = 2.5 - 0.3 x P t + 12 x M t Where Q

Farmer counting grasshoppers in his fields, Farmer counting grasshoppers in...

Farmer counting grasshoppers in his fields, probably not normally distributed due to growing conditions. After various rows the mean number of grasshoppers is 57 SD 12. What will b

Find third order partial derivatives, Question: Find all third order pa...

Question: Find all third order partial derivatives for the function   F(x,y)= log xy+ e (x+y) -x/y.

Area and perimeter, if perimeter is 300m length is 100m.find the breadth

if perimeter is 300m length is 100m.find the breadth

Triangles, In a triangle ABC, D &E is a are points on AB & AC ,if the one s...

In a triangle ABC, D &E is a are points on AB & AC ,if the one side of a triangle is 4cm & another side is 5 cm find that the ar(triangleABC):ar(BCDE)

Share and dividend, i want to get market value of 10 popular shares of all ...

i want to get market value of 10 popular shares of all working days in a week

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