Hopefully, anyone can help me with the following problem:
I have two functions: Inta=f(Pa) and ETa=f(Pna). Both can not be solved analytically, but I can plot them separately when I use >plot (evalf(Inta), Pa=0..100) and >plot(evalf(ETa), Pna=0..100). So far no problem.
Furthermore, Pna can be described as >Pna:=Pa-Inta, so in fact ETa=f(Pa). But I fail to plot ETa against Pa (or at least it takes ages..since it keeps 'evaluating'..... for days) . I tried several things like:
But I don't get a plot, only 'evaluating....'.
The equations for Inta and ETa are:
Inta := Pa*Int(Pm*(1-exp(-1.220/(Pm^.45-.2*Pm^.69+0.2e-1*Pm)))*exp(-Pm/(0.97e-1*Pa+8.57))/(0.97e-1*Pa+8.57), Pm = 0 .. infinity)/(0.97e-1*Pa+8.57)
ETa:= 9.285051068*Int(9.285051068*min(.5676676416*Pnm, 50.)*exp(-9.285051068*Pnm/Pna)/Pna, Pnm = 0. .. Float(infinity))
Does anybody know how I can do this?
Thanks in advance,