## 120 Reputation

12 years, 110 days

## Almost there!...

Thanks again Rouben.

The last couple lines, the value of C_2 varies with each input. in your code:

c_2(my_x, my_xbar,1) is 0.00069

c_2(my_x, my_xbar,9) is 0.00076

c_2(my_x, my_xbar,17) is 0.0007

I am looking the sum of these numbers from 1 to length of my input

Thanks a lot for helping me

## Xbar and x[i] are independent...

Thanks Rouben. Please consider Xbar and x[i] as independent variables. like x and y.

Xbar is obtained from a different series of measurments.

My question boils down to:

How can I subtitute x in c1 and c2.

Thanks

## No need to insult dude !...

English is my second language and Im learning it. I thought its clear! you can simply say you dont want to/cant/... help. These are some lines from a larger code and some of these "crap" will be used later.

Anyways thank you for your time :)

## still need help!...

I dont know where is the problem when I change your code. I have attached the modification. I simply multiply the integral by a function but cannot get the results.

integral.mw

## @Markiyan Hirnyk    Thanks m...

Thanks man. you saved me.

## @Preben Alsholm  I ran the code fo...

I ran the code for different functions and plotted the contour plots. Coppied and pasted in a new worksheet and then the next function and the next contour plots.

Doing this for years and first time facing this issue. Very disappointed Maple! Very!

Thank you for your time tho. Appreciate the time you put to help me.

## @Preben Alsholm  btv_3_contour_plo...

btv_3_contour_plots.mw

418_3_contour_plots.mw

## @Carl Love  many thanks for yo...

many thanks for your attention Carl

## @Kitonum  I dont know why it doesnt...

I dont know why it doesnt work for me! here is my function

PP:=0.3800179925e-3*exp(-0.6065722618e-3*(x-29.51704536)^2+(0.6650093594e-3*(x-29.51704536))*(a-12.94061928)-0.1106850312e-2*(a-12.94061928)^2)

I simply copied and pasted the procedure to a maple sheet and substitute x^2+y^2 with PP! Is it a wrong way? help me please

## @tomleslie  Thats why I love phoros...

Thats why I love phoroshop! lol

Thanks any way :)

Awesome. Problem solved.
Thanks a lot :)

## @Kitonum    Thanks for your re...

Thanks for your reply. Im taking the 3d graph. change rho and see what happens to the 3d graph. Maybe the coefficient of (x-constant)*(a-constant) gets large ( in the exponential) abd maple fails to plot. I have attached two pictures to show what I mean.

when rho increases, maple plots several pics!

thank you

## changing the colors...

How can I change the colors?

The colors appear are not different enough, it may be confusing. I need more distinct colors. I changed c1 to

c1:=[Color("#123456")[]]: still in a same color family. I am seeking for more obvious change in colors,

thanks again

## @Christopher2222    Thanks Chr...

I cant get the results. I would be grateful if you help me to know what is the problem. I despretly need some help from a pro. I have uploaded the file.

Thanks

heartview_-_Copy.mw

 >
 >
 (1)
 >
 (2)
 >
 (3)
 >
 (4)
 >
 (5)
 >
 (6)
 >
 (7)
 >
 >
 (8)
 >
 (9)
 >
 >
 > plot3d(PPP, x = 0 .. 200, y = 0 .. 110, labels=["x(m)","a(cm^2)","PDF(x,a)"]);
 >
 >
 >
 >
 >
 >
 >

## The code...

Thanks again and again ...

restart;

n:=20;
m:=1;
cc:=-200;
zzeta:=0.1;
sefr1:=0.3;
sefr:=0.2;
MM:=0;
lambda:=0.1;
Br:=1;
nn:=3;

u(tau):=421.7129935*tau-2217.587728*tau^2+8897.376593*tau^3-27612.59182*tau^4+64248.00336*tau^5-1.083977605*10^5*tau^6-10.57029600-1.080951714*10^6*tau^13+7.999517316*10^5*tau^14-4.788741005*10^5*tau^15+2.309563748*10^5*tau^16+26511.11102*tau^18-5959.001794*tau^19+1.148523882*10^5*tau^7-95.23809524*tau^21+4.545454545*tau^22-9435.563781*tau^8-2.587683745*10^5*tau^9+6.473880128*10^5*tau^10+948.0272727*tau^20-88660.41892*tau^17-1.008692404*10^6*tau^11+1.175504242*10^6*tau^12;

B := 1+(2*(1-zzeta))*Br*(int(tau*(diff(u(tau), tau))^2, tau = (1/2)*zzeta*(1-zzeta) .. 1/2*(1-zzeta)));
eq4 := 4*B*u(tau)-(1+zzeta)*(diff(tau*(diff(theta(tau), tau)), tau))/tau-(1+zzeta)*Br*(diff(u(tau), tau))^2;

theta(tau):=sum(p^ii*theta[ii](tau),ii=0..nn);
HH:= p*((4*(1+(2*(1-zzeta))*Br*(int(tau*(diff(u(tau), tau))^2, tau = (1/2)*zzeta*(1-zzeta) .. 1/2-(1/2)*zzeta))))*u(tau)-(1+zzeta)*(diff(theta(tau), tau)+tau*(diff(theta(tau), tau, tau)))/tau-(1+zzeta)*Br*(diff(u(tau), tau))^2)+(1-p)*(diff(theta(tau),tau\$2)):
eq5:=simplify(HH):
eq6:=collect(expand(eq5),p);

eq7:=
convert(series(collect(expand(eq5), p), p, nn+1), 'polynom');

for ii to nn do
ss[ii] := (coeff(eq7, p^ii)) ;
print (ii);
end do;

ss[0]:=diff(theta[0](tau), tau, tau);

icss[0]:=theta[0](0)=0, D(theta[0])(1)=1;

dsolve({ss[0], icss[0]});
theta[0](tau):= rhs(%);

for ii to nn do
ss[ii]:=(ss[ii]);
icss[ii]:=theta[ii](0)=0, D(theta[ii])(1)=0;
dsolve({ss[ii], icss[ii]});
theta[ii](tau):=value(rhs(%));
end do;

theta(tau):=evalf[10](sum((theta[jj])(tau),jj=0..nn));
> plot(, tau = .1 .. 1);

 1 2 Page 1 of 2
﻿