How to Integrate exp(x)*erfc(x) in MATLAB?

[bhartish]

I have doubt in MATLAB

How to integrate

exp(x)*erfc(x) ? in MATLAB

The warning msg it is displaying is :

EXPLICIT integral could not be found

 

[Unrest]

Are you trying to do numerical or symbolic integration? Maybe that's a function which has no analytic integral? If numerical, what are the limits?

 

[bhartish]

I have an expression involving seven parameters. I need to first integrate the entire expression from 0 to t and then get the values by substitution method. Matlab is not supporting symbolic integration. hence I am planning to Numerical .

 

[Unrest]

Oops I didn't see the title of you message. In that case there must be something wrong with your Matlab code since it should be able to integrate anything it can evaluate. Could you paste the code, including the limits of the integration?

Seems like a function that can easily be numerically integrated over any range of x.

Does it actually produce a solution, but you don't trust it because of the warning message?

 

[bhartish]

e^(a+bt) erfc(-c+rt)exp((t-q)/l)

within the limits o to t

 

[MATLABdude]

If you're using undefined a, b, c, r, q, and l, that's not numerical integration (especially if you're also attempting to integrate to an undefined 't').

 

[bhartish]

I have all the values of a, b, c, r, q, l and t. Then how to go about it ?

 

[MATLABdude]

Write a small script to do so. For instance, the following uses the trapezoid rule to find the integral between 0 and pi of the sine function.

	t=0:pi/100:pi;	% t values
	y=sin(1*t);		% calculation of function at all t
	A = trapz(y,t);	% implements trapezoid rule

Why did I explicitly use 1*t? While I could have just used t, it's to illustrate that there is no implicit multiplication in MATLAB (i.e. rt refers to the value / vector / matrix rt, while r*t refers to the product of r and t). I'd suggest writing up a small function (just as I've done above) in Notepad / Text Editor / EMacs / Vi, and then paste it here between the [code ][/ code] tags (no spaces in the tags).

 

[bhartish]

I will get back to you sir.

 

[bhartish]

I will clear one doubt. rt is r*t only. rt is definitely not a single variable.

I have one doubt can't we take some fictitous value and start integration of error function and exponential ?

 

[MATLABdude]

I will clear one doubt. rt is r*t only. rt is definitely not a single variable.

I have one doubt can't we take some fictitous value and start integration of error function and exponential ?

It's your work--you can put in whatever you want!