Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Mathematica looping

  1. May 26, 2012 #1
    pl help
    Here I have written a mathematica code. It is a small example of my actual problem. In this code, 't' is an arbitrary constant which is being fixed at the top of the programme by assigning a value. All the other variables 'a','b','c' depend upon 't' . and 'h1','h2' are the roots of quadratic equation with coefficients a,b,c.
    Further, 'f' is a function in an unknown variable 'p' and involves 'h1' and 'h2'.. We can determine 'p' by solving f(p)=0, giving two roots. Finally I calculated 'v' from the from assigned 't' and calculated second root of equation f(p)=0.
    Thus 'v' depends upon 't' both explicitly and implicitly.. With change in 't', "v" will change. I want to draw a graph between 't' and 'v'.. For this I need looping or some other command to assign different values to 't' at the top of the programme.. Me unable to do that.. Pl if possible , help me in this regard
    t := 3
    a := t^2
    b := (t + 1)*I
    c := t^2 + 2*t*I + 1
    h1 = (-b + Sqrt[b^2 - 4*a*c])/(2*a)
    h2 = (-b - Sqrt[b^2 - 4*a*c])/(2*a)
    a11 := Sqrt[p]*h1
    a12 := h1*2
    a21 := 3*h2^2
    a22 := 4 + p
    f[p_] = a11*a22 - a21*a12
    h = sols = NSolve[f[p] == 0, p]
    d = p /. sols[[2]]
    v = t/Re[d]
    Suppose I want to calculate my final answer for ten values of 't' (say 1,2,3.....10) and finally to plot a graph between the final values of 'v' and the ten assigned values for 't'. I have tried a lot in this regard and have defined a function, used array and table commands too.. But unable to pick those ten values by my for loop command and to draw graph...
     
  2. jcsd
  3. May 26, 2012 #2

    Dale

    Staff: Mentor

    You should not set t at the beginning. It seems that t is the argument to a function, so you should express everything in terms of functions of t:

    a[t_] := t^2
    b[t_] := (t + 1)*I
    c[t_] := t^2 + 2*t*I + 1
    h1[t_] := (-b[t] + Sqrt[b[t]^2 - 4*a[t]*c[t]])/(2*a[t])
    ...
     
  4. May 26, 2012 #3
    Now I have tried the following

    For[t = 1, t < 100, t = t + 10,
    a = t^2;
    b = (t + 1)*I;
    c = t^2 + 2*t*I + 1;
    h1 = (-b + Sqrt[b^2 - 4*a*c])/(2*a);
    h2 = (-b - Sqrt[b^2 - 4*a*c])/(2*a);
    a11 = Sqrt[p]*h1;
    a12 = h1*2;
    a21 = 3*h2^2;
    a22 = 4 + p;
    S = {{a11, a22}, {a21, a22}};
    f[p_] = Det;
    h = sols = NSolve[f[p] == 0, p];
    d = p /. sols[[2]];
    v = t/Re[d] ;
    data = {{v, t}} // TableForm;
    Print[data];
    ]
    The output is
    {{2.41743, 1}}
    {{-2.75, 11}}
    {{-5.25, 21}}
    {{-7.75, 31}}
    {{-10.25, 41}}
    {{-12.75, 51}}
    {{-15.25, 61}}
    {{-17.75, 71}}
    {{-20.25, 81}}
    {{-22.75, 91}}
    But problem is to plot.. In plotting command, it willl take the final value i.e. -22.75 and 91.. So still having the problem.. Want to plot a graph between 't' and 'v............'Pl help
     
    Last edited: May 26, 2012
  5. May 26, 2012 #4

    Dale

    Staff: Mentor

    Don't do it that way. Make everything into functions of t. Avoid explicit looping structures in interpreted languages whenever possible.

    The reason you are unable to plot the data is because you never set anything equal to the data, you simply printed it. Printing it displays it but does not store it in a variable.

    Instead of using For you should use Table to actually build up a table of values.
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook




Similar Discussions: Mathematica looping
Loading...