Discussion :

[herkitz]

Bonjour,

Je souhaite calculer l'intégrale (sur un triangle) suivante :

où la fonction phi est définie par :

J'utilise pour ce faire le code suivant :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
 
Phi=function(s,d){
	Phi= (exp(d)-exp(s))/(exp(d)-1)
	return(Phi)
}
Integrande2=function(s,d){
	 Res=max(0,(Phi(s,d)-1/s^5)^2)
	 return(Res)
}
INtegrale=function(){ integrate(function(y) { sapply(y, function(y) { integrate(function(x) Integrande2(x,y), 0, y)$value   }) }, 0, 200)}
INtegrale()

J'ai cependant cette erreur:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
 
Erreur dans integrate(function(x) Integrande2(x, y), 0, y) : 
  evaluation of function gave a result of wrong length

Auriez vous une idée?

Merci d'avance pour votre aide.