Answered
2D information to 1D colorplot
Simple example x = 0:10; y = x.^2; temp = rand(1,10); % temperature cm = jet(10); % colormap ind = 1+rou...

6 years ago | 0

Answered
Quiver with a system of ODEs.
For vector field u = x'; v = y'; To plot the curve [T,Y] = ode45 ... plot(Y(:,1),Y(:,2)) % (x,y) I suppose x and y depen...

6 years ago | 0

Answered
Solving nonLinear equation and Linearized Eqaution using ODE45
YOu have to express , and . Solve the system function du = nonlin(t,y) fi = y(1); psi = y(2); theta = y(3); ...

6 years ago | 0

| accepted

Answered
How to plot magnitude and directions between 2 points?
try this (not tested) data = load ('Input.txt'); [row, col] = size (data); Y2 = data (:,1); X2 = data (:,2); Y1 = data (:,3...

6 years ago | 0

| accepted

Answered
surface fixed point value Z
I want to know the value of Z from the surfc function Here is it

6 years ago | 0

Answered
Where are the bugs for this ODE finite difference problem that solve using Newton Raphson method?
Here is my attempt clc,clear N = 10; % 10 Set parameters L = 16*0.3048; %meter b = 19.625 * 0.0254; d = 1.625 * 0.0254; ...

6 years ago | 1

| accepted

Answered
I am receiving the error message Index exceeds the number of array elements (1) when trying to use ode45.
Try @(V,T)

6 years ago | 0

| accepted

Answered
I have a 14 equation system of second order differential equations which I am trying to solve using dsolve. The program just keeps running without finishing ?
alpha and E should be numerical function main m = 15000; k = 150000; n = 14; A = zeros(n); A(n,n) = -k/m; A(n,n-1) = k/m;...

6 years ago | 0

Answered
Euler 3D rotation between two vectors
Here is rotation using Rodriguez matrix Here is rotation about Y axis and Z. Here is the difference between them See ...

6 years ago | 0

| accepted

Answered
Simultaneous differential equations - derivative of y^3 wrt t
Just get dydt(1)=(-3*D*Cs/rho/r0^2)*y(1)*(1-y(2)); dydt(1)=nthroots(dydt(1),3);

6 years ago | 0

Answered
Numerically derive a continous, non-symbolic function
Derivative is (if you have numerical data) dy = (y(i)-y(i-1)) / (t(i)-t(i-1)); Maybe diff (if you have a function) ? syms x ...

6 years ago | 0

| accepted

Answered
3D Point Cloud - Gaussian Curvature
one way clc,clear % generate some data r = 3; t = linspace(0,10)'; x = r*cos(t); y = r*sin(t); z = sin(1*t); t1 = ...

6 years ago | 0

| accepted

Answered
Runge Kutta for system of eqs
I wrote function in the way similar to ode45 y_dot = @(x,dx)[dx; -4*dx-5*x]; %insert function to be solved for i = 1:n k1...

6 years ago | 0

Answered
The problem of finding intersection point of triangle mesh
To check every edge: for j = 1:size(face,2) i1 = 1; for i2 = [2 3 1] % checking edges [1 2], [2 3], [3 1...

6 years ago | 1

| accepted

Answered
How can i solve a system which contains matrices
Here is the best way i know: A = [5 3 1 1 1 1 3 1 2]; B = [c1(:) c2(:) c3(:)]'; C = inv(A) * B;

6 years ago | 1

| accepted

Answered
Using a for-loop to score values in different intervals
Correct form if 0 <= x(i) && x(i) < 1

6 years ago | 0

| accepted

Answered
Intersection between line and circle when line ends inside circle (using geom2d)
Someone should copy their comment to the answers section so it can be accepted. It should be me % Position of line: x1 = [0...

6 years ago | 1

| accepted

Answered
Removing points from meshgrid (mask)
Try this simple script [X,Y] = meshgrid(1:20); n = [1 -1]; n = n/norm(n); % normal of a line dist = -7; ...

6 years ago | 0

Answered
How can I find the surface normal of a 8 by 3 matrix?
Use griddata to create a surface Use surfnorm to find normal vectors [nx,ny,nz] = surfnorm(x,y,z); surf(x,y,z) hold on quiv...

6 years ago | 0

Answered
How do I increase the circle size and animate a stick figure?
Use for loop and pause t = linspace(0,2*pi); x = cos(t); y = sin(t); axis([-10 10 -10 10]) hold on for i = 1:10 h = p...

6 years ago | 0

| accepted

Answered
Finding specific data point x for given y
I need to find the CB and CS values at t=5 Use interp1 to find i need to find the t and CB values when CS = 0.001 Use fsolve ...

6 years ago | 0

Answered
Surface between 3D lines
You are concatenating in a wrong way XX = [D1 D1 D1]; YY = [F1 F2 F3]; ZZ = [z1 z2 z3]; surf(XX,ZZ,YY,'EdgeColor','none')

6 years ago | 2

| accepted

Answered
two body problem using ode45
- equation of equilibrium Is this condition fullfield?

6 years ago | 0

Answered
Numerically Integrating the differential equation below?
Any guidance on it is much appreciated! Helpfull page: ODE45 I do not know what I can do if I no tf. Try tf = 5. This should ...

6 years ago | 0

Answered
Streamline plot not complete
Works ok for me N = 15; xx = linspace(min(x(:)),max(x(:)),10); yy = linspace(min(y(:)),max(y(:)),10); [startx,starty] = mes...

6 years ago | 1

| accepted

Answered
2D Matrix revolution
Use cylinder for each curve separately: z = linspace(-1,1,20); % axis r = sqrt(1-z.^2); ...

6 years ago | 1

Answered
How do I align the vertices and the adjacency graph of a Buckyball?
Here is a way ind = full(B); % convert to ordinary matrix [rr,cc] = find(ind); % find non-zero elements ix = ...

6 years ago | 0

Answered
Find the coordinates of a point chosen on a plot
Short example x = linspace(0,2,20)'; % generate some data y = sin(x); plot(x,y,'.-b') xy = ginput(1); ...

6 years ago | 3

| accepted

Answered
How do I create a sphere without sphere function that will inflate like a balloon utilizing the getframe function?
Use surf to create a 3D surface rr = linspace(0,0.5,20); tt = linspace(0,2*pi,40); [R,T] = meshgrid(rr,tt); [X,Y] = pol2cart...

6 years ago | 0

Answered
Differential Algebraic Equations with vectors
Try f = @(t,u) inv(Z)*(J+p*W*pi*Q)*[u(3) u(4)]' + inv(Z)*(V+p*W*pi*Q)*[u(1) u(2)]'; F = @(t,u) [u(3) u(4) ...

6 years ago | 1

| accepted