Mathematica Mathematica Looping Help: Find 'v' for 10 Different Values of 't

  • Thread starter Thread starter suraj goyal
  • Start date Start date
  • Tags Tags
    Mathematica
AI Thread Summary
The discussion focuses on a Mathematica code issue where the user seeks to calculate and plot the variable 'v' for ten different values of 't'. The user initially fixed 't' at the beginning of the program, which led to difficulties in generating multiple outputs. Suggestions include defining all variables as functions of 't' and using the Table function instead of a For loop to store results. The user is advised to ensure that values are stored in a variable for plotting, rather than just printed. Ultimately, the goal is to successfully create a graph that represents the relationship between 't' and 'v'.
suraj goyal
Messages
2
Reaction score
0
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...
 
Physics news on Phys.org
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])
...
 
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:
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.
 

Similar threads

Replies
3
Views
2K
Replies
4
Views
2K
Replies
6
Views
10K
Replies
1
Views
3K
Replies
15
Views
4K
Replies
52
Views
12K
Replies
2
Views
4K
Back
Top