I have used data and the curve fitting procedure to find a curve and I now have to put error bars on the graph to complete my task? I know that there is an errorplot function but I cant seem to find a way to do it with the curvefitting thing?

I have attached the whole worksheet, I hope someone can help! Thanks Detectors_Tutorial.mwDetectors_Tutorial.mw

Even if there was some method of exporting the graph to another medium that would also be very helpful

Hello, 

I am pretty new to Maple and Im trying to find the parameters of an equation using some kind of fit routine but I can only find such a routine to fit an expression to data and not vice versa.

my equation is as follows

epsilon(E):= a0 + a1*ln(E) + a2*ln(E)^2 + a3*ln(E)^3

I do have data for E which I imagine I need. 

E:= {121, 244, 344, 411, 444, 778, 867, 964, 1085, 1112, 1212, 1299, 1408} all in keV :)

Any suggestions/help would be much appreciated, although I am new to the program so go easy on me :) 

Cheers

Hi,

I am trying to find the roots of Hankel function H1(2, z)?

j := 1;
for i from 0 to 10 do z0 := i*step+z_min; x[j] := fsolve(f = 0, z = z0, complex) end do;

for p to 10 do for j from p+1 to 9 do if `and`(Re(x[j])-Re(x[p]) < 0.1e-4, Im(x[j])-Im(x[p]) < 0.1e-4) then for i from j to 10 do x[i] := x[i+1] end do; p := p-1; break end if end do end do;

the first of these two processes works fine however the second does not. The second on is to get rid of same value solutions! I am not sure if I have missed anything and also is there a way to determine a max value in the complex domain and use it in a for loop?

any help would be great 

I'm have used a program to find the roots of a function 

f:=x*cos(x)-sin(x)*sin(x/1000);
/ 1 \
x cos(x) - sin(x) sin|---- x|
\1000 /

x_max:=50; x_min:=-50; step:=2; i_max:=(x_max-x_min)/step;
50
-50
2
50

j:=1:
for i from 0 to i_max by 1 do
x0:=x_min+i*step;
x[j]:=fsolve(f=0,x=x0);
j:=j+1;
end:

and my output was of the form of multiple "potential" roots and a bunch of which are the same. So I tried to get rid of the ones which were the same before actually finding the ones which ARE roots. To do that I done....

j := 1; for j to 50 do if x[j]-x[j+1] = 0 then ignore(x[j]) else print(x[j]) end if end do:

and it got rid of the ones which are of the above form but some roots are the same and seperated by more than 1 ... i.e x[ j ]= x[j + 2] or some other number. 

Basically I am trying to generalise the above for loop for all "numbers" instead of 1 but when I try some things the for loop doesnt like it. 

Any help would be good!