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

Vectors, |a.x|=1 where x = i-2j+2k then calculate a

|a.x|=1 where x = i-2j+2k then calculate a

Arthemetic progreession, ball are arranged in rows to form an equilateral t...

ball are arranged in rows to form an equilateral triangle .the firs row consists of one abll,the second of two balls,and so on.If 669 more balls are added,then all the balls canbe

The parallelogram, love is a parallelogram where prove that is a rectangle...

love is a parallelogram where prove that is a rectangle

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

Positive skewness-measure of central tendency, Positive Skewness - It ...

Positive Skewness - It is the tendency of a described frequency curve leaning towards the left. In a positively skewed distribution, the long tail extended to the right. In

Ratios, the ratio of boys to girls in the sixth grade is 2:3 if there are ...

the ratio of boys to girls in the sixth grade is 2:3 if there are 24 boys, how many are girls?

Principle of superposition, If y 1 (t) and y 2 (t) are two solutions to a...

If y 1 (t) and y 2 (t) are two solutions to a linear, homogeneous differential equation thus it is y (t ) = c 1 y 1 (t ) + c 2 y 2 (t )   ........................(3) Remem

Weighted mean-progression, Weighted mean - It is the mean which employ...

Weighted mean - It is the mean which employs arbitrarily given weights - This is a useful measure especially whereas assessment is being done yet the situation prevailing a

Geometry, triangular with base AB = 48cm and height CH=16cm is inscribed a ...

triangular with base AB = 48cm and height CH=16cm is inscribed a rectangle MNPQ in which MN: MQ = 9:5 Find MN and MQ

Prepare a bar diagram, Question Write a short note on the following: ...

Question Write a short note on the following: 1 The weekly salaries of a group of employees are given in the following table. Find the mean and standard deviation of the

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