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

Range of f(x) =4^x+2^x+1 is, Taking 2^x=m and solving the quadratic for get...

Taking 2^x=m and solving the quadratic for getting D>=0 we get range= [3/4 , infinity )

Solve the inequality |x - 1| + |x - 2|, Solve the inequality |x - 1| + |x -...

Solve the inequality |x - 1| + |x - 2|≤ 3. Working Rule:    First of all measure the expression to zero whose modulus happens in the given inequation and from this search the va

Explain multiples, Explain Multiples ? When a whole number is multiplie...

Explain Multiples ? When a whole number is multiplied by another whole number, the results you get are multiples of the whole numbers. For example,  To find the first four mult

What is a system of equations?, What is a System of Equations? And its Solu...

What is a System of Equations? And its Solution? Here is an example of a system of equations (also called a simultaneous system of equations) x 2 + y = 3

Pair of straight lines, how to solve the problems? methods to solve the que...

how to solve the problems? methods to solve the question of joint lines

Find interval of function, Find interval for which the function f(x)=xe x(1...

Find interval for which the function f(x)=xe x(1-x)   is increasing or decreasing function

Matrices, how to find inverse of matrix

how to find inverse of matrix

A fire in a building b is reported on telephone, A fire in a building B is ...

A fire in a building B is reported on telephone to two fire stations P and Q, 10km apart from each other on a straight road.  P observes that the fire is at an angle of 60 o to th

How many points did he score during his senior year, Michael scored 260 poi...

Michael scored 260 points during his junior year on the school basketball team. He scored 20% more points during his senior year. How many points did he score during his senior yea

Logarithm, I need help with one logarithm problem

I need help with one logarithm problem

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