thank you very much :) it works and this is the result do you know what this symbol indicates in the result ??? ~ w := (1/4)*(-K*ln(y1~+s)+ln(s-y1)*K)/Pi

it should be floating point but i have problem with the previous calculations that make it symbolic not a number ,i am trying to fix it now

THANKS its not a floating point

i have anotehr problem i removed the space but i have another error now printf(cat(`$`(`%12.5f`, 4)), vij(1, 1)[k]*`$`*k = 1 .. 4) printf(cat(`$`("% 12.5f ", 4)),`$`(vij(1, 5)[k], k = 1 .. 4)) ther error is Error, (in fprintf) number expected for floating point format thanks

thanks

thanks but my problem is that when i want to zoon out my figure ,and when i choose unconstrained ,the directions of the plot doesn't appear like what i want. is there any solution to zoom the figure???????????

this is the solution of the 1st step,hope it helps u aa := Vector(21); lat := Vector(21); lon := Vector(21); for i to 21 do lat[i] := i-11; lon[i] := 3.*(i-11); aa[i] := 1.4*cos((1/180)*Pi*lat(i))-2*sin((1/180)*Pi*lon(i)); end do what do you want plot?? aa versus long and lat OR x,y versus long and lat ,or what ??

thank you very much it was really helpfull mohamed

hello mr david did you get any solutoin to ur probelm????? i face the same problem now.

thank you very much this is very helpful way, i solved this problem by drawing it in a for loop,but i think ur way is easier and faster for b from 1 to Nmx do P[b]:=plot([x|U(..,b)]); end do; display([P[j]$j:=1..Nmx])